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INVESTIGATION  OF  COMPLEX  ANGLE  PROCESSING 
TO  REDUCE  INOUCED  ANGLE  POINTING  ERRORS 

David  G.  Burks  and  Edward  R.  Graf 
Abstract 

This  work  is  an  analysis  of  the  effect  of  a  tangent  ogive  radome  on 
the  pointing  accuracy  of  a  monopulse  radar  that  employs  an  aperture  an¬ 
tenna.  The  radar  is  assumed  to  be  operating  in  the  receive  mode  and 
the  incident  fields  at  the  antenna  are  found  by  a  ray-tracing  procedure. 
Rays  that  enter  the  antenna  aperture  by  direct  transmission  through  the 
radome  and  by  single  reflection  from  the  radome  interior  are  considered. 
The  theory  of  monopulse  radar  and  the  transmission  and  reflection 
properties  of  planar  dielectric  slabs  are  presented  first  to  form  a  back¬ 
ground  for  the  radome  analysis. 

Two  orthogonal  polarization  states  which  can  be  combined  to  pro¬ 
duce  an  arbitrarily  polarized  incident  field  are  considered.  The  anten¬ 
na  can  be  scanned  in  two  angular  directions  and  radar  pointing  error  is 
presented  for  both  these  angles  as  a  function  of  antenna  scan  angle 
and  polarization  of  the  incident  field.  Throughout  the  work,  compari¬ 
son  of  two  different  radome  wall  designs  is  made  in  order  to  both  il¬ 
lustrate  the  analytical  techniques  and  to  show  the  engineering  trade¬ 
offs  in  radonrs  design. 

A  method  that  can  be  used  to  compensate  for  radome- induced  error 
Is  presented.  Complex  angle  processing  for  the  purpose  of  pointing 
angle  reduction  was  investigated  based  upon  the  simulation  described 


above.  It  was  concluded  that  the  complex  angle  could  not  be  uniquely 
related  to  radar  pointing  error,  but  rather  appeared  to  correlate 
with  phase-front  curvature.  This  quantity  Is  not  considered  useful  for 
pointing  error  correction. 
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I.  INTRODUCTION 


The  Institute  of  Electrical  and  Electronics  Engineers  gives 

the  following  definition  [1] 

RADOME.  An  enclosure  for  protecting  an  antenna  from 
the  harmful  effects  of  its  physical  environment,  gen¬ 
erally  intended  to  leave  the  electrical  performance 
of  the  antenna  unaffected. 

Skolnik  [2]  divides  radomes  into  these  two  major  classifications, 
based  on  application  and  geometrical  constraints:  ground-based 
radomes  and  air-borne  radomes.  This  study  deals  with  radomes  used 
on  high-speed  missiles  which  require  a  high  degree  of  streamlining. 
Streamlining  places  severe  constraints  on  the  geometry  of  the  radome 
and  the  radome  material,  which  must  withstand  the  temperature,  pres¬ 
sure,  and  possibly  rain  or  lightning  encountered  during  high-speed 
flight.  The  antenna  considered  here  is  part  of  a  radar  system 
that  provides  information  concerning  target  position  and  velocity 
to  the  guidance  system  of  the  missile.  The  radome  may  be  considered 
a  "necessary  evil"  because  it  has  no  beneficial  effects  on  the  per¬ 
formance  of  the  radar  but  it  is  essential  for  protection  of  the  antenna 
and  for  good  flight  characteristics  of  the  missile.  Of  particular 
interest  is  the  angular  pointing  error  of  the  radar.  This  study  is 
addressed  mainly  at  that  problem. 

Airborne  radomes  first  had  extensive  application  on  aircraft 
in  World  War  II.  Many  design  criteria  and  construction  requirements 
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were  formulated  during  this  period.  Much  of  this  work  In  the 
United  States  and  England  is  summarized  in  a  volume  of  the  MIT 
Radiation  Laboratory  Series  [3].  In  another  volume  [4]  of  the 
same  series,  radomes  are  discussed  In  conjunction  with  antenna  in¬ 
stallation  problems.  These  two  works  list  the  essential  character- 
ics  of  radome  behavior  and  contain  references  to  most  of  the  war¬ 
time  work. 

All  early  radome  work  relies  on  ray  tracing  to  obtain  the  electro¬ 
magnetic  energy  reflected  from  or  transmitted  through  the  radome  wall. 
This  method  was  extended  by  Tricoles  [5]  who  treated  the  receiving 
problem  by  tracing  rays  through  the  radome  to  a  receiving  aperture, 
but  then  obtained  the  antenna  voltage  by  an  empirical  method  based  on 
aperture  response  to  Huygen's  sources.  In  a  later  paper  [6]  he  used 
the  same  methods  to  obtain  the  transmittance  of  an  axially  symmetric 
missile  radome.  Each  of  these  papers  indicate  that  scattering  from 
the  radome  vertex  may  be  responsible  for  some  of  the  measured  radome 
behavior,  and  the  second  paper  includes  a  crude  approximation  of  this 
scattering. 

Analysis  of  radome  effects  on  antenna  performance  has  been  carried 
out  for  both  transmitting  and  receiving  antennas.  The  most  frequently 
calculated  antenna  parameter  in  radome  analysis  is  the  far-field  anten¬ 
na  pattern.  Paris  [7]  has  presented  a  technique  for  determining  the 

transmitting  antenna  pattern  of  a  horn  antenna  covered  by  an  aircraft 
radome.  His  procedure  is  based  on  calculating  the  near-fields  incident 
on  the  radome  by  using  a  previously  developed  aperture  integration  pro¬ 
gram  [8],  treating  the  Incident  fields  as  local  plane  waves,  applying 
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plane-wave  transmission  coefficients,  and  then  treating  the  fields  ex¬ 
terior  to  the  radome  as  equivalent  sources  of  the  radiated  fields.  A 
similar  approacn  is  taken  by  Wu  and  Rudduck  [9]  to  obtain  transmitting 
boresight  errors  for  a  circular  aperture  behind  an  ogive  radome.  They 
simplify  the  determination  of  the  near  fields  incident  on  the  radome 
by  expressing  these  fields  in  a  plane  wave  spectrum.  For  the  uniformly 
illuminated  circular  aperture  they  consider,  this  results  in  an  ana¬ 
lytic  expression  for  the  near  fields. 

Since  virtually  all  radome  work  has  military  application,  few 
results  of  development  and  testing  appear  in  the  open  literature. 

Most  dissemination  of  radome  information  has  been  through  a  series 
of  symposia  begun  in  the  mid-1 950's  under  the  title  of  either  Radome 
Symposium  or  Symposium  on  Electromagnetic  Windows.  These  symposia 
were  originally  sponsored  jointly  by  Ohio  State  University  and  the 
U.S.  Air  Force  and  later  by  the  Georgia  Institute  of  Technology  and  the 
U.S.  Air  Force.  In  recent  years  the  Georgia  Institute  of  Technology 
has  been  the  host  and  sponsor  of  the  symposium.  The  proceedings  of 
these  symposia  provide  an  indication  of  the  work  performed  in  the  ra¬ 
dome  area,  but  due  to  the  unavailability  of  the  proceedings,  and  the 
abbreviated  nature  of  the  papers  in  them,  they  provide  little  concrete 
Information  for  other  researchers  to  draw  upon. 

The  importance  of  radomes  in  both  radar  and  antenna  systems 
is  accented  by  the  inclusion  of  full  chapters  on  radome  design  and 
characteristics  in  a  major  radar  handbook  [10]  and  a  major  antenna 
handbook  [11]. 
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The  most  extensive  treatment  of  radomes  Is  contained  In  the  two 
volume  report  edited  by  Tice  and  Walton  [12].  This  report  was  Intend¬ 
ed  to  survey  and  summarize  all  of  the  radane  work  done  up  to  that 
time.  The  report  covers  all  aspects  of  radome  design  and  testing, 
both  electrical  and  mechanical.  An  extensive  bibliography  is  in¬ 
cluded  with  references  to  both  the  open  and  classified  literature. 

It  Is  unfortunate  that  these  volumes  are  not  widely  available  due  to 
their  restricted  circulation. 

More  recently,  an  entire  handbook  has  been  devoted  to  radome 
research  [13].  This  volume  mainly  treats  materials,  but  there  are 
chapters  on  electrical  and  operational  considerations  and  a  historical 
Introduction.  A  report  on  radome  design  has  also  been  prepared  for 
NATO  [14].  This  report  contains  an  excellent  summary  of  how  environ¬ 
mental,  structural,  electrical,  and  material  considerations  are  com¬ 
bined  In  successful  radome  design.  A  subsequent  NATO  report  [15] 
gives  consideration  to  radome  materials. 

Simply  stated,  the  radome  problem  is  to  find  the  voltage  in  the 
antenna  feed  line  when  a  plane  wave  is  incident  on  the  radome.  An  exact 
solution  of  the  radome  problem  is  extremely  difficult.  The  reasons  for 
this  are  the  complicated  geometry  of  the  problem,  the  boundary  condi¬ 
tions  which  are  imposed  on  the  fields,  and  the  large  size  of  the  struc¬ 
tures  as  measured  in  wavelengths  at  the  operating  frequency. 

The  response  of  an  aperture  antenna  may  be  determined  from 
the  fields  incident  on  the  antenna  and  the  known  characteristics 
of  the  antenna.  The  radome  problem  can  thus  be  considered  as  that 
of  finding  the  incident  fields  at  the  antenna  in  the  presence  of  the 
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radome  and  the  antenna.  The  incident  fields  at  the  antenna  are 
composed  of  two  distinct  parts:  the  fields  passing  directly  through 
the  radome  to  the  antenna  and  the  fields  that  are  scattered  from 
the  various  objects  near  the  antenna  including  the  radome,  metal 
tip,  support  structure  and  various  feed  structures  on  the  antenna. 

The  second  category  implicitly  includes  multiple  scattering  between 
any  combination  of  the  mentioned  objects.  An  exact  solution  to  the 
radome  problem  would  consider  all  of  these  phenomena  simultaneously 
by  obtaining  a  solution  to  the  electromagnetic  field  equations  that 
would  satisfy  the  boundary  conditions  at  all  surfaces  of  material 
discontinuity.  Such  a  problem  could  be  formulated  in  terms  of 
integral  equations  for  the  tangential  field  components  over  a  closed 
surface  surrounding  the  antenna.  The  complicated  form  of  the  boundary 
conditions  and  the  typically  large  electrical  size  of  the  structures 
have  made  exact  techniques  intractable  and  approximate  solutions 
must  be  sought. 

As  previously  mentioned,  the  response  of  a  given  antenna  may 
be  determined  from  the  incident  fields.  The  importance  of  radome 
effects  should  thus  be  ordered  according  to  the  influence  on  the 
Incident  fields  at  the  antenna  and  an  approximate  analysis  will 
consider  only  the  most  important  of  these  effects.  A  primary 
consideration  is  the  effect  of  the  radome  wall  on  the  part  of  the 
incident  field  directly  incident  on  the  aperture.  A  well-designed 
radome  will  be  highly  transparent  and  this  component  will  constitute 
by  far  the  largest  amount  of  energy  incident  on  the  aperture.  Since 
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the  radome  wall  Is  a  curved  surface  with  a  transmission  coefficient 
that  varies  with  angle  of  incidence  and  polarization,  the  effect  of 
the  radome  wa'l  on  the  direct  wave  will  be  attenuation,  depolarization 
and  phase  shift. 

Of  secondary  priority  is  the  incident  field  component  that  is 
singly  reflected  from  the  internal  radome  wall.  This  component  may 
be  quite  significant  at  angles  of  incidence  several  degrees  off  bore- 
sight  where  specular  reflection  from  the  radome  wall  into  the  aperture 
is  favorable. 

Most  airborne  radomes  designed  for  high-speed  flight  will  be 
equipped  with  a  metal  tip  to  increase  the  erosion  resistance  of  the 
radome.  This  tip  is  a  source  of  indirect  scatter  at  all  angles  of 
incidence.  The  support  structure  behind  the  antenna  proper  is  thought 
to  contribute  little  to  the  antenna  voltage,  since  the  feed  is  shielded 
by  the  aperture  and  a  well -designed  feed  will  have  little  spillover  at 
the  aperture.  However,  the  structure  behind  the  antenna  can  contribute 
to  scatter  back  to  the  radome  and  in  turn  to  the  aperture:  multiple 
scatter.  Also  multiple  scatter  can  occur  from  the  feed  and  its  support 
structure. 

The  emphasis  of  this  study  is  on  the  analysis  of  existing 
radomes  as  opposed  to  design  of  radomes  having  desired  specifications. 
Of  primary  Interest  Is  the  prediction  of  the  pointing  error  of  a 
monopulse  radar  operated  in  the  radome.  For  the  type  radomes  of 
interest,  the  transmission  properties  of  the  radome  wall  and 
reflection  from  the  internal  wall  are  considered  to  be  the  major 


contributors  to  pointing  error  with  the  first  predominating  at  small 
gimbal  angles  and  the  latter  at  angles  well  off  the  radome  axis. 

Radar  pointing  error  can  be  a  significant  limitation  to  the 
accuracy  of  missile  guidance  systems.  Also,  the  rate  of  change  of 
pointing  error  with  aspect  angle,  commonly  called  boresight  error 
slope,  will  influence  the  radar's  measurement  of  target  angular 
velocity.  The  analysis  performed  in  this  work  is  directed  toward 
the  prediction  of  these  sources  of  error  in  order  to  determine  their 
impact  on  missile  guidance  and  to  evaluate  methods  for  error  compensa 
tion. 

This  work  is  organized  to  consider  the  fundamental  elements  of 
radome,  antenna,  and  radar  analysis  separately.  This  is  done  in  the 
following  chapter.  The  third  chapter  contains  an  analysis  of  these 
items  as  a  system  with  primary  emphasis  on  prediction  of  angular 
pointing  errors.  Methods  to  compensate  for  angular  errors  are  the 
topic  of  Chapter  IV. 


II.  FUNDAMENTAL  FACTORS  IN 
RADOME  ANALYSIS 


This  chapter  may  be  considered  a  survey  of  the  basic  Information 
needed  for  radome  analysis.  The  approach  taken  here  is  to  examine  the 
subsystems  encountered  in  the  radome-radar  problem  and  some  of  the 
models  used  to  obtain  tractable  solutions  to  realistic  problems.  The 
first  section  in  this  chapter  deals  with  the  theory  of  monopulse  radar 
operating  ideally  without  radome,  signal  distortion,  or  noise.  The 
second  section  introduces  the  radomes  that  will  be  considered  in  this 
work.  A  discussion  of  tangent  ogive  geometry  is  presented.  Next  are 
sections  on  the  transmission  and  reflection  properties  of  planar  slabs 
and  the  technique  of  ray  tracing.  These  basic  optical  methods  will  be 
used  to  model  the  local  properties  of  the  radome  surface.  The  final 
section  of  this  chapter  discusses  ray  tracing  which  forms  the  basis 
for  the  treatment  of  propagation  within  the  radome. 

MONOPULSE  RADAR 

Monopulse  is  a  technique  for  obtaining  target  direction  informa¬ 
tion  from  a  single  radar  return,  hence  the  origin  of  the  term.  The 
theory  of  monopulse  operation  is  well  known  and  discussed  in  many  in¬ 
troductory  radar  books  such  as  Skolnik  [13]  and  Barton  [17].  The  most 
definitive  treatment  of  the  mathematical  theory  of  monopulse  and  the 
underlying  principles  is  Rhodes  early  monograph  [18].  A  recent  volume 
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[19]  edited  by  Barton  contains  reprints  of  many  of  the  key  papers  on 
monopulse  and  is  recommended  reading  for  anyone  wishing  a  historical 
perspective  on  the  developments  in  this  area. 

This  section  is  designed  to  serve  as  an  introduction  to  mono- 
pulse  operation  under  ideal  conditions  and  to  form  a  foundation  for 
the  later  analysis  of  a  monopulse  system  operating  in  an  environment 
containing  a  radome.  In  the  main  body  of  this  work  it  will  be  shown 
that,  in  general,  a  radome  introduces  spatial  distortion  into  the 
wavefront  incident  on  the  receiving  antenna.  This  violates  one  of  the 
key  assumptions  concerning  the  received  wave  and  is  a  source  of  radar 
pointing  error. 

Fundamental  to  the  theory  of  monopulse  is  the  assumption  that 
the  wave  scattered  from  the  radar  target  is  a  uniform  plane  wave  over 
the  receiving  antenna  or  that  the  antenna  is  in  the  far-field 
of  the  target.  When  this  assumption  is  satisfied  the  response  of  a 
receiving  antenna  is  completely  described  by  its  far-field  antenna 
pattern.  The  far-field  antenna  pattern  of  a  radar  antenna  is  easily 
measured  and  analytic  expressions  for  such  patterns  are  not  difficult 
to  obtain. 

The  first  basic  problem  addressed  in  this  study  is  the  determi¬ 
nation  of  a  method  to  obtain  the  response  of  a  monopulse  system  to  a 
non-planar  non-uniform  wave.  This  problem  has  been  addressed  in  the 
literature  in  studies  on  radar  response  to  glint  [20],  multipath  [21], 
and  unresolved  targets  [22].  The  approach  taken  in  each  of  these  studies, 
however,  does  not  answer  the  basic  question  raised  here.  When  analyzing 
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the  above-mentioned  sources  of  a  distorted  wavefront,  the  approach  taken 
has  been  to  decompose  the  target  into  smaller  parts,  each  of  which 
scatters  a  wave  that  is  uniform  and  planar  at  the  receiving  aperture. 
Hence,  the  principle  of  superposition  can  be  applied  to  obtain  the  mono¬ 
pulse  response  in  these  environments.  Treating  a  general  wavefront  as 
a  superposition  of  plane  waves  traveling  in  different  directions  is  a 
well -documented  analytical  technique,  often  known  as  plane-wave  spectrum 
analysis  [23]. 

The  resolution  of  a  received  wave  that  is  distorted  by  a  radome 
into  a  plane-wave  spectrum  representation  is  a  difficult  problem  and  is 
not  justified  if  only  the  monopulse  response  is  desired.  Techniques 
are  presented  in  this  section  for  obtaining  the  monopulse  response 
directly  from  the  fields  incident  on  the  antenna,  be  they  planar,  dis¬ 
torted,  or  whatever.  However,  a  most  elementary  monopulse  system  in 
the  form  of  an  interferometer  is  examined  first.  This  provides  an  easy 
way  to  visualize  how  monopulse  works.  Later,  the  receiving 
properties  of  the  more  practical  parabolic  reflector  antenna  are  pre¬ 
sented. 

Consider  a  receiving  antenna  system  made  up  of  two  elements  sepa¬ 
rated  by  a  distance  D  as  shown  in  Figure  2-1.  Assume  a  time-harmonic 
plane  wave  traveling  in  free  space  is  incident  on  the  receiving  ele¬ 
ments.  Here  and  throughout  this  work  the  time  variation  of  all  time- 
varying  quantities  is  assumed  to  be  eJult  where  j  =  v^T,  u  is  the  fre¬ 
quency  in  rad/s,  and  t  is  the  time  in  seconds.  Furthermore,  the  time 
variation  will  usually  be  suppressed  and  the  conventions  of  ac  circuit 
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Figure  2-1.  Interferometer  ■ionopulse 
System. 
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analysis,  or  phasor  notation,  will  be  employed.  Thus  the  electric 
field  amplitude  of  the  incident  field,  ,  is  given  by 

pi  ,  pi  Jk(x  sine  +  y  cose)  (2.1) 

c  co  e 

where  E*  is  the  magnitude  of  the  field,  k  is  the  propagation  constant 
given  by  2n/x,  x  is  the  wavelength  of  the  wave,  and  e  is  the  incident 
angle  as  indicated.  The  phase  reference  for  the  incident  wave  and  the 
antenna  system  is  taken  to  be  the  origin. 

Induced  potentials  that  are  proportional  to  the  incident  electric 
field  will  exist  across  the  terminals  of  the  receiving  elements.  These 
signals  can  be  written  as 

V1  =  v0  SinQ  (2.2) 


at  element  1  and 


V0  *  V  e‘jl!  sine 
2  o 


(2.3) 


at  element  2  where  VQ  is  the  voltage  magnitude  which  is  identical  for 
element  1  and  2. 

The  distinguishing  feature  of  monopulse  is  the  method  of  com¬ 
bining  V-j  and  V ^  so  that  the  angle  of  incidence  can  be  determined. 
This  is  done  by  adding  and  substracting  the  antenna  voltages  to  form 
two  new  signals  known  as  sum  and  difference  signals.  Thus  the  sum 
signal ,  £,  i >  given  by 


l  ■  V1  *  V2 


(2.4) 
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and  the  difference  signal,  a,  is  given  by 


A  •  V1  -  V2.  (2.5) 

The  method  of  forming  these  signals  is  shown  in  Figure  2-1. 

It  is  convenient  to  introduce  a  new  variable  to  represent  the 
phase  of  and  V2  as  given  in  (2.2)  and  (2.3).  The  symbol,  u,  will 
be  used  to  denote  this  phase  angle  or, 

u  *  k?r  sine  =  sine  ;  (2.6) 

where  it  is  noted  that  u  contains  the  element  spacing  normalized  to 
the  wavelength.  Using  the  variable  u  and  substituting  (2.2)  and  (2.3) 
into  (2.4)  and  (2.5),  the  sum  and  difference  signals  can  be  written 
as 


l  -  v0  ejU  +  Vo  e'JU  ■  (2-7) 

4  *  vo  eJU  -  V0  e-ju  •  (2.8) 

or  by  using  Euler's  identity. 


l  «  2Vq  cos(u)  ,  (2.9) 

4  =  j2VQ  sin(u)  ,  (2.10) 

Signals  aid  V2  together  with  t  and  A  are  conveniently  displayed  on  a 
phasor  diagram  as  shown  in  Figure  2-2.  It  should  be  pointed  out  that 
the  choice  of  the  origin  as  the  phase  center  causes  z  to  be  purely 
real,  and  a  is  in  phase  quadrature  with  z. 
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Incidence  angle  detection  is  accomplished  by  taking  the  real  part 
of  the  ratio  of  -ja  to  l  or. 


-jA 

9a  Re  } 


2V0  sin(u) 
2Vq  cos(u) 


tan(u) 


(2.11) 


The  ratio  Re(-jA/£}  is  known  as  the  monopulse  ratio  and  is  the  output 
signal  of  the  monopulse  system.  Figure  2-3  is  a  plot  of  the  sum  and 
difference  response  and  the  monopulse  ratio.  For  tracking  radars,  the 
angle  of  incidence  will  be  near  the  antenna  axis,  and  thus  u  will  be 
small.  Using  the  small  angle  approximation,  tan(a)3a,  and  the 
definition  of  u,  equation  (2-11)  can  be  written  as 


Re  =  u  =  —  sine 

I  * 


(2.12) 


Again  using  a  small  angle  approximation  (sin(a)  =  a),  (2.12)  can  be 
solved  for  0, 


0  =  lf-1  1  Re  {Hp-}  . 


(2.13) 


Thus,  by  measuring  the  monopulse  ratio  and  multiplying  by  a  constant, 
the  angle  to  the  radar  target  may  be  determined.  This  information  may 
be  fed  to  a  gimbal  system  that  repositions  the  antenna  so  that  a  moving 
target's  position  is  continually  tracked. 

Some  comments  are  in  order  concerni  ig  this  monopulse  system. 


First,  the  system  will  have  zero  output  when  the  amplitude 

and  phase  of  the  voltages  at  elements  1  and  2  are  identical.  If  the 
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phase  of  the  signals  at  elements  1  and  2  are  different,  an  output  signal 
is  produced  that  is  linearly  (for  small  angles)  related  to  the  angle  of 
incidence  of  a  plane  wave  that  would  prodice  the  element  voltages.  In 
this  respect  the  two-element  intereferometer  monopulse  system  is  iden¬ 
tical  to  the  method  used  by  Sims  and  Graf  [24]  for  determining  wave- 
front  distortion  due  to  glint.  Their  method  involved  sampling  the  wave- 
front  at  two  points  separated  by  a  distance  representing  the  receiving 
aperture  diameter  and  using  the  phase  of  the  wave  at  these  points  to 
make  a  linear  approximation  of  the  wavefront  across  the  receiving 
aperture.  The  angular  difference  between  the  normal  to  this  wavefront 
and  the  true  angle  to  the  target  was  then  defined  as  the  radar  pointing 
error  due  to  glint. 

Note  also  that  the  sign  of  the  incidence  angle  is  determined  by 
the  phase  of  the  difference  channel  with  respect  to  the  sum  channel. 

When  the  incidence  angle  is  positive,  A  leads  £  by  90°  in  phase,  and 
lags  £  by  90°  when  the  incidence  angle  is  negative.  The  phase-sensitive 
detector  shown  in  Figure  2-1  indicates  the  sign  of  the  incidence  angle 
by  measuring  this  phase  relationship.  One  final  comment  concerning 
monopulse  systems  is  that  incidence  angle  determination  does  not  involve 
the  amplitude  of  the  incident  wave.  This  is  because  incident  angle 
determination  occurs  in  the  form  of  a  ratio.  Of  course  the  signal 
levels  must  be  significantly  above  the  noise  level  for  reliable  opera¬ 
tion  and  this  is  assumed  throughout. 

Next  the  receiving  properties  of  a  monopulse  system  employing  a  para¬ 
bolic  reflector  antenna  are  examined.  The  parabolic  reflector  is  a  commonly 
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used  antenna  and  has  received  considerable  attention  in  the  literature 
[25,  26].  As  previously  mentioned,  almost  all  antenna  analysis  treats 
the  transmitting  situation  and  obtains  the  fields  at  points  far  removed 
from  the  radiating  structure.  The  receiving  antenna  has  been  treated 
in  general  [27],  but  little  attention  has  been  given  to  the  aperture 
receiving  antenna.  The  purpose  of  the  remainder  of  this  section  is  to 
derive  expressions  for  the  response  of  a  monopulse  receiving  antenna 
excited  by  an  arbitrary  (non-planar,  non-uniform,  randomly  polarized, 
but  time-harmonic)  incident  field.  The  results  are  applicable  to  all 
aperture  antennas. 

The  Lorentz  reciprocity  theorem  will  be  employed  to  find  the 
receiving  characteristics  of  the  parabolic  reflector  antenna.  In  order 
to  use  this  theorem,  the  behavior  of  the  antenna  in  the  transmitting 
mode  must  be  known.  This  is  most  easily  done  through  the  use  of  the 
equivalence  principle,  that  is,  by  finding  a  set  of  sources  which  produce 
the  same  fields  as  are  produced  by  the  antenna  in  the  region  of  interest. 
As  will  be  seen,  the  equivalent  sources  play  an  important  role  in  the 
monopulse  receiving  system. 

The  parabolic  reflector  antenna  system  is  composed  of  two  princi¬ 
pal  parts,  a  feed  antenna  and  a  reflecting  surface  in  the  shape  of  a 
paraboloid  of  revolution.  As  shown  in  Figure  2-4,  the  feed  antenna, 
the  primary  radiator,  is  located  at  the  focus  of  the  parabola.  The 
spherical  wavefront  leaving  the  feed  antenna  will  be  reflected  at  the 
parabola  and  become  a  planar  wavefront  traveling  parallel  to  the 
antenna  axis. 
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Figure  2-4.  Parabolic  Reflectrr  Antenna  System 
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The  planar  circular  area  that  is  perpendicular  to  the  antenna  axis 
and  has  its  perimeter  at  the  reflector's  edge  is  called  the  aperture 
plane.  The  fields  over  the  aperture  plane  can  be  used  to  determine  the 
transmitting  and  receiving  characteristics  of  the  antenna.  This  is 
shown  first  for  the  transmitting  case  where  the  fields  produced  by  the 
antenna  in  region  2  >  0  are  to  be  found. 

The  region  z  >  0  is  source  free,  so  by  the  uniqueness  theorem  [28] 
the  fields  in  this  region  are  uniquely  related  to  the  tangential  field 
components  over  a  closed  surface  that  contains  the  actual  sources.  It 
is  convenient  to  enclose  the  sources  by  the  surface  formed  by  the  xy- 
plane  and  a  large  hemisphere  with  center  at  the  origin.  The  hemisphere 
may  be  thought  of  as  being  of  infinite  radius  or  of  finite  radius  but 
large  enough  so  that  the  fields  over  the  hemisphere  are  negligible. 

Thus  the  fields  in  the  region  2  >0  are  uniquely  related  to  the  tangen¬ 
tial  field  components  over  the  xy-plane. 

The  electric  and  magnetic  field  intensity  produced  by  the  antenna 
will  be  denoted  by  E3  and  H® ,  and  it  will  be  assumed  that  F3  and  H®  are 
zero  in  the  region  z  <  0,  that  is,  behind  the  reflector.  This  will 
occur  when  the  feed  does  not  produce  fields  in  this  region  (no  spill¬ 
over  in  the  terms  of  antenna  engineers),  and  when  scattering  from  the 
reflector  edge  is  neglected.  Thus,  all  the  antenna  fields  in  the  xy- 
plane  are  confined  to  the  aperture  plane. 

The  equivalence  principle  [29]  will  next  be  used  to  further  simplify 
the  antenna  analysis.  As  guaranteed  by  the  uniqueness  theorem,  the 
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fields  In  the  region  z  >  0  are  now  completely  determined  by  the  tangen¬ 
tial  field  components  In  the  aperture  plane.  According  to  the  equiva¬ 
lence  principle,  these  same  aperture  fields  may  be  supported  by  equiva¬ 
lent  surface  currents,  3d  and  M®,  in  the  aperture  plane,  and  all  the 
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fields  in  the  region  z  <  0  may  be  set  to  zero.  The  equivalent  currents 
are  found  from, 

Tsa  *  n  x  H®  ,  (2.14) 

and 

ffsa  s  E3  x  n  ,  (2.15) 

a 

where  E®  and  H®  are  the  fields  produced  by  the  antenna  and  n  is  a  unit 
vector  normal  to  the  aperture  plane  as  shown  in  Figure  2-4. 

The  equivalent  surface  currents  given  by  (2.14)  and  (2.15)  follow 
from  the  electromagnetic  boundary  condition  equations  which  state  that 
discontinuities  in  tangential  field  components  are  associated  with  surface 
current  densities.  At  a  smooth  interface  separating  regions  labeled 
1  and  2  these  equations  are: 

Js  »  n  x  [if,  •  H2]  ,  (2.16) 

=  [Tj  -  Eg]  x  n  ,  (2.17) 

where  subscripts  1  and  2  indicate  the  fields  at  the  interface  in  region 
1  and  2  respectively,  F  represents  the  electric  field  intensity,  H  the 
magnetic  field  intensity,  T$  is  the  electric  surface-current  density,  M$ 
is  the  magnetic  surface-current  density,  and  n  is  a  unit  vector  normal 


to  the  Interface  directed  into  region  1.  Thus  (2.16)  states  that  elec¬ 
tric  surface-current  densities  are  associated  with  discontinuities  in 
tangential  magnetic  fields  and  (2.17)  states  that  magnetic  surface- 
current  densities  are  associated  with  discontinuities  in  the  tangen¬ 
tial  electric  fields. 

By  choosing  region  2  to  be  the  region  z  <  0  and  by  setting  all 
the  fields  in  this  region  to  zero,  (2.14)  and  (2.15)  are  seen  to  follow 
from  (2.16)  and  (2.17).  The  fields  in  the  region  z  >  0  are  the  same 
as  before  but  now  they  are  supported  by  the  "equivalent  currents",  Ts 
and  ffs.  Thus  the  entire  problem  of  finding  the  antenna  fields  is  reduced 
to  the  problem  of  finding  the  fields  of  the  currents  given  by  (2.14)  and 
(2.15). 

It  should  be  noted  that  the  above  treatment  of  the  reflector  antenna 
results  in  the  same  formulation  as  the  classical  optics  problem  of 
diffraction  from  a  circular  aperture  in  a  conducting  screen.  All  such 
aperture  problems  may  be  solved  by  the  same  method  based  on  the  equiva¬ 
lence  principle. 

A  comnent  concerning  the  antenna  feed  is  in  order  before  proceeding. 
As  seen  in  Figure  2-4,  part  of  the  field  radiated  by  the  feed  will  be 
reflected  back  to  the  feed  and  interact  with  it.  This  is  often  referred 
to  as  aperture  blockage  and  may  have  a  significant  effect  on  the  perfor¬ 
mance  of  the  antenna  system,  particularly  when  the  size  of  the  feed  is 
an  appreciable  part  of  the  reflector  size.  Consideration  of  aperture 
blockage  is  beyond  the  scope  of  this  work  and  it  will  be  assumed  that 
the  antenna  system  can  be  represented  completely  by  the  equivalent 
currents  as  previously  discussed. 
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There  are  several  ways  of  determining  the  fields  of  the  antenna 
from  the  equivalent  currents.  These  include  direct  integration,  use  of 
potential  functions,  and  the  plane-wave-spe:trum  technique.  The  fields 
radiated  by  the  antenna  will  not  be  found  here  since  it  Is  the  antenna's 
receiving  properties  that  are  of  greatest  interest.  These  receiving 
characteristics  can  be  derived  from  the  Lorentz  reciprocity  theorem. 

The  Lorentz  reciprocity  theorem,  which  follows  directly  from 
Maxwell's  equations,  is  a  way  of  relating  two  sets  of  sources  and  their 
fields.  An  integral  form  of  this  theorem  is  given  by  [30], 

/os-  •  ^  -  H®  .  M*3)  dv  =  jh**  •  J®  -  H*3  •  H®)  dv  (2.18) 

vol  vol 

where  superscripts  indicate  a  set  of  sources  and  associated  fields  as 
used  previously,  and  the  integrations  extend  over  all  space.  Rumsey  [31] 
has  given  the  name  reaction  to  the  relation  between  sources  and  fields 
expressed  by  either  side  of  (2.18).  The  recripocity  theorem  thus  states 
that  the  reaction  of  the  field  a  on  source  b  is  equal  to  the  reaction  of 
the  field  b  on  the  source  a. 

In  order  to  apply  (2.18)  to  the  receiving  antenna  problem  the  a 
source  and  field  are  chosen  to  be  that  of  the  equivalent  currents 
representing  the  parabolic  reflector  antenna  and  the  b  source  is  chosen 
to  be  a  point  electric  current  element  (dipole)  given  by, 

O*3  *  JQb  6(|R  -  R'|)  , 


(2.19) 


where  3^  contains  the  magnitude  and  direction  of  the  current  element, 

R‘  Is  a  vector  from  the  origin  to  the  element's  location,  R  is  the  posi¬ 
tion  vector  (R  =  xx  +  yy  +  zz),  and  6  is  the  Dirac  delta  function.  M*5 
Is  set  to  zero. 

Substituting  (2.19)  into  the  left-hand  side  of  (2.18)  and  recog¬ 
nizing  the  sampling  property  of  the  Dirac  delta  function  yields. 


If  3"*1  is  chosen  to  have  unit  amplitude  and  to  be  in  the  direction  of  E0 
then  (2.20)  is  simply  the  magnitude  of  the  electric  field  at  the  point 
dipole  produced  by  the  sources  3d  and  M*.  From  the  right-hand  side  of 
(2.18),  it  is  seen  that  the  magnitude  of  E3  at  E' ,  as  given  by  (2.20), 
is  equivalent  to  the  reaction  of  the  fields  of  the  unit-amplitude  dipole 
on  the  equivalent  currents  of  the  parabolic  antenna. 

The  behavior  of  a  transmitting  antenna  at  distances  much  greater 
than  the  antenna  dimensions  (In  the  far-field)  is  completely  charac¬ 
terized  by  the  electric  field  it  radiates,  due  to  the  transverse  nature 
of  the  radiated  fields.  The  far-field  response  of  a  transmitting  anten¬ 
na,  defined  as  the  magnitude  of  the  electric  field,  is  obtained  from 
(2.20)  when  |R'|  is  large.  The  response  of  the  same  antenna  in  the 
receiving  mode,  defined  as  the  open-circuit  voltage  at  the  antenna 
terminals,  is  proportional  to  the  right-ham'  side  of  (2.18)  [31].  In 
the  remainder  of  this  section  the  receiving-mode  response  of  the  para¬ 
bolic  antenna  is  discussed. 
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The  right-hand  side  of  (2.18)  can  be  rewritten  using  the  aperture 
fields  as  given  by  (2.14)  and  (2.15)  yielding, 


yV6  .  (n  x  Ha)  -  H*5  •  (E®  x  n)  ds, 
A 


(2.21) 


where  A  indicates  that  tne  integration  is  only  over  the  antenna  aperture 
since  the  a  currents  are  zero  everywhere  else.  Throughout  this  study 
it  is  assumed  that  E®  and  H®  are  of  constant  phase  and  related  like  a 
plane  wave  in  the  aperture.  For  a  linearly  polarized  antenna  E®  will 
have  the  same  direction  over  the  aperture  and  H®  will  also  be  uni¬ 
directional  and  at  right  angles  to  E®. 

As  an  example,  assume  that  E®  in  the  aperture  is  given  by. 


1®  =  E  ®  X  (2.22) 

o 

and  H®  is 


where  zero  subscripts  indicate  magnitude.  The  second  part  of  (2.23) 
is  obtained  from  the  plane-wave  nature  of  E®  and  H®,  and  n  is  the  intrin¬ 
sic  impedance  of  free  space  (n  3  /y/e).  Figure  2-5(a)  illustrates 
these  aperture  fields. 
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Substituting  these  a  fields  into  (2.21)  yields, 

f •  (hQ*x)  +  H*  •  (E0a  y)  ds  ,  (2.24) 


where  the  cross  products  have  been  evaluated  for  n  =  z. 

EQa  is  common  to  both  terms  in  (2.24)  and  can  be  factored  to  yield 


[H6-  i  - 1  ?  •  i]  ds. 


(2.25) 


Thus  it  is  seen  that  Ea  acts  as  a  weighting  function  for  the  integration. 
Also,  since  Ea  has  constant  phase,  it  can  be  assumed  purely  real.  By 
adjusting  the  variation  of  Ea  over  the  aperture,  various  antenna  pat¬ 
terns  may  be  obtained.  In  practice  this  is  accomplished  by  design  of 
the  antenna  feed. 

The  importance  of  the  form  of  Ea  will  be  emphasized  by  giving  it 
the  name  aperture  illumination  function  and  using  the  symbol,  g,  for 
this  function.  The  monopulse  sum  and  difference  patterns  are  produced 
by  the  proper  aperture  illumination  functions.  To  see  how  this  is 
accomplished  it  is  first  necessary  to  find  the  antenna  response  to  an 
incident  plane  wave. 


To  obtain  the  antenna  response,  assume  that  a  distant  point  dipole 
produces  a  field  which  is  incident  on  the  aperture  as  in  Figure  2- 5(b). 
This  field,  which  is  denoted  by  b,  will  be  a  uniform  plane  wave  over 
the  aperture  and  can  be  written  as, 


jksinex 


e 


(-cosex  +  sinez) 


(2.26) 


2 


(a)  Aperture  fields 


Figure  2-5.  Fields  of  Parabolic  Antenna  and  Sampling  Source. 
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and 


H*5 


E  b  jksinex  « 
-9-e  (y) 

n 


(2.27) 


where  E  b  is  the  magnitude  of  the  electric  field, 
o 

Substituting  (2.26)  and  (2.27)  into  (2.25),  performing  the  dot  pro¬ 
ducts,  and  using  the  symbol  for  the  aperture  illumination  function 
yields, 

/E  b  i  .  jksinex 

g[^- +  1  E  b  cos  e]  e  ds  (2.28) 

o  no 
A 


By  removing  the  constants  from  the  integral  and  making  the  substitution, 
kx  =  k  sine,  (2.28)  becomes. 


-  E  1 
n  o 


(1  +  cose 


>/ 


jkyx 

g  e  ds. 


(2.29) 


The  integral  in  (2.29)  is  immediately  recognized  as  a  Fourier  transform 
of  F  and  expresses  a  well-known  relation  between  the  aperture  illumination 
function  and  the  far-field  antenna  pattern. 

The  integration  in  (2.29)  is  most  easily  performed  in  circular 
cylindrical  coordinates  for  the  circular  aperture  under  discussion.  Writing 
the  integral  completely,  (2.29)  becomes. 


2-tr  0/2 

^  E^b  (1  +  cose)  f  J  g(p .♦) 
o  o 


jk  pcos<(> 
e  pdpd<j> , 


(2.30) 


i 
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where  the  functional  dependance  of  g  is  now  shown  and  p  and  <t>  are  the 
coordinates  indicated  in  Figure  2-5.  In  order  to  normalize  this  expres¬ 
sion  as  was  done  for  the  interferometer,  the  following  variable  change 
is  introduced, 

r  *  2p/D  ,  (2.31) 

u  *  &  sine.  (2.32) 

Using  these  variables,  expression  (2.30)  is  written 

i  k  n2  r  r  jurcos* 

“  E0b(l  +  cose)  J  J  J  g(r,$)  e  rdrd<|>  (2.33) 

o  o 

where  g(r,<(>)  replaces  g(p,$)  according  to  (2.31). 

As  with  the  interferometer,  an  aperture  antenna  used  with  a  mono¬ 
pulse  system  will  have  two  far-field  response  patterns  which  are  again 
called  the  sum  and  difference  patterns.  These  patterns  are  deter¬ 
mined  by  the  aperture  illumination  function.  The  sum  pattern  is  an  even 
function  of  e  and  the  difference  pattern  is  an  odd  function  of  e.  It 
has  been  shown  that  the  far-field  pattern  of  the  aperture  involves  the 
Fourier  transform  of  the  aperture  illumination  function  (2.29). 

According  to  elementary  Fourier  transform  theory  [32],  the  transform  of 
an  even  function  is  even  and  the  transform  of  an  odd  function  is  odd. 
Thus  the  aperture  illumination  functions  for  the  sum  pattern  will  be 
even,  and  the  aperture  illumination  function  for  the  difference  pattern 
will  be  an  odd  function  of  x. 
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In  this  study  two  comnonly  used  sum-pattern  aperture  illumination 
functions  will  be  considered.  The  first  is  uniform  illumination,  or 

g(p>+)  *  1  ,  p  <  \  .  (2.34) 

This  illumination  is  sometimes  called  ideal  because  it  is  constant  to 
the  edge  of  the  aperture  but  zero  outside.  It  would  be  impossible  to 
create  in  practice  due  to  the  discontinuity  in  the  fields  at  the 
aperture  edge. 

The  second  sum  illumination  function  is  tapered  from  the  center  of 
the  aperture  to  the  edge  according  to  a  cosine  function, 

g(p,«)  •  cos(Ifr)  ,  p  <  |  .  (2.35) 

The  integral  part  of  (2.33)  is  called  the  aperture  pattern.  It  will 
be  shown  that  for  large  apertures  the  aperture  pattern  contains  all  the 
significant  variation  of  (2.33)  with  0.  The  aperture  patterns  for  the 
two  aperture  illumination  functions  are  shown  in  Figure  2-6.  These 
patterns  are  normalized  to  their  respective  maxima  and  are  plotted  in 
dB  versus  u.  It  is  noted  that  the  tapered  illumination  function  pro¬ 
duces  a  wider  main  lobe  but  has  lower  side-lobe  levels  than  the  uniform 
illumination.  These  aperture  patterns  could  be  obtained  analytically, 
by  using  a  two-dimensional  fast  Fourier  transform  algorithm,  or  by 
numerical  integration.  Figure  2-6  and  the  remaining  figures  in  this 
section  were  obtained  using  the  numerical  integration  technique 
that  is  used  to  find  the  response  of  the  antenna  with  radome. 

Figure  2-6  may  be  regarded  as  an  "universal"  aperture  pattern  since 
the  response  for  any  size  aperture,  any  wavelength,  and  any  incidence 
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Figure  2-6.  Normalized  Antenna  Pattern  of 
ar  Aperture  with  Uniform  and  Cosine  Illumination 
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angle  can  be  looked  up.  This  study  is  concerned  with  electrically 
large  apertures,  and  as  an  example  an  aperture  with  a  diameter  of  8  wave¬ 
lengths  will  be  used.  The  statement  that  :he  aperture  pattern  contains 
the  significant  angular  variation  of  the  antenna  response  is  now  justi¬ 
fied.  The  response  function  given  by  (2.33)  is  the  product  of  two 
terms  that  vary  with  9:  the  term  (1  +  cose)  and  the  aperture  pattern. 
Figure  2-7  is  a  plot  of  these  terms  in  dB  versus  incidence  angle,  8. 

Each  curve  is  normalized  to  its  maxima.  It  is  seen  that  the  term 
(1  +  cose)  is  only  3dB  down  at  e  =  65°  and  a  maximum  of  6dB  down  at 
90°,  while  the  aperture  pattern  is  greater  than  30dB  down  for  all 
angles  greater  than  about  30°.  Since  the  antenna  response  in  dB  is 
the  sum  of  these  two  terms,  little  error  is  introduced  by  ignoring 
the  (1  +  cose)  term. 

A  method  presented  by  Barton  and  Ward  [33]  will  be  used  to  obtain 
the  difference  pattern  illumination  functions.  They  note  that 
the  difference  pattern  of  monopulse  antennas  resembles  the  derivative 
of  the  sum  pattern.  If  this  relation  is  taken  to  be  exact,  the  difference 
illumination  function,  denoted  by  g^,  is  related  to  the  sum  illumination 
function,  g,  by. 


gd(x,y)  =  -p  g(x,y)  ,  (2.36) 

where  it  is  noted  that  gd  is  the  product  of  the  sum  illumination  and 
the  linear  odd  function,  x;  Lg  is  a  constant.  Since  g  is  even,  gd 


is  odd. 
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Figure  2-7.  Normalized  Antenna  Pattern  of  Circular 
Aperture  with  Uniform  Illumination  and  the  Term  (1  +  cos  e) 
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The  constant,  Ls>  is  the  effective  aperture  illumination  power 


width  in  the  terminology  of  [33]  and  is  given  by 

[f  (2irX)2  |g(x.y)l2  ds  *1  H 


f  |g(x,y) r 


(2.37) 


For  the  case  of  uniform  illumination  in  the  circular  aperture,  L$  is 
Dtt/2,  and  for  cosine  tapered  illumination  (2.37)  yields  Lg  =  1.071867D. 

When  the  difference  illumination  is  defined  by  (2.36)  with  L$ 
given  by  (2.37),  the  sum  and  difference  illumination  functions  contain 
equal  power.  This  would  be  expected  in  a  wel 1 -designed  monopulse 
antenna. 

Using  (2.36)  and  the  above  constants,  the  difference  illumination 
for  uniform  sum  illumination  is. 


4x  D 

9d  =  T  *  P  <  2  ; 

and  for  cosine  illumination  is, 


(2.38) 


5.861962X 


cos(^q-)  ,  p  <  2 


(2.39) 


where  a  mixed  coordinate  system  is  used  for  convenience. 


The  sum  and  difference  illumination  functions  and  the  far-field 
response  function  of  an  8x  diameter  aperture  are  shown  in  Figure  2-8 
for  uniform  illumination,  and  in  Figure  2-  3  for  cosine  tapered  illumina 
tion.  The  response  functions  are  again  normalized  to  the  maximum  of  the 


sum  patterns,  but  are  now  shown  with  a  linear  ordinate  scale. 


!-8.  Sum  and  Difference  Patterns  of  Circular 
Aperture  with  Uniform  Illumination. 
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Monopulse  angle  detection  is  accomplished  by  connecting  the  sum  and 
difference  outputs  of  the  parabolic  antenna  to  an  angle  detection  circuit 
as  was  shown  in  Figure  2-1.  The  angle  of  an  incident  plane  wave, 
is  determined  from  the  monopulse  ratio,  or 

9.  =  k"1  Re{^->  (2.40) 

where  K  i.s  the  slope  of  the  normalized  difference  pattern  evaluated  at 
e  *  0.  For  the  case  of  uniform  sum  illumination,  the  value  of  K  is  given 
by  Barton  and  Ward  to  be  1.573  0/a  for  measured  in  radians.  For  the 
cosine-tapered  illumination,  K  can  be  shown  to  be  1.475  D/A. 

The  incidence  angle  equations  for  the  8a  diameter  aperture  considered 
in  the  example  become, 

0i  *  (12.585)"1  Re{^>  ,  (2.41) 

for  the  uniform  illumination,  and 

ei  -  (11.800)"1  Re{^}  (2.42) 

for  the  cosine  tapered  illumination,  where  0.-  is  in  radians.  These  two 
equations  and  the  methods  outlined  above  will  be  used  to  find  the  indi¬ 
cated  target  angle  in  the  next  chapter. 

To  conclude  this  summary  of  monopulse  principles,  it  is  pointed 
out  that  the  interferometer  may  be  considered  as  an  aperture  illumination 
function  given  by, 

9  »  Eq  «[|R  -  (±§,  0,  0)|]  , 


(2.43) 


since  the  interferometer  samples  the  incident  wave  at  the  points  (±|-,  0,  0) 
as  described  earlier.  It  is  useful  to  keep  this  simple  idea  of  wavefront 
sampling  in  mind  while  considering  the  behavior  of  radomes. 

TANGENT  OGIVE  RADOME  GEOMETRY 

The  geometry  and  coordinate  system  of  the  tangent  ogive  radome  are 
discus -cd  in  this  section.  The  tangent  ogive  is  used  for  the  radome 
shape  throughout  this  study  because  it  well  represents  the  nose  section 
of  high-speed  missiles.  For  minimum  drag  the  high-speed  missile  must 
be  highly  streamlined.  With  proper  parameter  selection,  the  ogive  makes 
a  satisfactory  shape  for  the  missile  nose. 

The  tangent  ogive  is  a  surface  of  revolution  generated  by  the  arc 
of  a  circle  as  shown  in  Figure  2-10.  The  axis  of  the  ogive  is  chosen 
to  be  the  y-axis,  and  the  generating  arc  is  centered  on  the  z-axis  in 
order  to  obtain  tangency  along  the  contour  where  the  ogive  joins  the 
missile  body.  The  radius  of  the  ogive,  p,  may  be  expressed  as  a  func¬ 
tion  of  position  along  the  axis,  y,  and  the  generating  parameters  by 


where  R  is  the  radius  of  the  generating  arc,  and  W  is  the  base  diameter. 

A  measure  of  an  ogive's  streamlining  is  the  length-to-diameter  ratio, 
£/W.  This  ratio  is  encountered  in  aerodynamic  studies  where  it  was  given 
the  name  "fineness  ratio'1.  The  greater  the  fineness  ratio  the  more 
streamlined  is  the  ogive. 
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It  is  convenient  to  express  the  length  of  the  ogive  generating  arc 
in  terms  of  the  fineness  ratio  and  the  ogive  base  diameter,  W.  This  can 
be  done  by  considering  the  right  triangle  0T0‘  in  Figure  2-10.  The 
sides  of  0T0'  are  related  by, 

R2  =  l2  +  [R  -  j]2  ,  (2.45) 

which  can  be  solved  for  R  yielding, 

R  *  (£)2W  +  j  .  (2.46) 

This  expression  is  useful  for  calculating  R  for  an  ogive  with  specific 
fineness  ratio  and  diameter. 

The  shape  of  the  tangent  ogive  is  sometimes  specified  by  the  para¬ 
meter,  caliber.  The  caliber,  C,  is  defined  as  the  ratio  of  the  genera¬ 
ting  arc  length  to  the  base  diameter  or, 

C  =  R/W  (2.47) 

In  addition  to  knowing  the  ogive  shape,  the  unit  vector  normal  to 
the  ogive  surface  is  required  for  the  incidence  angle  calculations  dis¬ 
cussed  later.  It  is  most  convenient  to  express  this  unit  normal,  n, 
as  a  sum  of  three  rectangular  components,  or 

^  A  A  A 

n  =  nxx  +  n^y  +  nzz 


(2.48) 
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where  n  ,  n  ,  n  are  the  components  of  n  in  the  coordinate  directions 
x  y  z 

A  A  A 

x,  y,  and  z  respectively.  These  components  can  be  found  from  geometri¬ 
cal  considerations,  or  by  evaluating  vf/jvf|  with  f  given  by  the  surface 
equation.  The  component  of  n  in  the  axial  direction,  ny,  is  given  by 

ny  *  y/R  .  (2.49) 

The  component  of  n  perpendicular  to  y,  n^,  is  then  found  from  the  re¬ 
quirement  of  unit  magnitude,  or 

» 

nr  -  [1  -  ny2]\  (2.50) 

and  nx  and  nz  are  given  by, 

nx  =  nr  x/p  *  (2.51) 

nz  *  nr  z/p  .  (2.52) 

The  radomes  considered  in  this  study  have  uniform  thickness  and  fill 
the  region  between  two  concentric  ogives.  Thus  a  radome  with  thickness, 
t,  and  an  inner  surface  formed  by  an  ogive  with  generating  radius  R,  has 
as  an  outer  surface  an  ogive  with  generating  radius  (R  +  t). 

TRANSMISSION  AND  REFLECTION  PROPERTIES  OF  PLANAR  DIELECTRIC  SLABS 

In  this  section  the  transmission  and  reflection  coefficients  for 
plane  waves  incident  on  an  infinite  planar  dielectric  slab,  or  several 
slabs  sandwiched,  are  derived.  Virtually  all  analyses  of  radomes  have 
used  such  plane-wave  transmission  coefficients  to  relate  the  electro¬ 
magnetic  field  on  opposite  sides  of  a  radome  wall;  and  good  results  are 
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obtained.  The  technique  of  modeling  the  radome  wall  is  discussed  in 
the  following  section.  Although  transmission  coefficients  for  planar 
slabs  have  been  presented  in  the  literature  in  varying  degrees  of  com¬ 
pleteness  and  generality  (Collin  [34]  being  one  of  the  better  examples); 
they  are  derived  here  in  order  to  define  the  notation  and  demonstrate 
the  method  of  solution.  The  technique  used  here  is  similar  to  the 
equivalent  circuit  approach  used  by  Robinson  [35]. 

The  problem  of  transmission  and  reflection  of  a  plane  wave  inci¬ 
dent  on  a  planar  dielectric  slab  is  an  electromagnetic  boundary  condi¬ 
tion  problem.  Transmission  and  reflection  of  an  incident  wave  are 
necessary  results  of  the  boundary  conditions  at  surfaces  of  material 
discontinuity.  Complex  transmission  and  reflection  coefficients  can 
be  used  to  express  the  fractions  of  the  incident  fields  that  are 
respectively  transmitted  through  and  reflected  from  the  dielectric 
slab  as  well  as  giving  the  phase  of  these  fields.  The  general  equa¬ 
tion  for  a  plane  wave  is  used  to  express  all  the  fields. 

Electromagnetic  fields  in  homogeneous,  isotropic,  time-invariant, 
source-free  regions  satisfy  the  vector  Helmholtz  equation, 

V2F  -  ye  ^=0  (2.53) 

at 

where  F  is  either  the  electric  field  intensity  or  the  magnetic  field 
intensity,  v2  is  the  vector  Laplacian  operator,  time  dependence  is 
e^wt,  y  is  the  permeability  of  the  medium,  and  l  is  the  complex  permit¬ 


tivity  of  the  medium  given  by 


43 


e  *  e(l  -  j  tan  s). 


(2.54) 


with  e  the  permittivity  and  tan  6  the  loss  tangent  of  the  medium 
(tan  &  *  a/u>e).  The  use  of  the  complex  permittivity  in  (2.53)  accounts 
for  attenuation  of  the  fields  due  to  loss. 

The  solution  to  (2.53)  for  fields  that  vary  in  only  one  direction 
(plane  waves)  may  be  written  as. 


F=Fo  e 


jkar 


(2.55) 


where  an  is  a  unit  vector  in  the  direction  of  wave  motion;  r  =  xx  +  yy  + 
k  *  w^ue;  and  Fq  is  a  complex-vector  constant  representing  the  ampli¬ 
tude,  direction,  and  phase  of  F  or  FT  at  the  origin;  and  the  time  depen¬ 
dence  is  suppressed. 

It  will  be  useful  to  define  the  vector  k  as 


k  *  k  »n  '  kx*  +  V  *  kz*-  (?-56) 

This  vector  is  in  the  direction  of  the  Poynting  vector  and  is  sometimes 
called  the  wave  vector.  The  components  of  k  must  satisfy 

kx2  +  ky2  +  kz2  *  u2ye  '  (2.57) 

This  relation  provides  a  very  convenient  way  to  treat  the  boundary 
condition  problem  of  a  plane  wave  incident  on  a  planar  dielectric  slab. 
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The  geonetry  used  for  a  dielectric  si  ib  which  may  be  located 
between  two  dissimilar  regions  is  shown  in  Figure  2-11.  For  the 


Figure  2-1 1 .  Dielectric  Slab  and  Definition  of  Fields 
(Parallel  Polarization). 


45 


practical  case  of  a  single  layer,  the  regions  z  <  0  and  z  >  d  will 
contain  the  same  material,  usually  air. 

The  fields  in  the  various  regions  are  classified  according  to 
location  and  direction  of  travel  as  presented  in  Table  2-1.  Super 
scripts  identify  each  category. 


Table  2-1 

Classification  of  Dielectric  Slab  Fields 


Name 

Location 

Incident  Fields 

z  <  0 

i1,  i1 

Positive  Internal  Fields 

0  <  z  <  d 

Eu,  Hu 

Negative  Internal  Fields 

0  <  z  <  d 

> 

i  rc 

> 

1  LU 

Transmitted  Fields 

z  >  d 

Et,  H* 

Reflected  Fields 

z  <  0 

Er,  Hr 

These  classifications  are  straightforward  but  a  comment  is  in 
order  concerning  the  internal  fields.  When  a  wave  is  incident  at  z  =  0, 
in  general,  there  is  a  field  transmitted  through  the  interface  at 
angle  e2.  This  field  will  travel  to  the  z  =  d  interface  and  be 
partially  reflected,  return  to  z  =  0  and  again  be  partially  reflected 
and  so  on.  The  notation  Eu  and  Hu  means  the  superposition  of  all  fields 
in  the  slab  traveling  in  the  generally  positive  z  direction.  All  of 
these  fields  will  have  the  same  wave  vector,  ku,  as  will  seen  by  Snell's 
reflection  law.  Similarly  Ev  and  Hv  mean  the  sums  of  all  fields  travel¬ 
ing  in  the  generally  negative  z  direction.  These  definitions  are  neces¬ 
sary  to  solve  the  "steady-state"  problem. 
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Two  orthogonal  polarization  states  can  be  considered  for  the  inci¬ 
dent  field.  The  polarization  directions  are  with  respect  to  the  plane 
of  incidence:  the  plane  containing  the  normal  to  the  slab  and  k1.  Per¬ 
pendicular  polarization  means  E1  is  perpendicular  to  the  plane  of  inci¬ 
dence,  and  parallel  polarization  means  E1  is  contained  in  the  plane  of 
incidence.  An  arbitrarily  polarized  incident  field  can  be  expressed  as 
a  superposition  of  a  parallel  polarized  E1  and  a  perpendicular  polarized 
E1. 


Parallel  polarization  will  be  considered  first,  then  the  transmis¬ 
sion  and  reflection  coefficients  for  perpendicular  polarization  will 
be  obtained  using  the  principle  of  duality.  The  equations  for  the 


parallel  polarized  plane  waves  listed 


in  Table  2-1  are  as  follows: 


•cosei  y  + 

sinei  z) 

(2.58) 

:) 

(2.59) 

coser  y  - 

siner  z) 

(2.60) 

x) 

(2.61) 

(-coseu  y 

+  sineu  z) 

(2.62) 

(x) 

(2.63) 

(-:ose  y  -  sine  z)  (2.64) 
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k 


-j(k£y  +  k^(z-d)) 
-j(k£y  +  kj(z-d)) 


(-coset  y  -  sinet  z) 

00 


(2.66) 

(2.67) 


where  the  subscript,  "o“,  indicates  the  quantity  is  a  complex  constant 
and  Zn  is  the  intrinsic  impedance  of  the  ntfl  region  as  given  by 


(2.68) 


Z  is  the  ratio  of  transverse  components  of  E  and  H  in  a  uniform  plane 
n 

wave  traveling  in  the  medium  of  the  ntn  region.  Note  that  the  phase 
reference  of  the  transmitted  fields  is  taken  at  (0,0, d). 

Relations  between  the  fields  listed  above  are  obtained  from 
boundary  conditions  (2.16)  and  (2.17).  Since  no  surface  currents  will 
be  present  on  dielectric  sheets,  these  boundary  conditions  reduce  to 
continuity  of  tangential  fields  or. 


tl  =  Et2 

Htl  =  Ht2 

at  z  =  0  , 

(2.69) 

t2  *  Et3 

Ht2  =  Ht3 

at  z  =  d  , 

(2.70) 

where  the  subscript  t  indicates  total  tangential  field,  and  the  number 
refers  to  the  region. 

For  parallel  polarization  (2.69)  and  (2.70)  refer  to  the  y-  compo¬ 
nents  of  E  and  the  x-components  of  H.  In  terms  of  the  field  components 
(2.70)  can  be  written  as 
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+  V  ■  Eyl 

Iz  »  d 

(2.71) 

*  H*V  *  Hx‘ 

lz  *  d 

(2.72) 

Substituting  (2.62)  -  (2.67)  into  (2.71)  and  (2.72)  yields 

u  -j(k  Uy  +  k  Ud)  v  -j(k  Vy  -  k  vd) 

E„u  e  y  2  (-cose..)  +  Env  e  y  (-coseJ 


u#  o 


-  j  k  *y 

EQt  e  y  (-coset)  (2.73) 


q'  E  u  +  k“d)  -  Zj1  E  v  e'J(^  ‘  ^  •  Zj1  Ej  e-Jk> 

(2.74) 

Equations  (2.73)  and  (2.74)  must  hold  for  any  y,  given  the  follow¬ 
ing  relation. 


y  y  ky  * 


(2.75) 


which  can  be  written  in  terms  of  the  wave  angles  as 


ku  sineu  =  kv  siney  =  kt  sine^.  (2.76) 

Since  ku  equals  kv  (waves  in  the  same  medium),  the  first  equality  in 
(2.76)  expresses  Snell's  reflection  law:  angle  of  incidence  equals 
angle  of  reflection.  The  second  equality  in  (2.76)  is  Snell's  law  of 
refraction.  The  same  relationship  holds  at  the  z  =  0  boundary  and  can 
be  applied  co  isecutively  for  layered  media. 
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By  taking  advantage  of  (2.75),  equation  (2.73)  and  (2.74)  can 
be  simplified  to 


“jip  +  j 

E  u  e  cose,  +  E  v  e  cose,  =  E  *  cose* 
o  c  o  l  o  t 

_i  „  u  -j’l*  u  +j<|i  i  f 

Z2  tEo  •  -  E„  e  ]  =  Z3  E0 


(2.77) 

(2.78) 


where  eu  and  ey  have  been  replaced  by  e2  and  where  is  given  by 

*  -  (k“  d)  «  (k*  d)  =  kud  cose2  .  (2.79) 

A  matrix  equation  can  be  written  to  combine  (2.77)  and  (2.78), 


— —  — 

—ini 

—  — 

-J>  -j*  jlp 

e  +  e  Z  rn-0  ie  -  e  ) 

C  +  EoV>cos82 

r-  t 

Eq  coset 

-j*  j*  -j't' 

e  -  e  e  +  e 

zzK  -  Eov> 

Z"1  £  t 

J  fc0 

Z2cose2  2  2 

—  _ 

L 

— (2.80) — ! 

The  use  of  hyperbolic  functions  allows  this  to  be  written  as 


cosh(j>) 

•sinh(j^) 

Z2cose2 


-Z2co$e2sinh(j>) 
cosh  ( ) 


(Equ  ♦  E0v)cose2 


-  EoV) 


Eo  cos0t 


Z"1  E  1 
3  0 


M2. 81)  —l 

Multiplying  each  side  of  (2.81)  by  the  inverse  of  the  coefficient  matrix 
yields. 


(Equ  +  E0V)  cose2! 

Z2]  V  -  O 


f— 


[ cosh ( j4» )  Z2cose2  sinh(j^) 


.  i 


E„  cose* 
o  t 


Z'T  E  1 
3  0 


(2.82) 
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This  is  a  very  useful  result  since  it  relates  the  tangential  fields 
on  one  side  of  the  slab  to  those  on  the  other  side  or  in  general. 


EtN 

A 

B 

Et(N+l ) 

_HtN_ 

C 

0 

_Ht(N+l) 

(2.83) 


where  N  and  N+l  indicate  the  total  tangential  fields  at  the  Nth  and 
(N+l)th  interface  and  A,  B,  C,  and  D  are  identified  from  (2.82).  This 
result  is  analogous  to  the  ABCD  matrix  used  in  circuit  theory  to  relate 
the  input  voltage  and  current  of  a  two-port  network  to  the  output  volt¬ 
age  and  current. 

The  left-hand  side  of  (2.82)  contains  the  tangential  fields  at  the 
z  *  0+  boundary,  which  must  be  equal  to  the  tangential  components  of  the 
incident  fields.  Writing  this  in  matrix  form  yields. 


(Eo  +  cos01 

'(E“  *  E')  cose.,' 

.  (Eo  -  Eo >  . 

_Z-'  (E“  -  Ej)  _ 

(2.84) 


Since  and  E*"  are  both  of  interest  they  can  be  factored  out  of  the 
oo 

left-hand  side  of  (2.84)  and  the  right-hand  side  can  be  replaced  using 


(2.82) 


COS0i 

COS0^ 

1 

•r-  O 
Ui 

1 _ 

A 

B  . 

1 

m 

O  rt 

COS0 

-*;1 . 

t 

1 

u  o 

LU 

_ 1 

C 

D 

1 - 

r*i 

OJ  1 

Efc 

0 

(2.85) 


This  equation  contains  all  the  information  needed  to  find  the  transmis¬ 
sion  and  reflection  coefficients.  For  the  usual  case  of  a  slab  with 
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Transmission  and  reflection  coefficients  for  perpendicular  polari¬ 
zation  can  be  easily  obtained  from  the  above  by  duality.  The  symmetry 
of  the  electromagnetic  field  equations  permits  the  following  substitutions 
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to  be  made  in  any  equation  to  obtain  a  dual  equation:  H  replaces  E, 

-E  replaces  H  and  Z~^  replaces  Z.  These  and  other  duals  are  listed 
in  many  electromagnetic  texts  such  as  [36].  As  seen  in  Figure  2-13(a) 
and  (b),  these  substitutions  transform  parallel  polarization  to  perpen¬ 
dicular  polarization. 

The  ABCD  matrix  that  relates  tangential  field  components  at  the 
slab  boundaries  becomes. 


where. 


HtN 

A 

B 

Ht  (N+l ) 

.-Etf! 

C 

D 

-Et  (N+D_ 

• 

cose,  sinh(j>) 

A  B 

cosh(j>) 

C 

h 

l7  sinh(jtf)) 

C  0 

cos02 

cosh(jiji) 

(2.91) 


(2.92) 


and  ip  is  still  given  by  (2.79). 
The  dual  of  (2.85)  is 


cosei 

cosei 

H1 

0 

A 

B 

cose* 
o  t 

_  Z1 

-Z1  _ 

5. 

C 

D 

J  0 

(2.93) 


and  can  be  used  to  find  the  transmission  and  reflection  coefficients. 

The  perpendicular  polarization  transmission  coefficient,  T^,  is 
now  found  from. 


(2.94) 
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(c)  Perpendicular  polarization 


Figure  2-12.  Reference  Polarization 
Directions. 
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and  is  given  by. 


T,  = 


8  Z, 


C  cose 


A  + 


l 


cose^ 


+  D 


(2.95) 


The  perpendicular  polarization  reflection  coefficient,  r^,  is 
defined  as, 

jr 


r  .5 

1  «o  ’ 


(2.96) 


and  is  given  by  B  Z  C 


-A  - 


1 


cose 


cose.  Z1 
r±  B  Z,  C  cose. 


+  D 


A  + 


'I 


(2.97) 


cosei  Z1 


+  D 


The  minus  sign  in  (2.96)  is  chosen  to  conform  to  the  reference  direc¬ 
tions  of  Figure  2- 1 2 (c )  in  order  to  make  T^  =  T,  at  normal  incidence. 

The  insertion  phase  delay  for  perpendicular  polarization  is  found 
by  substituting  T^  for  T^  in  (2.90). 

Practical  radomes  are  usually  constructed  of  several  different 
materials  arranged  in  laminated  layers  in  order  to  achieve  design 
goals  for  strength,  weight,  rain  erosion  resistance,  etc.  The  above 
formulation  permits  the  quick  calculation  of  transmission  and  reflec¬ 
tion  coefficients  for  multiple  layered  slabs.  Recall  that  the  ABCD 
matrix  related  the  tangential  fields  at  opposite  sides  of  a  homoge¬ 
nous  layer.  When  several  layers  are  present  the  "output"  fields  of 
one  layer  become  the  "input"  fields  to  the  next  layer.  Thus  by 
cascading  the  transfer  matrices  for  each  layer  the  overall  transmission 
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and  reflection  properties  of  sandwiched  layers  can  be  determined. 

The  system  transfer  matrix  is  given  by 

(2.98) 

where  N  layers  are  present  and  the  terms  are  given  by  (2.82)  or  (2.92) 
depending  on  the  polarization.  Transmission  and  reflection  coeffi¬ 
cient  are  given  by  the  same  formulas  with  the  ABCD  parameters  deter¬ 
mined  from  (2.98).  Subroutine  TRANS  listed  in  the  appendix  computes 
the  transmission  and  reflection  coefficients  used  to  model  radome 
walls  in  this  study. 

Radomes  with  high  transmission  coefficients  are  desirable  in  most 
applications.  A  lossless  planar  slab  can  be  designed  to  have  a  trans¬ 
mission  coefficient  with  a  magnitude  of  one  for  a  particular  angle  of 
incidence.  This  is  seen  for  a  single-layer  slab  by  examining  the  denom¬ 
inators  of  the  transmission  coefficient  expressions  (2.87)  and  (2.95). 
When  the  slab  material  is  lossless,  i|>  is  purely  real  and  these  denomi¬ 
nators  will  equal  2  whenever  $  =  rnr,  with  n  an  integer.  For  a  slab 
constructed  of  a  given  material  this  is  achieved  by  making  the  thickness 
of  the  slab,  d,  according  to 
nx 

d  =  - 2 - -  ,  (2.99) 

2[urer  -  sin2ei] 

where  \Q  is  the  free-space  wavelength;  ur  and  £r  are  the  relative  perme¬ 
ability  and  permittivity  of  the  slab  respectively;  and  (2.76)  has  been 
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used  to  write  cos02  in  terms  of  sine^,  the  incident  angle.  For  normal 
incidence,  (2.99)  indicates  that  the  slab  thickness  must  be  an  integral 
number  of  half  wavelengths  in  the  slab  for  100%  transmission.  How¬ 
ever,  the  transmission  coefficient  for  this  case  is  unity  only  at 
normal  incidence  and  decreases  as  the  incidence  angle  increases.  Figure 
2-13  shows  the  transmission  coefficient  of  a  half-wave-thick  slab  as  a 
function  of  incidence  angle  for  both  polarizations.  Insertion  phase 
delay  is  also  plotted.  The  reflection  coefficients  for  this  slab  are 
shown  in  Figure  2-14. 

The  incidence  angles  for  the  radomes  in  this  study  are  quite  large, 
as  will  be  shown  in  the  next  chapter  and  a  slab  designed  for  complete 
transmission  at  an  angle  in  the  range  of  angles  actually  encountered 
will  make  a  better  radome  wall.  Figure  2-15  shows  the  transmission 
coefficients  for  a  slab  with  a  75°  design  angle  and  Figure  2-16  con¬ 
tains  the  reflection  coefficient  for  the  same  slab.  It  should  be 
noted  that  the  IPD  of  both  polarizations  are  equal  at  the  design 
angle  and  the  reflection  coefficient  vanishes  at  this  angle. 

PRINCIPLES  OF  RAY  TRACING 

This  section  discusses  the  application  of  ray  tracing  to  obtain 
electromagnetic  fields  that  are  transmitted  or  reflected  by  planar 
dielectric  slabs.  The  previous  section  dealt  with  infinite  planar 
slabs  and  tru»'  plane  waves  (infinite  in  extent).  Ray  tracing  fits 
this  theory  to  finite  structures. 

The  technique  of  ray  tracing  is  based  on  the  fact  that  all  electro¬ 
magnetic  fields  behave  locally  like  plane  waves  whenever  they  interact 
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Figure  2-13.  Transmission  Coefficient  of  Lossless 
Dielectric  Slab  with  Half-wavelength  Thickness. 
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Figure  2-14.  Reflection  Coefficient  of  Lossless 
Dielectric  Slab  with  Half-wavelength  Thickness. 
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Figure  2-15.  Transmission  Coefficient  of 
Lossless  Dielectric  Slab  with  75°  Design  Angle. 
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with  matter  that  is  large  with  respect  to  wavelength,  and  has  radii  of 
curvature  large  with  respect  to  wavelength.  Macroscopic  materials  have 
these  properties  at  optical  frequencies,  and  ray  tracing  is  frequently 
used  in  the  analysis  of  optical  systems. 

A  ray  is  a  directed  line  that  is  everywhere  perpendicular  to  a 
wave's  surfaces  of  constant  phase.  Rays  are  in  the  direction  of  wave 
propagation  which  is  the  direction  of  power  flow  or  the  Poynting  vector. 
When  dealing  with  plane  waves,  as  in  the  previous  section,  it  is  easy  to 
visualize  the  wave  traveling  in  the  direction  of  the  rays  in  the  various 
media. 

Snell's  transmission  and  reflection  laws  are  used  to  predict  the 
directions  of  the  "transmitted  ray"  and  "reflected  ray"  whenever  a  ray 
intersects  a  material  discontinuity.  Ray  tracing  of  a  transmitted  ray 
will  be  considered  here  first,  and  a  discussion  of  ray  reflection  will 
end  this  section. 

Consider  an  arbitrarily  polarized  wave  incident  on  a  planar  slab. 
This  wave  can  be  decomposed  into  a  sum  of  parallel  and  perpendicular 
polarized  waves  as  shown  in  Figure  2-17.  The  perpendicular  component  is 
found  by  first  forming  a  unit  vector  perpendicular  to  the  plane  of  inci¬ 
dence  (plane  containing  the  ray  and  the  surface  normal).  This  unit 
vector,  b,  is  often  called  a  binormal  vector  and  is  found  from, 

b  ,  (2.100) 

|rxn| 

A  A 

where  r  and  n  are  unit  vectors  in  the  directions  of  the  ray  and  surface 
normal  respectively  as  shown  in  Figure  2-17. 
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The  perpendicular  component,  E^,  of  the  Incident  electric  field, 

E^,  is  found  from  the  dot  product  as 

e{  -  V  •  b  .  (2.101) 

A 

Since  the  incident  field  is  transverse  to  r,  the  parallel  component  of 
the  incident  field  can  be  found  by  vectorial ly  subtracting  the  perpendic¬ 
ular  component  from  the  total  incident  field  or, 

E{,  -  E1  -  l[  ,  (2.102) 

where  E^ (  is  the  parallel  polarized  component  of  the  incident  field  and 
r]  is  given  by 

l[  =  (E1  •  b)b  .  (2.103) 

Of  course,  the  total  field  is  now  the  superposition  of  the  parallel 
and  perpendicular  parts  or, 

E1  =  I],  +  l[  .  (2.104) 

This  decomposition  is  useful  since  a  parallel  polarized  incident  wave 
produces  only  a  parallel  polarized  reflected  and  transmitted  wave.  Like¬ 
wise  there  is  no  cross-polarization  for  the  perpendicular  case. 

The  ray  angle  of  incidence,  9^,  is  important  since  the  transmission 
and  reflection  coefficients  are  functions  of  this  angle  as  shown  In  the 
preceding  section.  The  angle  of  incidence  can  be  found  from  the  vectors 

A  A 

r  and  n  by  use  of  the  dot  product,  or, 

1  A  A 

9i  =  cos"  (-r  •  n) 


(2.105) 
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The  ray  incident  on  the  slab  at  angle  ei  will  enter  the  slab  at 
angle  Og  given  by  Snell's  refraction  law,  (2.76),  and  emerge  on  the  op¬ 
posite  side  acain  at  angle  .  Although  the  transmitted  ray  exits  at 
the  same  angle  as  the  incident  ray,  there  is  an  offset  between  them 
because  1S  ^ess  t^an  9i  dielectric  slabs  in  air.  This  offset  is 

indicated  by  d  in  Figure  2-17.  The  amount  of  this  offset  will  depend  on 
the  slab  material,  slab  thickness,  and  the  angle  of  incidence.  Jenkins 
and  White  [37]  give  the  following  expression  for  d  produced  by  lossless 
dielectric  slabs, 

[n  cose.  I 

<2-'05> 

where  n  and  n^  are  the  indices  of  refraction  of  the  immersing  medium  and 
slab  respectively,  as  given  by  /u^  er  for  the  appropriate  medium. 

Plots  of  d  normalized  to  the  slab  thickness  versus  angle  of  inci¬ 
dence  are  shown  in  Figure  2-18  for  non-magnetic  slabs  with  various 
permittivities.  It  is  noted  that  for  intermediate  values  of  dielectric 
constant,  as  frequently  encountered  with  radome  materials,  the  deviation 
is  approximately  linear  with  incidence  angle. 

In  addition  to  the  ray  that  initially  emerges  from  the  dielectric 
slab,  other  rays  will  be  transmitted  after  multiple  reflection  from  the 
interior  slab  boundaries.  These  rays  exit  periodically  down  the  slab 
and  their  enerqy  diminishes  geometrically  with  the  number  of  interior 
reflections.  The  transmission  coefficients  previously  derived  consider 
all  these  reflections  in  order  to  satisfy  the  boundary  conditions  at  the 
interfaces. 


:  versus 
Various 
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The  approach  taken  in  this  work  is  to  ignore  the  ray  deviation  pro¬ 
duced  by  transmission  and  assume  the  incident  ray  continues  linearly 
through  the  slab.  The  radome  walls  under  consideration  here  have  thick¬ 
nesses  on  the  order  of  one-quarter  of  a  free-space  wavelength  and  even 
at  the  high  incidence  angles  encountered,  the  ray  deviation  is  approxi¬ 
mately  two-tenths  of  a  wavelength.  The  error  introduced  by  this  assump¬ 
tion  is  considered  to  be  in  line  with  other  approximations  made.  Multi¬ 
ple  reflections  are  taken  into  account  through  the  use  of  the  transmis¬ 
sion  coefficients  from  the  last  section. 

The  transmitted  electric  field  at  a  given  point  is  found  by  consid¬ 
ering  an  incident  ray  which  would  pass  through  the  point  with  the  slab 
removed  and  then  weighting  the  field  of  this  ray  with  the  appropriate 
transmission  coefficient,  or  in  equation  form, 

E1  •  E{|  T,,  +  e|t±  ,  (2.106) 

with  e| |  and  are  the  previously  derived  parallel  and  perpendicular 
components  of  the  incident  field,  f j (  and  fj_  are  the  insertion  transmis¬ 
sion  coefficients  given  by, 

T(|  *  |Tm1  exp  (-j  IPDj , )  ,  (2.107) 

and 

T±  --  |TJ  exp  (-j  IPD±).  (2.108) 

Since  the  two  transmission  coefficients  generally  have  magnitude 
less  than  one  and  differ  in  both  magnitude  and  phase,  equation  (2.106) 
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indicates  that  the  transmitted  field  is,  in  general,  attenuated,  phase 
shifted,  and  depolarized.  Furthermore  each  of  these  undesirable  charac¬ 
teristics  is  a  function  of  incident  angle.  All  these  phenomena  contrib¬ 
ute  to  radome-induced  radar  pointing  error. 

Reflection  may  also  be  treated  by  ray  tracing.  Consider  again  a 
plane  wave  incident  on  a  dielectric  slab  as  shown  in  Figure  2-19.  The 
incident  field  is  composed  of  parallel  and  perpendicular  components  as 
given  by  (2.102)  and  (2.103).  The  unit  vector,  r,  is  in  the  direction 
of  the  incident  ray  and  the  unit  vector,  p,  is  in  the  direction  of  the 
primary  reflected  ray. 

According  to  Snell's  reflection  law,  the  angle  of  incidence,  is 
equal  to  the  angle  of  reflection,  9r.  This  results  in  the  following 
equation  for  finding  the  reflected  ray  vector,  p: 

p  =  r  -  2(n  •  r)  n  ,  (2.109) 

A 

where  n  is  the  outward  unit  surface  normal. 

As  in  the  transmission  case,  additional  rays  emerge  from  the  slab 
due  to  multiple  internal  reflection  as  seen  in  Figure  2-19.  Again  for 
the  purpose  of  ray  tracing  these  rays  will  be  ignored;  however,  the 
reflection  coefficients  used  to  find  the  magnitude  and  phase  of  the  re¬ 
flected  ray  take  into  account  multiple  reflection  as  discussed  in  the 
previous  section. 

The  perpendicular-polarized  reflected  field,  Er,  is  given  by 

=  T±  e[  ,  (2.110) 
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where  is  the  reflection  coefficient  for  perpendicular  polarization, 
is  given  by  (2.103),  and  the  expression  is  evaluated  at  the  boundary. 

The  equation  for  the  parallel  polarizer!  reflected  field,  is  not 
as  simple  as  (2.110).  Recall  that  the  reference  direction  for  E((  was 
chosen  as  shown  in  Figure  2-19  so  that  parallel  and  perpendicular 
polarization  would  coincide  at  normal  incidence.  E1^  must  lie  in  the 
plane  of  incidence  and  be  normal  to  the  direction  of  propagation.  Thus, 
there  is  no  scalar  constant  which  gives  E1^  when  multiplied  by  | 

(except  at  normal  and  grazing  incidence). 

According  to  the  reference  directions  of  Figure  2-19,  the  tangential 
components  of  the  incident  and  reflected  field  are  related  by  the  reflec¬ 
tion  coefficient,  but  the  normal  components  are  related  by  the  negative  of 
the  transmission  coefficient.  The  normal  component  of  E| (  is  given  by 

ij  jn  *  (iu  •  n)  n  ,  (2.111) 


and  the  tangential  component  can  be  found  by  subtraction  as 


(2.112) 


where  n  and  t  in  the  subscript  indicates  normal  and  tangential  respective¬ 
ly.  The  reflected  components  at  the  boundary  are, 


r  e1 
ii  bnt 


(2.113) 


and 


(2.114) 


for  the  reference  directions  used  here. 
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The  total  reflected  field  is  the  sum  of  (2.113)  and  (2.114)  which 
simplifies  by  substitution  of  (2.111)  and  (2.112)  to, 

lru  -  r,,  [i],  -  2(e],  •  n)  n]  .  (2.115) 

This  equation  gives  the  magnitude,  phase,  and  direction  of  the  reflected 
field  at  the  interface. 

As  discussed  in  conjunction  with  the  transmission  problem,  an  inci¬ 
dent  field  will,  in  general,  have  parallel  and  perpendicular  field  compo¬ 
nents.  The  reflection  coefficients  for  these  two  polarization  states 
differ  in  magnitude  and  phase.  Thus,  as  was  the  case  with  transmission, 
the  fields  reflected  from  a  planar  dielectric  slab  will  be  attenuated, 
phase  shifted  and  depolarized.  One  might  hastily  conclude  that  these 
effects  must  degrade  radome  performance,  but  it  will  be  seen  that  indeed 
any  reflection  from  a  radome  wall  is  undesirable. 

This  chapter  has  discussed  the  basic  elements  used  in  radome  analy¬ 
sis.  The  following  chapter  will  connect  these  elements  to  analyze  radome 
system  performance. 
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III.  RADOME- INDUCED  RADAR 
POINTING  ERRORS 

In  this  chapter  the  effect  of  an  ogive  radome  on  an  incident  uni¬ 
form  plane  wave  is  considered.  It  will  be  shown  that  the  radome  causes 
amplitude,  phase  and  polarization  distortion  in  the  wave  transmitted 
through  the  radome.  This  distortion  causes  the  radar  inside  the  radome 
to  indicate  erroneous  pointing  directions. 

Two  mechanisms  for  wavefront  distortion  are  considered.  The  first 
is  due  to  the  variation  of  the  radome  wall  transmission  coefficient  with 
incidence  angle  and  the  second  is  due  to  reflection  of  energy  from  the 
interior  of  the  radome  wall.  The  methods  for  treating  these  sources  of 
distortion  are  based  on  techniques  presented  in  Chapter  II.  Radar 
response  to  wavefront  distortion  is  obtained  by  integration  of  the  inci¬ 
dent  wave  over  the  radar  antenna  aperture. 

This  chapter  begins  with  a  short  section  discussing  pointing  error. 
Sections  on  the  two  sources  of  wavefront  distortion  and  the  composite 
effects  of  distortion  follow.  The  chapter  is  concluded  by  considering 
radar  response  to  radome- induced  wavefront  distortion. 

CONCEPT  OF  POINTING  ERROR 

Radar  pointing  error  is  the  difference  between  the  angle  that  a 
radar  indicates  to  a  given  target  and  the  true  angle  to  that  target.  A 
similar  term  given  by  the  Institute  of  Electrical  and  Electronics 
Engineers  (IEEE)  is  boresight  error  which  is  defined  as  "the  angular 

deviation  of  xhe  electrical  boresight  of  an  antenna  from  its  reference 
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boresight"  [38].  Although  the  IEEE  gives  no  definition  for  the  word 
"boresight",  it  obviously  originated  from  the  practice  of  using  optical 
instruments  mounted  on  the  antenna  to  calibrate  the  electrical  output 
of  the  radar  with  the  mechanical  position  of  the  antenna. 

Sources  of  radar  pointing  error  may  be  broadly  divided  into  two 
categories:  internal  sources  and  external  or  environmental  sources. 
Internal  sources  of  error  refer  to  imperfections  in  the  radar  hardware 
(or  software)  and  would  include  amplitude  and  phase  unbalance  in  signal 
forming  networks,  noise  in  signal  channels,  nonlinearities,  and  quantiza¬ 
tion  errors.  External  sources  of  error  are  due  to  environmental  factors 
which  would  cause  even  a  "perfect"  radar  to  indicate  an  erroneous  direc¬ 
tion  to  the  radar  target.  Some  external  sources  of  error  are  glint, 
multipath,  non-uniform  propagation  media,  clutter,  jamming,  and  of  course 
radomes.  (A  different  point  of  view  will  be  taken  in  the  next  chapter 
where  the  radar  and  radome  are  considered  a  single  system.) 

As  a  simple  example  of  how  pointing  error  can  arise  consider  a 
uniform  plane  wave  incident  on  a  dielectric  prism  as  illustrated  in 
Figure  3-1 (a).  According  to  the  familiar  principles  of  the  prism,  the 
waves  entering  and  leaving  the  prism  are  in  different  directions  due  to 
the  unequal  optical  path  lengths  of  the  rays  through  different  parts  of 
the  prism.  A  radar  measuring  the  direction  to  the  source  of  the  wave 
leaving  the  prism  will  indicate  that  the  source  is  located  in  the  direc¬ 
tion  perpendicular  to  the  wavefront  at  the  antenna  as  was  discussed  in 
the  previous  chapter.  The  difference  between  the  direction  to  the  target 
indicated  by  the  radar  and  the  true  direction  to  the  target  is  the  point¬ 
ing  error. 
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(b).  Plane  wave  incident  on  radome. 


Figure  3-1.  Two  Illustrations  of  Wavefront  Distortion. 
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While  this  example  may  have  no  practical  counterpart,  it  does 
illustrate  an  environmental  source  of  pointing  error  that  behaves  much 
like  the  radone.  (This  example  is  also  similar  to  propagation  through 
a  stratified  atmosphere.) 

The  separation  of  colors  when  white  light  strikes  an  optical  prism 
illustrates  still  another  property  of  environmental  sources  of  pointing 
error.  Color  separation  is  due  to  variation  of  the  prism  dielectric 
constant  with  frequency.  This  causes  the  exit  angle  from  the  prism  to 
vary  with  frequency  (color).  For  white  light  incident  at  a  given  angle, 
the  angle  to  the  source,  as  indicated  by  an  optical  sensor  behind  the 
prism,  would  depend  on  the  frequency  to  which  the  sensor  is  sensitive. 
All  of  the  previously  mentioned  environmental  sources  of  pointing  error 
are  frequency  dependent  (some  quite  strongly). 

As  a  preview  to  the  remainder  of  the  chapter,  consider  a  uniform 
plane  wave  incident  on  a  radome  as  illustrated  in  Figure  3-1 (b) .  As  was 
seen  in  the  discussion  of  plane  waves  incident  on  planar  dielectric 
slabs,  the  wave  transmitted  through  the  radome  will  have  amplitude  and 
phase  that  varies  with  the  angle  between  the  radome  surface  normal  and 
the  incident  rays.  Variation  in  the  amplitude  of  the  transmission 
coefficient  causes  the  transmitted  wave  to  be  non-uniform.  Variation  in 
the  phase  (time  delay)  of  the  transmission  coefficient  causes  the  wave 
to  be  non-planar. 

Both  of  these  forms  of  distortion  lead  to  radar  pointing  error  in 
the  monopulse  radar  which  is  designed  to  operate  with  uniform  plane 
waves  incident  at  the  antenna.  Throughout  this  work  it  is  assumed 
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that  radome- induced  distortion  is  the  only  source  of  error  present;  that 
is,  all  radar  elements  and  other  propagation  considerations  are  assumed 
to  be  ideal. 

WAVEFRONT  DISTORTION  DUE  TO  TRANSMISSION 

In  the  receiving  formulation  of  the  radar  problem,  the  radar 
response  is  calculated  from  the  fields  incident  on  the  antenna  aperture 
and  known  characteristics  of  the  antenna.  A  radome  around  the  antenna 
introduces  distortion  into  the  wave  that  is  incident  on  the  antenna. 

The  first  distortion  mechanism  to  be  considered  here  treats  the  radome  as 
locally  planar  and  modifies  the  incident  wave  with  a  transmission  coeffi¬ 
cient  to  obtain  the  transmitted  wave.  Ray  tracing  is  then  carried  out 
to  account  for  propagation  to  the  antenna  aperture. 

The  radome  and  incident  field  are  depicted  in  Figure  3-2.  The 
incident  wave  is  regarded  as  being  made  up  of  parallel  rays  traveling 
perpendicular  to  the  wavefront.  The  radome  is  treated  as  being  locally 
flat  in  the  region  near  the  intersection  of  a  ray  and  the  radome  surface. 
This  planar  approximation  to  the  radome  surface  is  oriented  to  coin¬ 
cide  with  the  plane  tangent  to  the  radome  surface  at  the  ray-radome 
intersection  point,  and  the  thickness  of  the  planar  slab  approximation 
is  the  same  as  the  radome  wall  thickness  at  the  intersection  point. 

The  incidence  angle,  e^,  for  a  given  ray  is  determined  by  equation 
(2.105).  Parallel  and  perpendicular  compot ents  of  the  incident  field 
are  formed  with  respect  to  the  plane  of  incidence  as  outlined  in  Chapter 
II.  The  parallel  and  perpendicular  components  of  the  incident  field  are 
then  weighted  by  the  appropriate  plane  wave  transmission  coefficients 


1 


6 


Figure  3-2.  Tangent  Plane  Approximation  of  Radome  Wall. 
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for  the  planar  slab.  These  weighted  fields  now  make  up  the  field 
transmitted  to  the  radome  interior. 

The  ray  thus  transmitted  to  the  radome  interior  is  traced,  without 
deviation,  straight  to  the  antenna  aperture.  This  is,  of  course,  an 
approximation  as  was  discussed  in  the  previous  chapter.  The  field  over 
the  entire  antenna  aperture  is  obtained  by  considering  each  of  the  rays 
in  the  incident  wave  that  passes  through  the  aperture.  The  appropriate 
tangent  plane  approximation  of  the  radome  wall  must  be  determined  for 
each  ray,  since  the  direction  of  the  surface  normal  varies  over  the 
radome.  With  this  discription  of  the  ray  tracing  procedure  as  back¬ 
ground,  a  more  explicit  discription  of  the  details  is  now  given. 

The  tangent  ogive  has  been  previously  discussed,  and  the  same  con¬ 
ventions  and  notations  will  apply  to  the  present  discussion.  The  axis 
of  the  ogive  will  be  taken  as  the  y-axis.  The  y  coordinate  of  a  point 
on  the  ogive  surface  with  x  and  z  coordinates  given  by  x  and  z  is  found 
from, 

y  -  [R2  -  (V*2  +  z2  +  R  -  if)2]*3  ,  (3.1) 

V~2  2~ 

x  +  z  is  recognized  as  the  perpendicular  distance  from  the 
y-axis,  and  other  variables  are  as  defined  for  equation  (2.44). 

A  separate  coordinate  system  will  be  used  for  the  antenna  aperture. 
The  rectangular  coordinate  axes  of  this  system  are  designated  x',  y', 
and  z',  and  the  axis  of  the  aperture  is  taken  to  be  the  y'  axis.  The 
origin  of  the  x'y'z'  system  will  lie  on  the  ogive  axis.  Figure  3-3 
shows  an  outline  of  the  ogive  and  the  antenna  aperture  with  their 
respective  coordinate  axes. 
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The  origin  of  the  aperture  is  translated  with  respect  to  the  ogive 
origin  so  that  clearance  with  the  missile  body  is  achieved  for  large 
aperture  scan  angles.  The  symbol  y  is  used  for  this  translation. 

The  aperture  can  be  rotated  about  its  x'  and  z'  axes  to  direct  the 
antenna  in  any  direction  ir»  the  forward  half-space.  The  order  and 
directions  of  these  rotations  are  significant,  and  the  convention  adopted 
here  will  be  to  rotate  first  about  the  x'-axis  and  then  about  the 
z'-axis.  Positive  rotation  will  mean  that  the  rotation  of  the  primed 
coordinate  system  with  respect  to  the  unprimed  coordinate  system  is  in 
the  direction  that  the  fingers  of  the  right-hand  curl  when  the  right 
thumb  is  directed  in  the  positive  coordinate  direction  of  the  axis  of 
rotation.  The  symbol  a  will  be  used  to  denote  the  positive  rotation 
angle  about  the  x'-axis,  and  s  will  denote  the  positive  rotation  angle 
about  the  z'-axis.  The  angle  a  is  shown  in  Figure  3-3  and  e  is  zero. 

It  is  important  to  be  able  to  locate  points  with  primed  coordinates 
in  the  unprimed  coordinate  system.  The  transformation  that  accomplishes 
this  is  given  by. 


x 

y 

z 


COSB 

coset  sins 
sina  sins 


-sins 

COSa  COSS 

sina  coss 


0 

x' 

0 

sina 

y' 

4- 

yg 

COSa 

z' 

0 

(3.2) 


Without  loss  of  generality,  the  incident  rays  will  be  assumed  to  be 
parallel  to  the  yz-plane  since  the  coordinate  axes  of  the  radome  can 
always  be  rotated  to  achieve  this  orientation  due  to  the  symmetry  of 
the  radome.  A  rotation  of  the  antenna  aperture  about  its  y'-axis  would 
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be  necessary  to  consider  waves  incident  on  the  aperture  from  arbitrary 
directions,  however,  since  the  antenna  does  not  have  rotational  sym¬ 
metry  due  to  its  polarization  and  aperture  illumination  function. 

Points  on  an  incident  ray  obey  the  equation, 

y  ■  (z  -  z0)  tan  9  +  y0  ,  (3.3) 

where  y  and  z  are  coordinates  on  the  ray,  yQ  and  zq  are  the  y  and  z 
coordinates  respectively  of  some  point  on  the  ray,  and  0  is  the  angle 
of  the  ray  as  measured  from  the  xz-plane  (the  same  for  all  rays  since 
they  are  parallel)  as  shown  in  Figure  3-3  for  a  ray  passing  through  the 
aperture  center.  Of  course,  all  z-coordinates  on  any  particular  ray 
are  the  same. 

In  order  to  obtain  the  field  over  the  antenna  aperture,  rays  will 

be  traced  to  the  intersection  points  of  a  rectangular  grid  over  the 

aperture.  Suppose  (xfl,  y  ,  za)  is  a  point  in  the  aperture  expressed  in 

the  primed  coordinate  system.  The  coordinates  of  (x  y  ,  z  )  in  the 

unprimed  coordinate  system  are  obtained  by  applying  transformation  (3.2). 

The  equation  for  the  ray  passing  through  (x  y  z  )  is  given  by  (3.3) 

dad 

with  yQ  and  zQ  being  replaced  by  the  appropriate  transformed  coordinates 

<V  V  za>- 

Next  the  intersection  of  the  ray  and  the  radome  must  be  found.  The 
intersection  point  must  satisfy  the  equation  of  the  ray,  (3.3),  and  the 
equation  of  the  ogive,  (3.1).  Simultaneous  solution  of  these  two  equa¬ 
tions  results  in  finding  the  roots  of  a  (rather  messy)  fourth-order  equa¬ 
tion.  Since  the  x  coordinate  of  the  ray  is  constant,  the  intersection 
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point  must  have  this  x  coordinate.  Equation  (3.1)  thus  reduces  to  a 
function  of  one  variable  that  expresses  the  height  of  the  radome  along 
a  slice  parallel  to  the  yz-plane  as  a  function  of  z. 

The  intersection  point  is  the  point  where  the  ray  has  the  same 
height  as  the  ogive.  This  point  is  found  numerically  by  a  procedure 
based  on  Mueller's  method  [39]  for  finding  the  roots  of  functions  of  one 
variable.  Subroutine  INTRSC  listed  in  the  Appendix  performs  these 
computations. 

The  incidence  angle  of  the  ray  at  the  radome  is  found  from  a  unit 
vector  along  the  ray  and  the  unit  vector  normal  to  the  radome  at  the 
intersection  point.  The  unit  vector,  a,  pointing  out  of  the  radome  at 
the  intersection  point  and  along  any  incidence  ray  is  given  by, 

a  *  sine  y  +  cose  z  .  (3.4) 

A  method  for  finding  the  unit  vector  normal  to  the  radome  surface  has 
been  previously  discussed.  Subroutine  NORM  listed  in  the  Appendix  makes 
the  calculations  necessary  to  find  this  vector.  The  incident  angle  of 
the  ray  is  then  given  by, 

9i  =  cos"1  [a  •  n]  ,  (3.5) 

where  is  the  incident  angle,  and  n  is  the  unit  surface  normal  at 
the  intersection  point. 

The  polarization  of  the  incident  wave  is  transverse  to  the  direc- 
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tion  of  propagation  since  this  field  is  a  uniform  plane  wave.  Two 
principal  polarization  directions  will  be  considered  here.  These 


polarization  directions  are  defined  with  respect  to  the  yz-plane  and 
form  a  basis  to  express  an  arbitrarily  polarized  incident  field.  An 
incident  field  will  be  said  to  be  TE  polarized  (for  Transverse  Electric) 
if  its  electric  field  lies  transverse  to  the  yz-plane,  that  is  in  the 
x  direction,  for  any  incident  angle  9. 

An  incident  field  is  TM  polarized  (for  Transverse  Magnetic)  If  its 
electric  field  is  parallel  to  the  yz-plane  for  any  incident  direction. 
Unit  vectors  in  the  directions  of  TE  and  TM  polarized  fields  are  given 
by, 

aTE  =  *  *  (3.6) 

for  the  TE  field,  and 

aTM  *  cose  y  -  sine  z  ,  (3.7) 

for  the  TM  field,  where  9  is  the  previously  defined  angle  of  the 
incident  rays.  These  unit  vectors  are  indicated  in  Figure  3-3. 

Any  incident  polarization  can  be  expressed  as  a  sum  of  TE  and  TM 
components  with  appropriate  amplitude  and  phase.  Waves  polarized  in 
only  the  TE  or  TM  directions  will  first  be  considered.  Later  waves  with 
diagonal  linear  polarization  and  circular  polarization  are  discussed. 

The  procedure  for  obtaining  the  fields  transmitted  through  the 
radome  and  directly  incident  on  the  aperture  is  summarized  as  follows. 
Beginning  at  a  point  in  the  aperture  with  known  aperture  coordinates, 
the  radome  coordinates  of  the  point  are  found.  Next  a  ray  is  con¬ 
structed  through  the  aperture  point.  The  intersection  of  this  ray  with 
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the  radome  wall  is  found  numerically.  At  the  ray-radome  intersection 
the  radome  wall  is  approximated  by  a  planar  slab,  and  the  ray  angle  of 
incidence  is  *ound  from  the  dot  product  of  the  unit  surface  normal  and 
a  unit  vector  in  the  ray  direction. 

The  incident  field  at  the  ray-radome  intersection  is  decomposed 
into  parallel  and  perpendicular  components  with  respect  to  the  local 
plane  of  incidence  as  prescribed  by  equations  (2.101)  and  (2.102).  The 
transmitted  field  is  then  found  by  multiplying  the  incident  field  compo¬ 
nents  by  the  appropriate  plane-wave  transmission  coefficients  for  the 
planar  slab  approximation  to  the  radome  wall.  The  transmitted  field, 

E*,  is  given  by 

£*  =  t1  ^  +  f^E^  ,  (3.8) 

where  E^  and  E^  are  the  perpendicular  and  parallel  components,  respec¬ 
tively,  of  the  field  incident  at  the  aperture  point  in  the  absence  of 
the  radome;  and  and  are  the  insertion  transmission  coefficients 
for  perpendicular  and  parallel  polarization,  respectively,  as  given  by 
(2.108)  and  (2.107). 

It  has  been  previously  noted  that  and  T^  are,  generally, 
unequal  functions  of  incident  angle.  Because  the  radome  surface  is 
curved  the  local  ray  incident  angle  will,  in  general,  be  different  for 
any  two  different  rays.  Thus,  equation  (3.8)  indicates  that  the  phase 
and  amplitude  of  the  transmitted  field  will  vary  over  the  antenna  aper¬ 
ture.  Furthermore,  since  the  radome  surface  is  curved,  the  parallr“ 
and  perpendicular  components  of  the  incident  field  will  be  different  for 
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different  rays  in  the  aperture.  The  above  phenomena  are  summarized  by 
the  expression:  wavefront  distortion  due  to  transmission. 

An  example  is  presented  to  illustrate  the  ray-tracing  technique. 

The  radome  considered  is  a  tangent  ogive  with  length  of  30aq  and  a  base 
diameter  of  10aq  for  a  fineness  ratio  of  3.  The  radome  wall  is  taken  to 
be  a  lossless  non-magnetic  homogeneous  dielectric  with  relative  permit¬ 
tivity  of  4.  The  wall  has  a  constant  thickness  of  0.2855a  which  is  the 

o 

thickness  of  a  planar  slab  designed  for  100%  transmission  at  an  incidence 
angle  of  75°.  The  transmission  curves  for  this  dielectric  were  presented 
in  Figure  2-15.  A  design  angle  of  75°  was  selected  since  it  is  near  the 
mean  of  the  actual  ray  incidence  angles  of  a  wave  that  is  incident  along 
the  axis  of  an  ogive  with  fineness  ratio  of  3.  Most  of  the  rays  in  this 
case  are  incident  at  angles  ranging  from  70°  to  80°. 

Primary  emphasis  in  this  section  will  be  placed  on  the  phase  of  the 
field  over  the  antenna  aperture.  As  can  be  seen  from  the  transmission 
curves,  the  amplitude  of  the  transmission  coefficients  remain  near  unity 
over  the  entire  range  of  the  angles  encountered;  however,  the  phase  of 
the  transmission  coefficients  varies  about  20°  over  this  same  range. 

The  amplitude  of  the  field  over  the  aperture  is  thus  approximately  con¬ 
stant,  but  the  phase  is  not.  Since  the  monopulse  radar  indicates  the 
direction  to  the  target  as  perpendicular  to  the  weighted  phase  slope, 
the  phase  of  the  fields  incident  on  the  aperture  is  of  extreme  impor¬ 
tance. 

As  noted  earlier,  only  linear  polarization  in  the  TE  and  TM  direc¬ 
tions  will  be  considered  initially  and  the  phase  of  the  incident  field 
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across  the  antenna  aperture  will  be  found  for  a  given  polarization  and 
incidence  angle.  The  slope  of  the  phase-front  across  the  aperture  leads 
to  radar  pointing  errors  and  this  phase  front  distortion  varies  with 
incidence  angle.  Radar  pointing  errors  have  a  corresponding  angular 
dependence-  Thus,  considerable  insight  into  radome-induced  r3dar 
pointing  error  is  provided  by  the  phase  front  across  the  aperture. 

Figures  3-4  through  3-8  contain  the  insertion  phase  delay  for  TE 
polarized  rays  traced  through  the  radome  to  the  aperture,  plotted  as  a 
function  of  position  in  the  aperture.  This  series  of  figures  has  a 
progressively  lower  angle  of  incidence  (e)  with  the  incidence  angle 
noted  on  each  figure.  The  small  diagram  in  the  upper-right-hand 
corner  of  each  figure  is  similar  to  Figure  3-3  and  shows  the  incidence 
angle  for  a  ray  to  the  aperture  center  as  well  as  the  coordinate  axes. 

Also  indicated  on  the  small  figure  is  the  polarization  of  the  incident 
field. 

For  each  of  these  figures  the  aperture  was  rotated  to  be  normal  to 
the  incident  rays.  Thus,  each  figure  can  be  considered  as  the  phase  front 
of  the  distorted  wave  across  the  aperture  when  the  antenna  is  adjusted 
for  zero  pointing  error. 

The  aperture  used  for  these  figures  has  a  diameter  of  8x  and  the 
IPD  of  a  ray  is  plotted  above  the  base  plane  only  for  rays  incident  on 
the  aperture.  The  height  of  the  antenna  gimbal  point  above  the  radome 
base  (y  in  Figure  3-3)  is  2xq,  and  the  small  figure  shows  this  to  scale. 
The  grid  spacing  in  this  series  of  figures  is  xq/2. 

Several  attributes  of  these  figures  should  be  pointed  out.  First 
the  IPD  is  symmetric  about  the  x'  *  0  plane.  This  is  due  to  the  mirror 
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Figure  3-4.  Insertion  Phase  Delay  of  TE  Polarized  Field 
with  90°  Incidence  Angle. 


Figure  3-5.  Insertion  Phase  Delay  of  TE  Polarized  Field 
with  85°  Incidence  Angle. 
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synmetry  of  the  radome  and  incident  polarization  about  this  plane.  The 
shadow  of  the  radome  tip  is  observed  In  Figure  3-4  as  the  point  of  mini¬ 
mum  phase  delay  since  the  lowest  angle  of  incidence  occurs  for  the  ray 
passing  through  the  tip.  The  phase  front  is  observed  to  increase  from 
the  center  due  to  the  increase  in  the  angle  of  incidence  of  the  corre¬ 
sponding  rays  and  the  monotonic  Increase  in  the  IPD  with  incidence  angle 
for  this  radome  wall. 

In  Figure  3-5  the  shadow  of  the  radome  tip  is  seen  about  midway 
between  the  aperture  center  and  edge  as  the  point  of  minimum  phase  delay. 
At  80°  incidence  the  tip  shadow  has  moved  out  of  the  aperture,  but  the 
phase  front  is  still  observed  to  slope  down  toward  that  point.  At  the 
lower  angles  of  incidence  the  phase  front  is  seen  to  have  less  variation 
over  the  aperture  because  the  incident  rays  are  now  passing  through  the 
(relatively  flatter)  side  of  the  radome. 

Also  it  will  be  noted  that  as  the  incidence  angle  decreases  the 
average  height  of  the  IPD  across  the  aperture  decreases.  Again  this  is 
due  to  the  monotonic  variation  of  the  IPD  with  incidence  angle  over  the 
range  of  angles  encountered. 

Next  a  TM  polarized  incident  wave  is  considered.  Figures  3-9 
through  3-13  show  the  IPD  for  this  polarization  with  all  other  parameters 
as  in  Figures  3-4  through  3-8.  Again  the  IPD  is  noted  to  be  symmetric 
about  the  x'  =  0  plane  due  to  the  symmetry  of  the  incident  field  and  the 
radome.  At  on-axis  incidence  TE  and  TM  IPD  differ  only  by  a  90°  rota¬ 
tion  about  the  y  axis.  However,  at  non-axial  incidence  the  TE  and  TM 
IPD  are  different.  It  will  be  seen  later  that  radar  pointing  error  can 
be  a  strong  fjnction  of  polarization  at  certain  incidence  angles. 
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Figure  3-10.  Insertion  Phase  Delay  of  TM  Polarized 
with  85°  Incidence  Angle. 
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Figure  3-11.  Insertion  Phase  Delay  of  TM  Polarized  Field 
with  80°  Incidence  Angle. 


Figure  3-13.  Insertion  Phase  Delay  of  TM  Polarized  Field 
with  60°  Incidence  Angle. 
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Another  example  will  be  considered  to  illustrate  some  additional 
properties  of  wavefront  distortion  due  to  transmission.  In  the  previous 
example,  it  was  noted  that  the  insertion  phase  delay  of  the  planar  slab 
designed  for  100%  transmission  at  75°  is  a  monotonic  increasing  function 
of  incidence  angle  over  the  range  of  angles  encountered  (less  than  80°). 
This  feature  is  responsible  for  the  concave  phase  front  over  the  aperture 
at  normal  incidence.  The  effect  of  changing  the  wall  thickness  is  now  to 
be  considered. 

In  this  example  the  radome  wall  is  considered  to  be  polyimide 
quartz,  a  common  material  used  for  radome  construction  due  to  its  good 
mechanical  and  thermal  properties  along  with  good  electrical  character¬ 
istics.  Polyimide  quartz  is  non-magnetic,  has  a  relative  permittivity 
of  about  3.2,  and  a  loss  tangent  of  about  0.008  at  x-band.  A  wall 
thickness  was  selected  for  a  60°  design  angle  using  equation  (2.99) 
and  ignoring  losses.  The  transmission  coefficients  for  this  wall  are 
shown  in  Figure  3-14. 

The  first  thing  that  should  be  noted  about  these  transmission 
curves  is  that  the  magnitudes  are  not  unity  at  the  design  angle.  This  is 
due  to  the  non-zero  loss  tangent  of  the  material  which  prevents  100% 
transmission.  Next,  the  parallel  polarization  IPD  is  seen  to  be 
increasing  for  angles  up  to  about  82°,  but  the  perpendicular  IPD 
increases  for  angles  up  to  only  75°  and  then  begins  to  decrease.  This 
will  be  seen  to  have  serious  effects  on  phase  front  distortion  and  radar 
pointing  errors. 
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The  insertion  phase  delay  through  a  radome  constructed  of  this  new 
wall  material,  but  with  all  other  parameters  as  before,  is  presented  for 
the  incidence  angle  e  =  80°.  Figure  3-15  is  for  TE  polarization  and 
Figure  3-16  is  for  TM  polarization  plotted  as  before.  A  striking 
difference  now  exists  between  the  two  polarization  states.  The  TE  I PD 
now  has  an  average  negative  slope  in  the  z'  direction  but  the  slope  is 
positive  near  the  z'  axis.  The  slope  of  the  TM  I PD  is  positive  in  the 
z'  direction. 

Since  the  pointing  direction  indicated  by  the  radar  is  related  to 
the  slope  of  the  phase  front  across  the  aperture.  Figures  3-15  and  3-16 
show  that  a  reversal  in  the  direction  of  pointing  error  can  occur  at  a 
given  incidence  angle  due  to  a  change  in  the  polarization  of  the  incident 
field. 

The  large  differences  in  the  slopes  of  the  phase  fronts  above  are 
due  to  the  differences  in  the  IPD  of  the  planar  slab  for  parallel  and 
perpendicular  polarization.  For  minimum  phase  distortion  the  IPD  of  the 
two  polarizations  should  be  nearly  equal  over  the  entire  range  of  inci¬ 
dence  angles  encountered.  By  this  criterion  the  75°  design  wall  is 
better  than  the  60°  design  wall  for  the  tangent  ogive  radome  with  fine¬ 
ness  ratio  of  three.  Radar  pointing  error  curves  confirm  this. 

So  far,  incident  polarization  states  that  are  symmetrical  about  the 
radome' s  x'  *  0  plane  have  been  presented.  This  symmetry  is  not  present 
for  linear  polarization  in  any  but  the  pure  TE  or  TM  directions.  Figure 
3-17  illustrates  what  may  be  referred  to  as  diagonal  polarization  for  a 
ray  traced  to  the  aperture  center.  The  incident  polarization  in  this 
case  has  components  in  both  the  TE  and  TM  directions. 
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Figure  3-17.  Outline  of  Tangent  Ogive  Radome 
Showing  Both  Scan  Angles  and 
the  Polarization  Angle. 
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A  unit  vector,  aE,  in  the  direction  of  the  incident  field  can  be 
expressed  as 

A  A  A 

aE  *  cos(y)  aTE  +  sin  (y)  aTM  ,  (3.9) 

A 

where  y  is  the  polarization  angle  as  shown  in  Figure  3-17,  and  a^E  and 

A 

*TM  are  the  previously  defined  unit  vectors  in  the  TE  and  TM  directions 
respectively. 

Since  diagonal  polarization  is  not  symmetric  with  respect  to  the 
x'=  0  plane,  there  will  be  differences  in  the  fields  transmitted  to  the 
x'  >  0  and  x'  <  0  halves  of  the  antenna  aperture.  The  aperture  is  free 
to  rotate  about  its  z1  axis  to  track  the  incident  field.  (Figure  3-17 
shows  the  aperture  rotated  through  the  angle  b  positively  about  the  z* 
axis.)  Radar  response  to  diagonal  polarization  will  be  a  pointing  error 
in  the  B  rotation  direction  as  well  as  the  previously  mentioned  pointing 
error  in  the  a  rotation  direction.  These  pointing  errors  are  related  to 
the  phase-front  distortion  in  the  x'  and  z'  directions  respectively. 

To  illustrate  this  phase-front  distortion,  the  insertion  phase  delay 
over  the  antenna  aperture  is  plotted  for  the  polyimide  quartz  wall  for 
several  polarization  angles.  Figures  3-18,  3-19  and  3-20  show  the  phase 
front  for  polarization  angles  of  22.5°,  45°  and  67.5°  respectively. 

This  series  of  figures  have  the  same  parameters  as  were  used  in  Figures 
3-15  and  3-16.  The  polarization  angle  is  noted  on  the  small  inset  figure 
as  before. 

This  series  of  figures  shows  a  continuous  change  in  the  phase  front 
from  TE  polarization,  as  was  shown  in  Figure  3-15,  to  TM  polarization  as 
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Figure  3-20.  Insertion  Phase  Delay  of  67.5°  Polarized  Field  Incident 
at  80°  on  Polyimide  Quartz  Radome. 
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was  shown  in  Figure  3-16.  In  fact,  all  of  these  figures  may  be  taken 
together  to  represent  the  transition  from  0°  polarization  (pure  TE)  to 
90°  polarization  (pure  TM).  Note  also  that  the  maximum  average  phase 
slope  in  the  x'  direction  appears  at  a  polarization  angle  of  45°. 

The  phase-front  plots  presented  in  this  section  provide  considerable 
insight  into  radome- induced  radar  pointing  error.  They  will  be  referred 
to  when  pointing  errors  are  considered,  but  before  doing  this  another 
radome  distortion  mechanism  is  considered. 

WAVEFRONT  DISTORTION  OUE  TO  INTERNAL  REFLECTION 

This  section  concerns  the  analysis  of  the  electromagnetic  energy 
that  arrives  at  the  antenna  aperture  after  passing  through  the  radome 
and  being  reflected  from  the  interior  of  the  radome  wall  a  single  time. 
Ray  tracing  is  used  to  account  for  propagation,  and  the  radome  wall  is 
approximated  by  planar  dielectric  slabs  as  before. 

The  geometry  for  internal  reflection  is  illustrated  in  Figure  3-21. 
A  direct  ray  is  again  shown  passing  through  the  radome  wall  to  the 
aperture  center.  The  second  ray  shown  cannot  enter  the  aperture 
directly,  but  enters  the  aperture  after  being  reflected  from  the  inside 
of  the  radome  wall.  The  term  internal  reflection  will  be  used  to  refer 
to  illumination  of  the  aperture  by  rays  that  are  reflected  from  the 
interior  of  the  radome  wall. 

Internal  reflection  contributes  to  racar  pointing  error  and 


increased  antenna  sidelobe  levels,  because  the  internal  reflected  ener 
gy  does  not  arrive  at  the  aperture  from  the  direction  to  the  target. 
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Figure  3-21.  Outline  of  Tangent  Ogive  Radome 
Showing  Internal  Reflection. 
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The  total  field  over  the  aperture  is  the  superposition  of  the 
direct  and  reflected  fields.  For  this  reason  internal  reflection  is 
sometimes  referred  to  as  the  Lloyd's  mirror  effect  due  to  its  similarity 
with  an  optical  experiment  in  which  an  interference  pattern  is  obtained 
with  a  single  light  source  by  using  the  light  direct  from  the  source 
and  light  reflected  from  a  planar  mirror  [40]. 

It  has  also  been  pointed  out  that  the  internal  reflection  phenomena 
is  similar  to  holography  [41].  Recall  that  holography  involves  recording 
the  interference  pattern  of  a  reference  beam  and  an  object  beam.  In  the 
radome  case,  the  reference  beam  is  the  direct  wave  and  the  object  beam 
is  the  internal  reflected  wave.  The  two  waves  interfere  over  the  aper¬ 
ture  which  is  analogous  to  the  photographic  plate  used  in  optical  holo¬ 
graphy.  Use  of  this  concept  may  provide  a  powerful  tool  for  optimizing 
radome  performance. 

In  this  section,  only  the  internal  reflected  field  that  is  incident 
on  the  aperture  will  be  obtained.  The  total  field  in  the  aperture  and 
the  radar  response  to  internal  reflection  are  considered  in  the  follow¬ 
ing  sections.  The  beginning  point  is  to  obtain  the  field  that  is 
incident  on  the  inside  wall.  Note  that  internal  reflection  is  impossible 
for  on-axis  incidence.  Thus,  internal  reflection  need  not  be  considered 
until  the  incident  field  is  far  enough  off  axis  to  make  ray  reflection 
into  the  aper';ure  possible. 

In  this  study,  it  was  decided  to  consider  only  the  internal  reflect¬ 
ed  rays  that  enter  and  are  reflected  at  points  on  opposite  halves  of 


the  radome.  The  reason  for  this  can  be  seen  by  considering  a  ray  that 
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enters  and  reflects  from  the  same  side  of  the  wall  (ray  entry  and 
reflection  at  points  where  z  <  0  in  Figure  3-21).  Under  this  condition 
the  ray  incidence  angle  is  extremely  high  at  the  entry  point.  Since 
the  wall  transmission  coefficient  is  rapidly  going  to  zero  toward 
grazing  incidence,  little  energy  can  enter  the  radome  at  these  angles. 

It  must  be  mentioned,  however,  that  these  are  precisely  the 

■\ 

conditions  favorable  for  excitation  of  a  surface  wave  on  the  radome. 
Consideration  of  surface  waves  are  beyond  the  scope  of  this  analysis, 
but  some  of  their  possible  effects  will  be  mentioned  in  the  next  section. 

Under  the  conditions  of  opposite-wall  transmission  and  reflection, 
the  shadow  of  the  radome  tip  in  the  y  =  0  plane  lies  outside  the  radome. 
The  angle  at  which  this  occurs  is  used  in  the  computed  calculations  to 
signal  the  possibility  of  internal  reflection.  When  internal  reflec¬ 
tion  is  possible,  the  intersection  of  the  ray  and  radome  must  be  found 
at  two  points:  the  transmission  point  and  the  reflection  point. 

The- ray- radome  intersection  points  are  found  by  the  numerical 
search  procedure  previously  described.  The  transmission  point  is  found 
from  the  root  of  the  ray-radome  difference  equation  with  z  coordinate 
between  zero  and  the  radome  base  radius.  The  reflection  point  is 
associated  with  the  root  to  the  left  of  zero,  and  inside  the  radome  as 
shown  in  Figure  3-21.  The  radome  generating  parameters  for  the  radome 
interior  surface  are  used  in  finding  the  re 'lection  point. 

At  the  transmission  point,  the  transmitted  field  is  obtained  from 
the  incident  field  in  exactly  the  same  manner  as  for  a  direct  ray.  That 
is,  the  wall  is  modeled  as  a  planar  slab;  the  ray  incidence  angle  is 
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found  from  the  surface  normal;  the  field  is  decomposed  into  parallel 
and  perpendicular  components;  the  transmission  coefficients  are  found 
for  the  incident  angle;  and  the  transmittec  field  is  given  by  the  pro¬ 
duct  of  the  incident  field  and  the  transmiision  coefficient.  This 
transmitted  field  becomes  the  incident  field  at  the  reflection  point. 

At  the  reflection  point,  the  radome  wall  is  again  modeled  as  a 
planar  slab.  The  reflection  coefficient  and  direction  of  the  reflected 
ray  are  of  interest  at  this  point.  The  reflection  coefficient  is  a 
function  of  incidence  angle  and  polarization,  thus  these  must  be  known. 
The  incidence  angle  is  determined  in  the  same  manner  as  for  a  transmitted 
ray;  except  now  the  inward  unit  vector  normal  to  the  radome  wall  is  used. 

The  reflected  field  at  the  reflection  point  is  obtained  by  decompos¬ 
ing  the  incident  field  into  parallel  and  perpendicular  components  and 
then  multiplying  by  the  appropriate  plane  wave  reflection  coefficient. 
This  reflected  field  then  becomes  the  incident  field  at  the  aperture 
after  propagation  from  the  reflection  point. 

The  reflected  ray  leaves  the  reflection  point  in  the  direction 
predicted  by  Snell's  reflection  law.  Equation  (2.109)  gives  an  expres¬ 
sion  for  a  unit  vector  in  the  direction  of  the  reflected  ray  when  the 
vectors  that  apply  to  the  present  case  are  substituted. 

Not  all  reflected  rays  will  intersect  the  antenna  aperture.  Thus, 
it  is  importart  to  know  If  a  given  ray  will  be  incident  on  the  aperture. 
Indeed,  there  is  no  point  in  performing  the  above  calculations  for  rays 
that  do  not  contribute  to  the  radar  response  —  thus  this  condition  is 
checked  first.  The  method  used  is  to  ignore  all  rays  that 
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reflect  at  points  in  back  of  the  aperture  and  to  ignore  all  rays  that 
intersect  the  aperture  plane  at  points  outside  the  aperture  peri¬ 
meter. 

The  point  at  which  a  reflected  ray  intersects  the  aperture  is  im¬ 
portant,  because  it  determines  the  value  of  the  aperture  illumination 
function  to  be  applied  to  that  ray  during  aperture  integration.  The 
distance  from  the  ray  reflection  point  to  the  aperture  intersection 
point  is  also  important,  because  phase  retardation  occurs  over  this 
distance.  The  method  used  to  find  these  two  quantities  is  described 
next. 

Suppose  (xr,  y  ,  zr)  are  the  rectangular  coordinates  of  a  ray 

A 

reflection  point  and  c  is  a  unit  vector  in  the  direction  of  the  reflected 
ray  as  obtained  from  the  incident  ray  direction  and  the  radume  surface 
normal  by  equation  (2.109).  The  coordinates  of  a  point  (x,  y,  z), 
on  the  ray  trajectory  are  given  by. 


X  =  CS  +  X  , 
x  r 

(3.10) 

y  =  cys  +  yr. 

(3.11) 

2  =  C2S  +  Zr, 

(3.12) 

where  s  is  the  distance  from  (xp,  yr,  zr)  to  (x,  y,  z)  and  cx,  cy,  cz 
are  the  rectangular  components  of  c  which  are  simply  the  direction 
cosines  of  the  ray.  The  problem  is  to  find  s  such  that  (x,  y,  z)  lies  in 
the  antenna  aperture  plane. 
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The  equation  for  the  antenna  aperture  plane  can  be  obtained  in 
radome  coordinates  by  inverting  the  matrix  in  (3.2)  and  solving  for 
the  points,  y;  =  0,  since  these  points  describe  the  aperture  plane.  The 
y  coordinate  of  the  aperture  plane  is 


y 


tan  8 
cos  a 


tan  a  z  +  yg, 


(3.13) 


where  a  and  8  are  the  antenna  rotation  angles. 

The  distance,  s,  from  the  ray-reflection  point  to  the  ray-aperture 
plane  intersection  point  is  obtained  by  substituting  (3.10),  (3.11),  and 
(3.12)  into  (3.13)  and  solving  for  s.  This  yields 


s 


r  cos  a 


-c 


tan  8 
x  cos  a 


zr  tan  a  + 
+  cz  tan  a 


(3.14) 


The  ray-aperture  plane  intersection  point  is  obtained  from  (3.10),  (3.11), 
and  (3.12). 

Before  considering  some  examples,  the  procedure  for  finding  the 
internally  reflected  field  at  the  aperture  is  briefly  summarized.  The 
process  begins  by  tracing  the  incident  ray  through  the  radome.  The 
intersection  points  are  found  by  a  numerical  technique.  At  the  trans¬ 
mission  point  the  incident  ray  is  modified  by  the  transmission 
coefficient.  The  transmitted  ray  becomes  the  incident  ray  at  the 
reflection  poiit  on  the  opposite  side  of  the  radome.  At  the  reflection 
point  the  ray  is  modified  by  the  reflection  coefficient  and  redirected 
according  to  Snell's  reflection  law.  The  ray  then  travels  to  the 
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antenna  aperture.  The  phase  of  the  internal  reflected  field  at  the 
aperture  is  determined  by  the  transmission  and  reflection  coefficients 
as  well  as  the  propagation  distances  involved. 

A  radome  with  constant  wall  thickness  (designed  for  100%  trans¬ 
mission  at  75°  incidence)  will  again  be  used  for  an  example.  The 
transmission  and  reflection  coefficients  of  this  wall  were  presented 
in  Figures  2-15  and  2-16.  The  radome  is  a  tangent  ogive  with  the  same 
dimensions  (length  30xo,  base  diameter  1 0xo )  used  before,  and  the 
antenna  aperture  has  a  radius  of  4x0.  Incident  fields  with  TE  and  TM 
polarization  are  considered. 

In  Figure  3-22  the  magnitude  of  the  internal  reflected,  TE  polar¬ 
ized  field  over  the  antenna  aperture  is  shown  for  the  incident  angle, 

0  =  70°.  Recall  that  internal  reflection  does  not  occur  for  incidence 
angles  near  axial.  Seventy  degrees  (or  20°  off  axis)  represents  the 
beginning  of  favorable  conditions  for  internal  reflection. 

The  magnitude  of  the  incident  field  outside  the  radome  is  unity 
in  Figure  3-22  and  all  similar  figures  to  follow.  Thus,  the  neignt 
of  the  reflected  field  surface  provides  a  comparison  of  the  relative 
magnitudes  of  the  direct  and  reflected  fields. 

It  should  be  noted  in  Figure  3-22  that  the  reflected  field  is 
confined  to  the  left-rear  edge  of  the  aperture.  The  reason  for  this  is 
apparent  in  F.gure  3-21  where  it  is  seen  tt a t  favorable  conditions  for 
ray  reflection  exist  only  for  rays  to  the  left  edge  of  the  aperture  at 
high  incidence  angles.  As  the  incidence  angle  decreases,  the  reflected 
field  spreads  over  the  aperture. 
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Figures  3-23  through  3-25  show  the  reflected  field  for  incidence 
angles  beginning  at  60°  and  decreasing  in  10°  increments.  At  50° 
Incidence,  the  reflected  field  covers  the  entire  aperture.  At  incidence 
angles  less  than  40°  internal  reflection  is  no  longer  possible  because 
the  reflected  rays  either  pass  over  the  aperture  cr  intersect  the 
aperture  from  the  rear.  However,  multiple  internal  reflection  becomes 
a  possibility  at  these  low  incidence  angles. 

Figures  3-22  through  3-25  were  made  by  considering  a  set  of  rays 
spaced  \0/2  in  a  square  grid  perpendicular  to  the  direction  of  propaga¬ 
tion.  Each  ray  was  traced  through  the  radome,  reflected,  and  checked 
for  intersection  with  the  rotated  aperture.  If  a  ray  intersected  the 
aperture,  the  complex  value  of  its  vector  components  were  associated 
with  the  aperture  grid  point  nearest  the  intersection  point.  Whenever 
two  or  more  rays  intersected  the  aperture  in  the  same  grid  cell,  they 
were  coherently  sunmed  to  produce  the  total  field  in  that  cell.  Of 
course,  a  given  grid  cell  in  the  aperture  might  have  several  or  no  rays 
intersecting  it.  The  peaked  appearance  of  some  of  these  figures  is  a 
result  of  the  ray  spacing.  Use  of  a  finer  grid  would  tend  to  fill  the 
region  between  peaks. 

The  above  approach  is  like  considering  the  rays  as  centers  of  ray 
tubes  which  direct  the  flow  of  energy.  If  two  or  more  ray  tubes  inter¬ 
sect  at  the  aperture,  the  total  field  is  the  coherent  sum  of  the  fields 
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Figure  3-24.  Magnitude  of  Internal  Reflected,  TE  Polarized  Field 

with  50°  Incidence  Angle. 
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in  each  ray  tube.  The  total  field  magnitude  may  take  on  any  value  from 
zero  to  the  sum  of  the  individual  ray  magnitudes,  depending  on  the 
amplitude,  phase,  and  direction  of  the  individual  fields. 

The  left-rear  edge  of  the  aperture  is  seen  to  have  a  high  reflected 
field  magnitude  for  all  angles  of  incidence.  This  is  due  to  a  focusing 
effect  of  the  radome.  The  concave  shape  of  the  radome  surface  reflects 
rays  to  this  region  and  the  path  lengths  are  such  that  the  fields  there 
interfere  constructively. 

Reflected  fields  for  TM  polarization  are  considered  next.  Figures 
3-26  through  3-29  contain  these  data  presented  in  the  same  manner  as  the 
TE  case.  By  comparing  these  figures  with  the  corresponding  figures  for 
the  TE  case,  it  is  seen  that  the  spatial  distributions  of  these  reflected 
fields  are  the  same  (since  this  is  determined  by  geometry  alone),  but 
the  magnitudes  are  considerable  less. 

The  difference  in  the  TE  and  TM  reflected  field  magnitudes  can  be 
explained  by  the  wall  reflection  coefficients  and  the  polarization  of 
the  field  at  the  reflection  surface.  At  the  points  favorable  for  ray 
reflection,  a  TE  polarized  field  is  mostly  perpendicularly  polarized 
with  respect  to  the  wall;  while  a  TM  polarized  field  is  mostly  paral- 
lelly  polarized.  Upon  referring  to  the  wall  reflection  coefficients 
(Figure  2-6),  one  sees  that  the  magnitude  of  the  reflection  coefficient 
for  perpendicular  polarization  is  larger  tfan  that  for  parallel  polari¬ 
zation  at  the  high  incidence  angles.  This  points  out  the  need  for  a 
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wall  with  a  high  transmission  coeffient  (low  reflection  coefficient) 
at  the  ray  reflection  points  as  well  as  the  transmission  points. 

This  last  point  is  further  emphasized  by  considering  the  internal 
reflection  properties  of  the  polyimide  quartz  radome  introduced  in  the 
last  section.  Recall  that  this  wall  was  designed  for  100%  transmission 
at  60°  incidence,  but  the  design  goal  is  not  achieved  due  to  loss.  The 
reflection  coefficients  for  this  wall  are  given  in  Figure  3-30. 

The  magnitude  of  the  internal  reflected  field  over  the  aperture 
will  be  presented  as  in  the  last  example  with  no  changes  other  than 
the  wall  construction.  Figure  3-31  shows  the  magnitude  of  the  reflected 
field  for  TE  polarization  and  70°  incidence.  Note  that  a  different 
scale  is  used  for  the  magnitude  and  the  peak  is  nearly  ten  times  larger 

than  the  peak  in  Figure  3-22.  Recall  that  the  incident  field  outside 
the  radome  has  magnitude  of  only  one.  The  increased  reflection  coeffi¬ 
cient  of  the  wall  together  with  the  focusing  effect  have  produced  a 
field  with  magnitude  of  about  six  at  the  aperture  edge. 

The  internal  reflected  field  at  60°  incidence  is  shown  in  Figure 
3-32  again  for  TE  Polarization.  This  figure  also  shows  higher  field 
magnitudes  than  the  corresponding  75°  wall  case  (Figure  3-23). 

The  50  incidence  angle  case  is  shown  in  Figure  3^33.  Comparison 
with  Figure  3-24  shows  that  the  internal  reflections  are  less  than  for 
the  case  of  Figure  3-24.  This  is  because  the  incidence  angle  has  now 
decreased  to  values  where  the  wall  reflection  coefficient  is  small. 
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Figure  3-30.  Reflection  Coefficients  of  Polyimide  Quartz 
Wail  Oesigned  for  60°. 
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Figure  3-31.  Magnitude  of  Internal  Reflected,  TE  Polarized  Field 
Incident  at  70°  on  Polyimide  Quartz  Radome. 


Figure  3-33. 
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These  last  few  figures  indicate  the  reason  why  radome  wall  thick¬ 
ness  is  frequently  tapered  in  practice.  By  varying  the  wall  thickness, 
the  radome  can  be  made  highly  transparent  ii  the  region  that  fields 
enter  and  highly  transparent  in  the  areas  that  would  contribute  unwanted 
reflections. 

The  effect  of  internal  reflection  on  radar  pointing  error  will  be 
presented  after  the  next  section  which  concerns  the  composite  effects 
of  radome- induced  wavefront  distortion. 

COMPOSITE  WAVEFRONT  DISTORTION 

The  previous  two  sections  have  considered  two  factors  that  produce 
radome- induced  wavefront  distortion:  transmission  and  internal  reflec¬ 
tion.  The  operating  radome-radar  system  will  have  pointing  errors  due 
to  these  two  sources  and  several  others  all  acting  together  to  produce 
a  composite  effect. 

Primary  emphasis  in  this  study  is  placed  on  transmission  and 
internal  reflection;  however,  other  factors  may  produce  greater  pointing 
error  in  some  systems  or  at  certain  angles  in  any  system:  In  this 
section  some  of  the  composite  effects  of  transmission  and  internal 
reflection  are  discussed.  Also,  some  of  the  possible  sources  of  radome- 
induced  pointing  error  that  have  been  neglected  in  this  study  will  be 
mentioned. 

A  radome  produces  amplitude  and  phase  distortion  in  the  wave 
transmitted  directly  through  the  radome  to  the  antenna  aperture. 
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Furthermore,  for  scan  angles  sufficiently  far  off  axis  (greater  than 
about  14°  for  the  tangent  ogive  with  fineness  ratio  of  three),  the 
aperture  is  also  illuminated  by  internal  reflected  fields  that  do  not 
arrive  from  the  general  direction  to  the  radar  target.  Both  of  these 
sources  of  aperture  illumination  contribute  to  the  radar  response. 

The  principle  of  superposition  is  used  to  determine  the  total 
Incident  field  over  the  aperture.  Thus  the  total  field  over  the  aper¬ 
ture  is  the  coherent,  vector  sum  of  the  direct  and  reflected  fields. 

As  such,  the  relative  phase  and  direction  of  the  fields  are  important 
along  with  the  field  amplitudes. 

At  low  scan  angles,  the  internal  reflected  energy  is  incident  only 
at  the  edge  of  the  aperture  near  the  reflecting  surface.  This  energy 
arrives  in  phase  from  the  various  reflecting  points,  so  a  large 
reflected-field  amplitude  is  observed  at  the  aperture. 

At  large  scan  angles,  the  internal  reflected  field  spreads  over 
the  entire  aperture.  There  is  rapid  phase  variation  of  the  reflected 
field  over  the  aperture  in  this  case  because  the  arrival  direction  is 
far  off  the  antenna  axis.  Later  it  will  be  seen  that  the  radar  has 
little  response  to  these  fields,  since  their  sum  over  the  aperture  is 
small . 

In  addition  to  the  internal  reflected  field,  other  factors  contri¬ 
bute  to  the  ^otal  field  incident  at  an  aperture  in  a  real  radome.  One 
such  factor  is  tip  scattering.  Even  if  the  radome  wall  is  homogeneous, 
the  vertex  of  an  ogive  is  a  source  of  scattered  energy.  This  phenomena, 
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which  Is  much  like  edge  diffraction,  is  due  to  the  discontinuity  in 
surface  derivative  at  the  vertex.  Frequently  a  missile  radome  will 
contain  a  metal  or  ceramic  insert  in  the  tip  to  reduce  rain  errosion 
and  improve  strength.  Such  a  tip  insert  can  further  enhance  tip  scat¬ 
ter  because  of  the  material  discontinuity  at  the  tip.  Also  some  mis¬ 
sile  radomes  contain  a  Pitot  tube  at  the  tip  which  is  another  source 
of  unwanted  scatter.  Tne  effects  of  any  of  these  factors  are  diffi¬ 
cult  to  estimate.  Empirical  investigation  is  likely  to  provide  the 
greatest  rewards  in  this  area.  Indeed,  even  the  canonical  problem  of 
scatter  from  the  vertex  of  a  hollow  dielectric  cone  has  not  been  solved 
completely. 

Another  way  energy  may  arrive  at  the  antenna  aperture  is  via  a 
surface  wave.  A  surface  wave  is  a  field  that  is  bound  to,  and  propa¬ 
gates  along,  a  material's  surface.  Such  fields  attenuate  exponential¬ 
ly  away  from  the  surface.  Surface  waves  may  be  excited  on  a  body  by  a 
wave  that  is  incident  near  grazing,  such  as  is  often  the  case  for  a 
radome.  Since  the  surface  wave  attenuates  with  distance  away  from  the 
surface,  it  is  likely  that  any  effect  on  the  radar  system  would  be  due 
to  interaction  with  the  edge  of  the  antenna  near  the  radome  wall. 

Additional  scatterers  inside  the  radome  can  contribute  unwanted 
fields  to  the  receiver.  Such  scatterers  could  include  the  antenna 
feed  horn  and  supporting  structure,  tubing  to  Pitot  tube,  heating  wires 
for  a  radome  de-icer,  metal  behind  the  antenna,  etc.  Any  of  these 
items  may  contribute  to  distortion  of  the  received  wave. 
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Multiple  scattering  between  the  radome  internal  wall,  or  the  wall 
and  one  of  the  above  items,  or  between  any  two  items  mentioned  above,  or 
between  the  antenna  and  any  other  item,  or  any  combination  of  these, 
may  affect  the  received  wave  at  the  aperture. 

In  order  to  determine  the  radar  response  to  all  of  the  above 
factors,  and  others  not  mentioned,  it  would  be  necessary  to  solve  a  set 
of  integral  equations  for  the  fields  everywhere  in  and  around  the 
radome,  taking  into  account  all  of  the  matter  present  and  enforcing 
boundary  conditions  everywhere.  Such  equations  are  not  likely  to  be 
formulated  in  the  near  future  and  can  not  be  solved  for  practical  size 
radomes  because  of  the  limitations  of  even  the  largest  and  fastest 
computers.  For  this  reason,  high-frequency  approximations  have  proved 
most  useful  in  radome  analysis.  They  provide  considerable  insight  into 
the  physical  phenomena  involved,  and  they  will  likely  continue  to  be 
used  and  improved. 

Just  as  the  case  with  antenna  analysis  and  design,  engineering 
practice  and  simplifying  approximations  must  be  used  in  radome  analysis 
and  design  with  systematic  improvement  being  based  on  experiment. 
Transmission  aberrations  and  internal  reflection  are  considered  to  be 
the  greatest  source  of  radome-induced  radar  pointing  error,  but  it 
must  always  be  remembered  that  the  radar  response  is  the  composite 
effect  of  many  factors. 
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RADAR  RESPONSE  TO  WAVEFRONT  DISTORTION 

The  effect  of  a  tangent  ogive  radome  on  the  performance  of  a 
monopulse  radar  Is  considered  in  this  section.  Primary  emphasis  is 
given  to  the  most  undesirable  effect  of  the  radome— radar  pointing 
error.  Wavefront  distortion  produced  by  the  radome  has  been  discussed, 
therefore  this  section  begins  by  considering  how  monopulse  signals  are 
obtained  from  the  incident  field  at  the  antenna.  Throughout  this  sec¬ 
tion,  reference  is  made  to  the  aperture  phase  and  amplitude  plots  of 
earlier  sections  in  order  to  relite  the  phenomena  of  wavefront  distor¬ 
tion  to  radar  pointing  error.  The  effects  of  various  wall  thicknesses, 
aperture  illumination  functions,  and  polarizations  are  also  compared. 

In  Chapter  II  the  receiving  properties  of  an  aperture  antenna  were 
discussed.  It  was  assumed  that  a  uniform  plane  wave  was  incident  on 
the  aperture.  The  received  signals  were  determined  from  the  incident 
field  and  equivalent  currents  representing  the  antenna  by 

means  of  the  reciprocity  theorem.  The  same  procedure  will  now  be  used 
to  find  the  received  signals  when  the  antenna  is  enclosed  by  a  radome. 

As  explained  in  Chapter  II,  the  voltage  at  the  terminals  of  an 
aperture  antenna,  operating  in  the  receive  mode,  is 


V  =  (1  +  cos  e 


ds, 


(3.15) 


where  a  uniform  plane  wave  with  electric  field  intensity,  F1',  is  inci¬ 
dent  at  the  aperture;  e  is  the  angle  between  the  normal  to  the  antenna 
aperture  and  the  direction  to  the  source  of  the  incident  field;  Js  is 
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the  equivalent  electric  current  for  the  antenna  operating  in  the 
transmitting  mode;  C  is  a  complex  constant;  the  integration  is  over  the 
antenna  aperture;  and  V  is  either  the  sum  or  difference  voltage,  de¬ 
pending  on  whether  the  equivalent  currents  used  are  for  the  sum  or 
difference  pattern. 

When  radome- induced  wavefront  distortion  is  present,  the  incident 
field  at  the  aperture  is  not  planar  and  the  angle,  9,  in  (3.15)  becomes 
meaningless.  However,  since  the  distorted  wavefront  is  well  behaved 
over  the  aperture,  a  local  angle  of  incidence  can  be  found.  The 
incident  angle  can  then  be  brought  inside  the  integration  and  considered 
as  a  variable,  yielding 

v  =  c  f  (1  +  cos  e)!1  -X  ds.  (3.16) 

Jk  s 

In  order  to  obtain  numerical  results,  this  integral  is  approximated 
by  a  finite  sum  over  the  aperture  with  T1  determined  by  the  previously 
discussed  ray  tracing  procedures.  A  rectangular  grid  of  sample  points 
is  used  in  the  aperture  to  find  the  voltage  produced  by  the  direct 
wave.  This  yields 

N  N 

V  =  ^  +  cos  9mn^  *  Js(nAx',nAx')  as,  (3. 17) 

m=-N  n=-N 

with  m2  +  n2  <  N2 

where  9mn  is  the  incident  angle  at  the  mn- :h  sample,  the  functional 
dependence  of  E1  •  J  is  indicated,  as  is  ax'ax|  and  N  is  given  by 


N  -  [D/(2ax' )]  , 


(3.18) 


with  D  the  aperture  diameter,  ax' the  sample  point  spacing,  and  the 
brackets  take  the  largest  integer  less  than  or  equal  to  the  quantity 
inside. 

The  contribution  of  internal  reflected  energy  to  the  sum  and  differ¬ 
ence  signals  is  determined  from  the  reflected  rays  that  intersect  the 
antenna  aperture.  As  previously  described,  a  set  of  regularly  spaced 
rays  are  traced  through  the  radome,  reflected,  and  traced  to  the  aper¬ 
ture  plane.  Those  rays  actually  entering  the  aperture  contribute  to 
the  sum  and  difference  signals  according  to 

V  =  C  EE  +  cos  0mn)ri  *  As  •  (3*l9) 

with  (x')2  +  (z1)2  <  D/2 

where  x/and  z'are  the  aperture  coordinates  of  the  ray-aperture  intersec¬ 
tion  point,  emn  is  the  angle  between  the  antenna  normal  and  a  reflected 
ray,  and  the  summation  extends  over  all  possible  singly  reflected  rays. 

The  surface  current  density,  J$,  above  serves  the  dual  role  of 
sampling  the  tangential  component  of  the  incident  field  (through  the 
dot  product)  and  weighting  the  summation.  The  name,  aperture  illumina¬ 
tion  function,  was  given  to  the  magnitude  of  Js>  and  two  examples  “.were 
presented  in  Chapter  II:  uniform  and  cosine  taper.  Radar  pointing 
error  for  these  two  aperture  illumination  functions  will  be  compared 
in  the  examples  that  follow. 
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Throughout  this  section  the  receiving  antenna  is  assumed  to  be 
circular  polarized.  This  is  done  so  that  the  antenna  will  be  sensitive 
to  incident  fields  with  any  linear  polarization.  Actual  missile  radar 
systems  may  transmit  and  receive  either  linear  or  circular  polarization, 
but  the  field  scattered  by  the  target  will,  in  general,  have  time- 
varying  polarization  due  to  the  complex  geometry  of  the  target.  The 
assumption  of  a  circular-polarized  receiving  antenna  is  made  in  order  to 
compare  the  effects  of  various  incident  field  polarizations. 

A  circular-polarized  receiving  antenna  can  be  represented  by  an 
equivalent  surface  current  density  that  rotates  at  an  angular  velocity, 
u  (the  radian  frequency  of  the  fields),  in  the  aperture  plane.  Jordan 
and  Balmain  [42]  describe  how  a  rotating  unit  vector  in  the  direction  of 
this  current  can  be  constructed.  The  rotating  unit  vector  is  composed 
of  two  orthogonal  vectors  having  quadrature  time  dependence  with  left 
or  right  sense  determined  by  their  phase  relationship  (which  vector 
leads) . 

For  the  aperture  antenna  under  discussion  here,  two  suitable  ortho¬ 
gonal  vectors  in  the  aperture  plane  are  x1  and  z'  (see  Figure  3-3). 

These  vectors  are  combined  to  produce  the  rotating  unit  vector  in  the 
right-hand  direction,  ar,  as  given  by 

a  3  —  (x1  -  jz1),  (3.20) 

r  y? 

where  propagation  in  the  -y'  direction  (receiving)  has  been  assumed. 

The  square  root  of  two  factor  is  necessary  to  normalize  ar>  that  is 
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make  ar  •  a  ■  1.  Recall  that  one  vector  is  conjugated  when  taking  dot 
products  of  complex  vectors.  (The  "unit"  vectors  given  by  Jordan 
and  Balmain  have  not  been  normalized.) 

The  left-hand  unit  vector  is  simply  the  conjugate  of  (3.20).  It 
is  not  difficult  to  show  that  left  and  right  circular  polarizations  are 
orthogonal.  Left-circular  polarization  will  not  be  considered  here. 

The  equivalent  surface  current  density,  7S,  for  the  aperture 
antenna  can  now  be  written  as 

Js  »  g(r,$)  ar,  (3.21) 

where  g  is  the  (scalar)  aperture  illumination  function,  either  uniform 
or  cosine  taper  here. 

A  computer  program  was  written  to  perform  the  numerical  aperture 
integration  described  in  this  section.  The  fields  at  the  antenna  are 
calculated  by  ray  tracing.  The  radome  wall  transmission  and  reflection 
coefficients  are  calculated  and  stored  at  3°  intervals  before  any 
aperture  integrations  are  performed.  Whenever  a  transmission  or 
reflection  coefficient  is  required,  the  three  closest  stored  values 
are  interpolated  with  a  second-order  polynomial  to  provide  the  desired 

coefficient.  This  interpolation  scheme  was  found  to  provide  better 
than  five  decimal-place  accuracy,  which  is  well  within  the  accuracy  of 
other  approximations  used. 

The  complete  program  used  for  radome  error  calculation  is  listed 


in  the  Appendix. 
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This  program  has  two  options  after  monopulse  signals  have  been 
computed.  The  first  is  to  reposition  the  antenna  to  the  angles  indi¬ 
cated  by  the  computed  monopulse  signals.  Indicated  errors  are  then 
recomputed  and  the  antenna  repositioned  according  to  the  new  indicated 
angles.  This  procedure  is  repeated  until  the  antenna  is  repositioned 
less  than  a  user-specified  amount.  Radar  pointing  error  is  then  the 
difference  between  the  pointing  angle  of  the  antenna  and  the  angle  of 
the  incident  field. 

The  second  option  is  to  define  radar  pointing  error  as  the  angle 
indicated  by  the  monopulse  signal  when  the  antenna  has  no  pointing 
error.  This  option  was  included  in  the  program  when  it  was  discovered 
that  after  the  initial  repositioning  of  the  antenna  according  to  the 
first  option,  the  next  indicated  errors  were  consistently  zero  within 
numerical  accuracy.  Thus,  both  options  provide  the  same  results,  but 
the  second  requires  only  about  half  the  computer  run  time.  The  pointing 
errors  .presented  throughout  this  work  are  the  indicated  errors  according 
to  the  second  option. 

Some  examples  of  radar  pointing  error  are  now  considered.  For 
each  of  these  examples,  the  radome  is  considered  to  be  a  tangent  ogive 
with  constant  wall  thickness  and  fineness  ratio  of  three.  The  radome 
base  diameter  is  1 0ao  and  the  antenna  aperture  diameter  is  8\0.  The 
antenna  gimbal  point  is  2\0  above  the  radome  base  plane  (y  =  2x0  in 
Figure  3-3).  All  these  conditions  are  the  same  as  for  previously 
presented  examples. 
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The  antenna  aperture  can  be  rotated  in  two  orthogonal  angular 
directions:  the  a  and  0  scan  angles  of  Figure  3-17.  Radar  pointing 
error  for  a  given  scan  direction  is  defined  as  the  scan  angle  of  the 
antenna  minus  the  scan  angle  of  the  incident  wave.  Pointing  errors  in 
the  a  scan-angle  direction  are  called  in-plane  pointing  errors  since 
they  are  in  the  plane  formed  by  the  incident  field  direction  and  the 
radome  axis.  Pointing  errors  in  the  6  scan-angle  direction  are  called 
cross-plane  pointing  errors. 

A  radome  with  the  wall  design  for  100%  transmission  at  75°  is 
considered  first.  The  transmission  and  reflection  coefficients  for 
this  wall  were  given  in  Figures  2-15  and  2-16.  The  incident  field  is 
polarized  in  the  principal  polarization  directions,  TE  and  TM,  but 
the  right-hand  circular  component  of  these  fields  is  detected  at  the 
aperture. 

The  in-plane  pointing  error  for  this  case  is  shown  in  Figure  3-34. 
The  aperture  illuminiation  function  is  uniform.  Pointing  error  is  shown 
as  a  function  of  the  a  scan  angle  and  the  effect  of  including  internal 
reflection  is  illustrated.  Note  that  internal  reflection  has  no  effect 
on  pointing  error  for  scan  angles  up  to  about  14°  since  internal 
reflection  is  not  possible  for  these  angles. 

In  the  TM  case,  the  internal  reflected  fields  cause  almost  no 
change  in  pointing  error.  But  in  the  TE  case,  internal  reflection 
<  causes  a  oscillation  in  the  pointing  error  about  the  no-reflection 
value.  As  was  seen  in  Figures  3-22  through  3-29,  TE  polarization 
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produces  much  larger  reflections  to  the  aperture  than  TM  polarization. 
The  perturbation  of  the  TE  pointing  error  is  thus  related  to  the  inter¬ 
ference  between  this  larger  reflection  and  the  direct  wave. 

It  is  also  noted  in  Figure  3-34  that  the  pointing  errors  are  nega¬ 
tive  for  all  scan  angles  (except  TE  at  0°  to  6°).  This  is  due  to  the 
slope  of  the  phase  front  of  the  incident  field  at  the  aperture.  As 
was  shown  in  Figure  3-5  through  3-8  and  3-10  through  3-13,  the  average 
slope  of  the  phase  front  over  the  aperture  is  positive  in  the  z'  direc¬ 
tion  when  the  aperture  has  zero  pointing  error.  This  means  that  the 
aperture  must  scan  in  the  negative  o  direction  in  order  to  become 
aligned  with  the  average  phase  slope,  so  the  radar  pointing  error  is 
negative. 

Next,  consideration  is  given  to  the  effects  of  substituting  a 
different  wall  design  in  the  above  example.  The  wall  design  for  60° 
incidence  previously  presented  will  be  used.  The  transmission  and 
reflection  coefficients  for  this  wall  were  given  in  Figures  3-14  and 
3-30.  In-plane  radar  pointing  error  for  this  wall  design,  but  with  all 
other  parameters  the  same  as  for  Figure  3-34,  are  presented  in  Figure 
3-35. 

Note  the  differences.  Beginning  with  zero  pointing  error  on 
axis,  the  TE  error  curve  increases  while  the  TM  error  curve  decreases. 
Maximum  error  is  larger  and  occurs  at  a  lower  scan  angle  than  for  the 
75°  wall  design.  It  was  shown  in  Figures  3-15  and  3-16  that  for  the  60c 
wall  design  and  10°  scan  angle,  the  average  phase  slope  across  the 
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aperture  in  the  z*  direction  had  different  signs  for  TE  and  TM  polariza¬ 
tion.  These  different  phase  slopes  lead  to  radar  pointing  errors  in 
different  directions.  This  could  have  serious  implications  for  the 
tracking  radar,  since  the  indicated  direction  to  the  target  could 
change  by  a  large  amount  simply  because  of  a  change  in  the  incident 
polarization. 

The  effect  of  internal  reflection  in  Figure  3-35  again  begins 
at  a  scan  angle  of  about  12°.  But  now  instead  of  causing  a  small  per¬ 
turbation  on  the  no-reflection  error  curve,  internal  reflection  causes 
large  deviations  from  this  curve.  TE  polarization  is  again  seen  to  have 
the  greater  sensitivity  to  internal  reflection  effects.  The  large 
reflected  fields  that  illuminate  the  aperture,  as  seen  in  Figures  3-31, 
are  responsible  for  these  large  errors. 

At  large  scan  angles,  say  greater  than  45°,  internal  reflection 
produces  little  effect  on  radar  pointing  error.  As  previously  discussed, 
the  reflection  surface  becomes  highly  transparent  at  the  larger  scan 
angles,  thus  less  reflected  energy  is  entering  the  aperture.  Also  at  the 
larger  scan  angles,  the  reflected  field  has  several  cycles  of  phase 
variation  over  the  aperture;  thus,  the  integrated  response  to  internal 
reflection  is  smaller.  This  last  effect  is  similar  to  the  response  to 
a  plane  wave  that  enters  the  aperture  through  a  side-lobe  of  the  antenna 
pattern. 

The  effect  of  the  aperture  illumination  function  on  radome- induced 
pointing  error  is  considered  next.  In  Chapter  II  it  was  noted  that 
uniform  aperture  illumination  produces  a  far-field  antenna  pattern 
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with  a  narrower  main  lobe  than  cosine  tapered  illumination,  but  higher 
sidelobe  levels  are  associated  with  uniform  illumination.  When  receiv¬ 
ing,  the  aperture  illumination  function  serves  to  weight  the  incident 
fields  during  aperture  integration. 

In-plane  radar  pointing  error  for  the  two  principle  polarizations 
and  both  aperture  illuminations  are  shown  in  Figure  3-36  for  the  75° 
wall  design  and  in  Figure  3-37  for  the  60°  wall  design.  Internal 
reflection  is  included  for  both  figures.  In  each  figure  it  is  seen  that 
cosine  taper  reduces  the  effect  of  internal  reflection.  This  is 
because  cosine  taper  reduces  the  contribution  of  fields  near  the  aper¬ 
ture  edge  to  the  aperture  integration. 

However,  at  low  scan  angles,  the  radar  pointing  error  is  greater 
for  the  cosine-tapered  illumination.  This  is  thought  to  be  due  to  the 
effect  of  fields  near  the  shadow  of  the  radome  tip.  In  this  shadow 
region  the  field  is  highly  distorted  because  of  the  large  variation  of 
incidence  angle  in  the  tip  region.  At  low  scan  angles,  the  tip  shadow 
is  near  the  aperture  center  (the  region  given  highest  weight  by  cosine 
taper),  thus  the  aperture  is  responding  to  the  most  distorted  part  of 
the  incident  field. 

So  far,  only  polarization  in  the  pure  TE  and  TM  directions  has  been 
considered  for  the  incident  field.  Radar  response  to  a  diagonal,  linear¬ 
ly-polarized  incident  field  is  now  considered.  Such  a  field  has  both  TE 

and  TM  components  as  given  by  equation  (3.9).  The  polarization  angle, 
y,  indicated  in  Figure  3-17  will  be  used  as  a  parameter.  Diagonal 
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Figure  3-36.  In-Plane  Radar  Pointing  Error  for  Radome  with  75°  Wall  Design 
Illustrating  Effect  of  Aperture  Illumination  Function. 
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polarization  leads  to  both  in-plane  and  cross-plane  radar  pointing 
errors,  as  would  be  expected  from  the  phase-front  plots  presented 
earlier. 

In-plane  radar  pointing  error  versus  scan  angle,  a,  is  shown  in 
Figure  3-38  for  the  75°  wall  design  and  in  Figure  3-39  for  the  60° 
wall  design.  Error  curves  are  plotted  for  the  polarization  angles 
indicated.  Uniform  aperture  illumination  is  used  for  both  figures;  and 
the  effect  of  internal  reflection  is  included.  In  Figure  3-38,  it  is 
seen  that  little  variation  in  pointing  error  occurs  as  the  polarization 
is  changed.  This  is  due  to  the  near-equal  insertion  phase  delay  of  the 
75°  wall  design  for  parallel  and  perpendicular  polarization.  Thus,  the 
phase  delay  in  passing  through  the  radome  wall  will  be  about  the  same 
whatever  the  incident  polarization. 

Something  very  different  happens  with  the  60°  wall  design.  The 
opposite-sign  errors  for  TE  polarization  {y  ~  0°)  and  TM  polarization 
[y  *  90°)  have  been  pointed  out.  Pointing  error  is  now  observed  to  take 
on  intermediate  values  between  these  two  as  the  polarization  angle 
changes.  The  phase  front  across  the  aperture  for  a  =  10°  was  shown  in 
Figures  3-15  through  3-20  for  each  of  the  linear  polarization  angles 
considered.  It  is  interesting  to  see  how  the  phase  slopes  relate  to 
the  radar  pointing  errors. 

In  addition  to  an  average  slope  in  the  z'  direction  across  the 
aperture,  diagonal  polarization  leads  to  an  average  slope  in  the  x' 
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direction.  This  average  phase  slope  in  the  x'  direction  produces  radar 
pointing  error  in  the  8  scan  angle  or  cross-plane  pointing  error.  Mono¬ 
pulse  signals  are  formed  for  the  8  scan  direction  in  the  same  manner  as 
for  the  a  scan  direction.  The  tracking  radar  measures  the  angle  to  the 
target  in  both  the  a  and  8  directions  simultaneously  and  repositions  the 
antenna  to  obtain  zero  indicated  error  signals  in  both  angles. 

The  antenna  has  zero  cross-plane  error  when  e  -  0  according  to  the 
convention  adopted  for  the  incident  field  (see  Figure  3-17).  The  cross¬ 
plane  pointing  error  is  thus  taken  to  be  the  8  scan  angle  of  the 
antenna  when  the  cross-plane  monopulse  signal  is  zero. 

Cross-plane  pointing  error  versus  a  scan  angle  is  presented  for 
the  two  radomes  previously  used  for  examples.  Figure  3-40  is  for  the 

75°  wall  design  and  Figure  3-41  is  for  the  60°  wall  design.  Note  that 

both  radomes  have  zero  error  at  a  =  0,  but  the  cross-plane  error  for  the 
75°  wall  design  rises  slow  and  almost  monotonically,  while  that  of  the 

60°  wall  design  rises  to  a  peak  at  about  a  =  10°  and  then  decreases  to 

small  errors  past  a  =  30°. 

At  first  thought,  it  would  seem  that  Figures  3-40  and  3-41  contain 
mistakes  since  the  cross-plane  error  is  not  zero  for  0°  and  90°  polariza¬ 
tion,  the  symmetrical  polarization  states.  The  aperture  phase  front 
plots  did  show  these  fields  as  being  symmetrical  about  x'  *  0.  But  that 
which  was  plotted  in  these  figures  was  the  component  of  the  field  at  the 
aperture  with  the  same  polarization  as  the  incident  field.  Radome- 
induced  distortion  has  mirror  symmetry  about  the  x  *  0  plane.  A  TE  or 
TM  polarized  antenna  also  has  symmetry  about  the  %’  =  0  plane.  However, 
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a  TE  polarized  antenna  would  have  no  cross-plane  error  at  any  scan  angle 
when  a  TE  polarized  incident  field  is  present;  likewise  a  TM  antenna  has 
no  cross-piano  error  when  a  TM  field  is  incident. 

The  errors  plotted  in  Figures  3-40  and  3-41  are  for  a  right-hand 
circular-polarized  receiving  antenna.  The  circular-polarized  anten¬ 
na  does  not  have  mirror  symmetry  about  the  x‘  =  0  plane.  Thus  the 
circular-polarized  antenna  does  not  have  zero  pointing  error  even 
when  symmetrically  distorted  fields  are  present. 

A  circular-polarized  antenna  has  been  used  here  so  that  the  same 
antenna  will  receive  all  incident  linear  polarizations.  If  a  1  inear- 
polarized  antenna  were  used,  the  received  signal  level  would  drop  sub¬ 
stantially  for  cross-polarized  incident  fields.  The  radar  pointing 
error  due  to  a  poor  signal -to-noise  ratio  is  likely  to  be  much  greater 
than  the  radome- induced  pointing  errors. 

Pointing  error  curves  for  incident  fields  having  the  same  polariza¬ 
tion  as  the  receiving  antenna  (RHC)  are  shown  in  Figures  3-38  through 
3-41.  For  this  case,  no  cross-plane  error  is  present  (again  due  to 
symmetry  considerations).  Both  Figures  3-40  and  3-41  show  that  the 
pointing  error  for  0°  polarization  and  90°  polarization  are  almost 
equal  and  opposite.  Since  circular  polarization  is  made  up  of  equal 
parts  of  these  two  polarizations,  the  net  pointing  error  is  zero. 

It  is  also  noted  in  the  cross-plane  error  curves  that  maximum 
error  is  for  45°  polarization  at  all  scan  angles.  This  is  the  polariza¬ 
tion  angle  with  maximum  asymmetry. 
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In  this  chapter  some  fundamental  concepts  have  been  applied  to 
obtain  approximations  for  radome- induced  distortion.  Although  the 
approach  taker  has  been  quite  simple,  it  provides  good  insight  into 
radome  behavior.  Throughout  the  chapter  the  relation  of  wavefront 
distortion  to  pointing  error  was  stressed.  Several  of  the  factors 
necessary  for  a  good  radome  have  also  been  pointed  out.  Next,  a  tech¬ 
nique  to  compensate  for  pointing  error  and  some  additional  information 
available  from  monopulse  are  considered. 


IV.  RADOME- INDUCED  RADAR  POINTING  ERROR  CORRECTION 

In  this  last  major  chapter,  an  error  compensation  technique  is 
presented.  The  technique  is  easily  implemented  and  should  be  quite 
effective  in  the  digital-computer-controlled  radar  systems  of  today  and 
the  future.  The  second  section  of  the  chapter  treats  the  relation  of 
the  complex  indicated  angle  and  radome-induced  wavefront  distortion. 

The  relation  of  the  complex  indicated  angle  to  some  other  sources  of 
distortion  has  been  presented  in  the  literature,  but  treatment  of  the 
radome  problem  is  thought  to  be  new. 

ERROR  COMPENSATION 

The  concept  of  error  compensation  is  very  simple:  if  the  error 
characteristics  of  a  system  are  known  a  priori  as  a  function  of  system 
parameters,  then  the  error  can  be  removed  from  the  measured  output  to 
yield  an  error-free  measurement.  Of  course,  a  sufficiently  simple 
system  is  required  so  that  each  output  can  be  uniquely  associated  with 
a  given  input.  This  is  the  case  for  radome-induced  radar  pointing 
errors. 

The  traditional  approach  to  radome-induced  error  reduction  is  trial 
and  error.  Beginning  with  a  wall  design  optimized  for  other  system 
considerations  (strength,  errosion  resistatce,  thermal  properties,  etc.), 
the  pointing  errors  are  then  measured.  Next  the  radome  is  "patched-up" 
by  applying  layers  of  dielectric  tape  to  parts  of  the  radome  interior. 

The  errors  are  measured  again  and  the  process  is  repeated  until  a  con- 
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figuration  having  acceptable  error  characteristics  is  found.  It  is 
obvious  that  this  procedure  requires  considerable  experimental  time 
spent  by  very  experienced  personnel.  All  this  adds  up  to  high  costs. 

In  1976  Loyet  [43]  presented  a  scheme  for  digitally  compensating  an 
angle  tracking  system  for  radome- induced  error.  The  technique  involves 
measuring  the  error  characteristics  of  the  system  only  once.  The 
measured  errors  are  then  stored  as  a  function  of  scan  angle.  In  flight 
the  radar  is  allowed  to  track  normally,  but  the  angle  sent  to  the  guid¬ 
ance  computer  is  the  actual  angle  of  the  antenna  minus  the  pointing 
error  for  that  angle  as  obtained  by  interpolating  the  stored  error 
values. 

There  are  several  advantages  of  this  technique  over  the  traditional 
method.  The  system  need  be  measured  only  once,  thus  range  time  is  saved. 
No  high-level  expertise  is  required,  as  for  patching  The  technique  is 
easy  to  implement  with  a  microprocessor,  which  adds  negligible  weight 
to  the  system  and  can  be  tucked  away  in  any  unused  space.  By  increasing 
the  number  of  stored  error  values,  the  pointing  error  can  be  reduced  to 
arbitrarily  low  levels;  while  some  error  always  remains  due  to  patching. 

The  remainder  of  this  section  presents  a  compensation  scheme  based 
on  the  error  calculations  of  the  previous  chapter.  An  angle  compensator 
will  be  easier  to  implement  if  little  error  is  already  present;  thus  a 
radome  design  having  low  initial  error  characteristics  should  be  selected 
as  a  starting  point.  For  the  example  to  b«;  considered  here,  the  tangent 
ogive  with  constant  wall  thickness  designed  for  75°  incidence  will  be 
used.  The  in-plane  pointing  error  for  this  radome  with  the  previously 
used  antenna  were  presented  in  Figure  3-38. 
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Compensation  for  pointing  error  is  accomplished  by  subtracting  the 
pointing  error  at  a  given  scan  angle  from  that  scan  angle  to  obtain  a 
new  estimate  of  the  target  direction.  The  new  estimate  is  given  by, 

a1  =  a  -  ae,  (4.1) 

where  a1  is  the  new  estimate  of  the  a  scan  angle  to  the  target,  a  is  the 
scan  angle  of  the  tracking  null,  and  ag  is  the  in-plane  pointing  error 
at  a. 

The  scan  angle  to  the  target  takes  on  a  continuous  range  of  values, 
but  the  errors  can  only  be  stored  for  a  finite  number  of  angles.  To 
obtain  the  error  for  an  arbitrary  scan  angle  some  sort  of  interpolation 
is  necessary.  Two  methods  suggested  by  Loyet  for  doing  this  are  1.  to 
fit  a  polynomial  to  the  nearest  data  points  and  2.  to  fit  an  analytic 
expression  to  the  entire  set  of  data.  The  approach  taken  here  is  to 
simply  perform  linear  interpolation  between  the  nearest  two  data  points 
to  obtain  a  value  for  ag.  A  scheme  using  second-order  interpolation 
of  the  three  nearest  data  points  was  found  to  offer  no  significant 
improvement. 

The  in-plane  pointing  error  curve  for  45°  polarization  was  selected 
for  the  reference  error  curve.  As  seen  in  Figure  3-38,  this  error  curve 
is  near  the  average  overall  polarization  angles.  Reference  data  was 
stored  for  6°  increments  beginning  at  a  =  0.  Figure  4-1  shows  the  error 
in  a'  obtained  by  (4.1)  as  outlined  above.  This  figure  should  be 
compared  with  the  original  error  curve.  Figure  3-38.  The  difference  in 
ordinate  scales  is  over  six. 


POLARIZATION  ANGLE 


In-Plane  Radar  Pointing  Error  After  Compensation 
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Again  error  curves  are  shown  for  five  polarization  directions  (y  in 
Figure  2-17).  The  error  for  45°  polarization  is  zero  every  6°,  and  the 
other  curves  are  generally  on  opposite  sides  of  the  45°  curve  depending 
on  whether  y  is  more  or  less  than  45°. 

The  new  pointing  errors  are  smaller  than  the  original  at  all  scan 
angles  except  near  zero  where  the  compensation  scheme  actually  increases 
the  error.  This  is  because  linear  interpolation  of  the  error  between 
0°  and  6°  gives  a  negative  error  curve,  while  the  actual  error  is  posi¬ 
tive.  A  better  interpolation  scheme  might  begin  at  a  =  4°  and  not 
apply  error  correction  to  angles  less  than  4°.  Details  such  as  this 
would  need  to  be  worked  out  for  an  actual  system.  This  example  is 
presented  only  to  illustrate  the  method. 

Another  benefit  of  the  error  compensator  is  that  the  residual 
error  curve  is  more  rapidly  varying  and  noise-like.  This  actually 
improves  the  system  performance.  Recall  that  one  of  the  primary  func¬ 
tions  of  the  radar  is  to  provide  angle  time-rate  of  change  to  the  guid¬ 
ance  computer.  In  the  tracking  scenario,  the  target  is  initially  far 
off  axis  and  moves  toward  the  axis  throughout  the  flight. 

When  the  slope  of  the  error  curve  is  positive  the  radar  will 
over  estimate  the  angle  rate  of  the  target.  Oscillatory  flight  charac¬ 
teristics  may  result  as  the  missile  trajectory  is  changed  to  intercept 
what  appears  to  be  a  target  with  a  high  angle  rate.  When  the  error  slope 
is  negative  the  target's  angle  rate  is  under  estimated.  This  results 
in  a  sluggish  missile  response  to  target  motion.  Thus,  the  uncorrected 
system  will  over  estimate  the  angle  rate  at  the  beginning  of  flight, 
while  near  the  end  of  flight  the  rate  will  be  under  estimated. 
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Kuehne  and  Yost  [44]  have  shown  through  computer  simulation  of 
missile  flights  that  greatest  accuracy  is  achieved  when  the  radar 
pointing  error  changes  slope  frequently.  Under  these  conditions, 
periods  of  target-angle-rate  over  estimation  are  rapidly  followed  by 
periods  of  under  estimation  and  vice-versa.  Filtering  by  the  guidance 
loop  and  missile  dynamics  can  reduce  the  effects  of  this  type  error. 

Thus  the  error  tends  to  average  to  zero  Instead  of  having  effects  that 
accumulate  with  time. 

The  simple  error  correction  technique  presented  here  has  great 
promise  and  should  receive  further  attention  in  both  system  simulation 
studies  and  experimental  verification. 

This  section  will  conclude  with  some  comments  on  cross-plane  errors. 
Cross-plane  errors  result  from  asymmetry  in  polarization  and  aperture 
illumination  function.  After  the  above  discussion  of  in-plane  error 
correction,  the  next  question  is  logically--can  the  same  thing  be  done 
for  cross-plane  errors? 

Recall  from  the  cross-plane  error  curves.  Figures  3-40  and  3-41, 
that  at  a  given  a-scan  angle  the  error  is  a  function  of  polarization 
angle  (y  in  Figure  3-17).  The  cross-plane  error  is  a  zero-mean  function 
of  y  at  any  fixed  scan  angle.  It  can  be  shown  that  the  cross-plane 
error  curve  versus  y  has  half-wave  symmetry  with  period,  n,  that  is, 

Sg ( y )  ■  -Se(y  +  n/2),  (4.2) 

where  sg  is  the  cross-plane  error,  y  is  the  polarization  angle,  and  a 
is  fixed.  This  is  for  a  linear-polarized  incident  field  and  a  circular- 
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polarized  receiving  antenna.  For  a  1  inear-polarized  receiving  antenna 
the  cross-plane  error  would  be  odd. 

If  the  polarization  of  the  received  wave  were  linear  with  a  rapidly 
changing  radome  direction  with  zero  mean,  the  cross-plane  error  indicated 
by  the  radar  would  be  zero.  Also  if  the  received  wave  were  same-sense 
circular  polarized  there  would  be  no  cross-plane  error,  as  seen  before. 

Further  consideration  of  the  cross-plane  error  problem  must  be  coupled 
with  some  knowledge  of  the  scattering  properties  of  the  target. 

As  a  last  comment,  it  is  noted  that  a  very  "smart"  radar  system 
could  resolve  the  polarization  of  the  incident  field  by  use  of  a  dual 
polarized  antenna  and  necessary  signal  processing.  Having  determined 
this,  error  correction  could  be  made  as  outlined  before. 

COMPLEX  INDICATED  ANGLE  AND  WAVEFRONT  DISTORTION 

In  this  section,  a  somewhat  more  general  measure  of  radar  pointing 
angle  is  discussed.  This  measure  is  known  as  the  complex  indicated  angle 
(CA).  In  Chapter  II,  the  incidence  angle  of  a  plane  wave  impinging  on  a 
radar  antenna  was  determined  from  the  ratio  of  the  sum  and  difference 
output  voltages  by 

0  *  K'1  Re{^-}  ,  (4.3) 

where  9  is  the  incidence  angle  measured  from  boresight,  K  is  the  slope  of 
the  normal izec  difference  pattern  at  boresijht,  and  a  and  z  are  the  dif¬ 
ference  and  sum  voltages  respectively,  {see  (2.40)} 
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The  real  operator  (Re)  In  (4.3)  is  not  necessary  when  a  plane  wave 
Is  incident  on  the  antenna,  since  it  was  shown  in  Chapter  II  that  a  and 
Z  are  in  phase  quadrature  in  this  case.  However,  when  the  incident  field 
is  not  a  plane  wave,  the  phase  of  a  and  e  generally  differs  by  some¬ 
thing  other  than  90°.  In  this  case,  the  angle  indicated  by  the  radar 
may  be  considered  as  the  complex  quantity, 

CA  «  K'1^]  *  (4.4) 

Since  a  complex  indicated  angle  is  due  to  a  non-planar  wave  at  the  antenna, 
the  term  "direction  to  the  target"  may  have  no  meaning.  However,  the  real 
and  imaginary  parts  of  the  CA  are  due  to  specific  attributes  of  the  in¬ 
cident  wave  and,  therefore,  do  provide  information  about  the  incident 
field. 

Sherman  [45]  was  the  first  to  study  the  complex  indicated  angle. 

This  work  considered  a  distorted  wave  at  the  antenna  due  to  the  inter¬ 
ference  of  two  far-field  targets  with  different  amplitude,  phase  and 
direction.  He  showed  that  by  processing  the  CA  from  several  radar  re¬ 
turns,  the  direction  to  each  of  the  two  targets  could  be  determined 
although  the  targets  were  unresolvable  with  an  ordinary  monopulse  radar. 
Sherman  [46]  further  extended  the  theory  of  the  CA  to  include  the  effects 
of  multipath  providing  the  second  target  return  signal. 

Peebles  and  Goldman  [47]  also  considered  the  CA  method  for  targets 
with  multipath  by  examining  the  effects  of  the  terrain  surface  reflection 
coefficient  and  system  noise.  They  showed  that  the  CA  method  provides 
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significant  improvement  in  target  location  when  the  signal-to-noise 
ratio  is  greater  than  20dB  and  the  surface  is  almost  planar.  Sub¬ 
sequent  experimental  work  [48]  indicated  that  the  assumption  of  a 
planar  terrain  was  a  severe  limitation  to  the  CA  method  for  multipath 
because  the  roughness  of  the  actual  terrain  led  to  diffuse  scattering 
which  caused  ambiguities  in  the  locations  of  the  target  and  image. 

To  see  how  the  complex  angle  method  works  for  two  targets  or  a 
target  and  its  image  over  a  planar  terrain,  consider  two  unequal  am¬ 
plitude  plane  waves  incident  on  a  monopulse  antenna  near  boresight  at 
angles  0-|  and  0^  as  measured  from  boresight.  By  the  principle  of  super 
position,  the  total  sum  and  difference  output  voltages  are  the  super¬ 
position  of  the  voltages  produced  by  the  two  plane  waves  acting  alone. 
Thus,  the  complex  angle  is 


CA  =  K 


-1 


~J^1  ”3^2 


Z1  +  Z2 


(4.5) 


where  subscripts  identify  the  source  of  the  signals. 

Since  the  signals  arrive  from  near  boresight,  the  difference 
voltages  may  be  written  in  terms  of  the  sum  voltages  and  the  angles  of 
arrival  by  the  use  of  (2.40). 

This  yields. 


CA 


0 0 2^2 


Z1  +  Z2 


(4.6) 
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Dividing  numerator  and  denominator  of  (4.6)  by  Z-|  simplifies  the  ex¬ 
pression  for  CA  to 

,fe,  +  e,a  exp  (j<j>) 

CA  =  K'l— - - -  (4.7) 

1  +  a  exp  (j4>) 

where  has  ^een  exPressed  as  an  exP  (ji>)-  This  the  general  form 

for  the  complex  angle  produced  by  two  plane  waves. 

Several  things  are  to  be  noted  about  CA.  First,  if  signal  2  van¬ 
ishes,  (a=0)  then  the  radar  indicates  the  direction  to  the  source  of  1, 

9-| .  Likewise  if  1  vanishes,  the  indicated  angle  is  in  the  direction  of 
source  2.  However,  if  both  signals  are  present,  the  CA  is  truly  complex 
and  CA  depends  on  the  ratio  of  amplitudes  and  the  phase  difference  as  well 
as  the  angles  of  arrival.  The  real  part  of  CA  is  no  longer  an  indicator 
of  the  direction  of  arrival  of  either  wave  and  in  fact,  it  may  indicate 
the  "source"  to  be  well  outside  the  angular  region  subtended  by  the  two 
targets,  a  fact  confirmed  by  glint  and  multipath  studies. 

The  complex  angle  technique  for  resolution  is  based  on  the  ob¬ 
servation  that  the  CA  (4.7)  is  a  bilinear  transformation  [49]  of  exp  (j4>). 
Hence,  as  exp  (j$)  traces  a  circle  in  the  complex  plane  as  <j>  varies 
through  360°,  the  CA  also  traces  a  circle  in  the  complex  plane.  The 
radius  of  the  CA  circle  is  determined  by  the  constant  a. 

The  method  proposed  by  Sherman  was  to  display  the  CA  on  a  properly 
calibrated  oscilloscope  screen.  Since  the  relative  phase  of  the  two 
target  returns,  <j>,  will  be  rapidly  varying  due  to  the  motion  of  the  tar¬ 
gets,  the  spot  on  the  screen  would  trace  a 

circle.  The  constant  a  is  then  determined  from  the  circle  radius. 
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Equation  (4.7)  could  then  be  solved  for  0^  and  by  decomposing  the 
complex  equation  into  two  real  equations. 

Although  the  above  approach  gives  the  mechanics  of  the  complex 
angle  technique,  it  is  only  a  signal  analysis  approach  and  gives  little 
insight  into  the  structure  of  the  fields  that  are  incident  on  the  aper¬ 
ture.  Radome- induced  wavefront  distortion  is  the  source  of  a  complex 
indicated  angle  which  is  far  more  complicated  than  the  two  plane  wave 
case.  In  order  to  understand  how  this  arises,  it  is  necessary  to  begin 
at  the  aperture. 

To  illustrate  how  the  CA  arises  for  radome  distortion,  the  aperture 
integrations  used  in  Chapter  III  will  be  reduced  to  one  dimension  and 
the  incident  field  will  be  assumed  to  have  only  one  vector  component. 

With  these  restrictions,  the  aperture  integration  (3.15)  reduces  to 

w 

j  f(x)g(x)dx,  (4.8) 

-w  J 

w 

A  -  /  f(x)gd(x)d\,  (4.9) 

-w-' 


where  f(x)  is  the  incident  field;  g(x)  is  the  sum  illumination  function, 
gd(x)  is  the  difference  illumination  function,  and  2w  is  the  aperture 
width.  Recall  that  g  must  be  even  and  gd  must  be  odd.  Also  g  and  gd  are 
assumed  to  be  real . 

The  incident  field,  f(x),  is  a  comple:  function  and  can  be  decom- 


(4.10) 
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where  subscripts  e  and  0  are  for  even  and  odd,  subscripts  r  and  i  indicate 
real  and  imaginary,  and  each  f  is  real.  The  sum  and  difference  signals 
can  be  written  as 


z  =«  rr  +  j z. 

a  =  ar  +  ja. 


(4.11) 

(4.12) 


where 


A. 


'r  ~f  rre  g  d*>  (4.13) 

^  ~J f1e  g  dx,  (4.14) 

,  =/fro  9d  (4'15) 

Ai  s/f io  9d  dx*  (4*16) 

the  limits  of  integration  are  as  for  (4.8)  and  full  advantage  has  been 
taken  of  the  properties  of  even  and  odd  functions  integrated  over  sym¬ 
metrical  limits.  The  complex  angle  is  now  expressed  by 


CA 


.-1  V^r 

V^i 


(4.17) 


which  can  be  rewritten  as 


„  .  K-1  *  frVl 


(4.18) 


+  Z? 

r  i 


Recall  that  in  the  case  of  the  single  off-axis  plane  wave  the  even  and  odd 
components  am  in  phase  quadrature  and  (4. ’8)  predicts  a  purely  real  CA. 

Examination  of  the  imaginary  part  of  (4.18)  reveals  that  CA  becomes 
complex  when  the  sum  and  difference  voltages  are  no  longer  in  phase  qua- 
rature.  This  may  arise  in  the  case  of  two  or  more  plane  waves  incident  on 
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the  aperture  or  because  of  radome  distortion.  One  of  the  goals  of  this 
study  has  been  to  determine  if  the  imaginary  part  of  the  CA  would  pro¬ 
vide  additional  information  that  could  be  used  to  reduce  radar  pointing 
error  due  to  radomes.  After  examination  of  the  simulation  results  of 
the  previous  chapter,  it  was  concluded  that  the  imaginary  part  of  CA  could 
not  be  uniquely  related  to  radar  pointing  error  but  appeared  to  correlate 
with  phase  front  curvature.  Thus,  this  quantity  is  not  considered  useful 
for  pointing  error  correction. 


V.  CONCLUSION 


A  high-frequency  analysis  of  the  tangent  ogive  radome  has  been 
presented.  Geometric  optics  has  been  used  to  account  for  propagation 
of  the  electromagnetic  field  inside  the  radome  that  is  due  to  a  plane 
wave  incident  on  the  exterior.  Transmission  and  reflection  of  the 
field  at  the  radome  surface  was  obtained  by  treating  the  wall  and  the 
wave  as  locally  planar  and  applying  plane-wave  transmission  and  reflec¬ 
tion  coefficients.  The  principles  of  monopulse  radar  were  also  re¬ 
viewed. 

It  was  shown  that  radar  pointing  error  results  from  distortion  in 
the  incident  field  at  the  antenna.  Such  distortion  is  caused  by  the 
transmission  properties  of  the  wall  and  by  reflections  from  the  wall 
interior.  Both  these  distortion  mechanisms  were  seen  to  be  functions 
of  incident  angle  and  polarization.  Pointing  errors  were  found  for  both 
scan  planes  of  the  antenna  aperture.  It  was  seen  that  asymmetry  of  the 
incident  polarization  leads  to  cross-plane  errors  while  in-plane  errors 
exist  for  any  polarization  state. 

Complex  angle  processing  for  the  purpose  of  pointing  angle  reduc¬ 
tion  was  investigated  based  upon  the  simulation  described  above.  It 
was  concluded  that  the  complex  angle  could  not  be  uniquely  related  to 
radar  pointirg  error,  but  rather  appeared  ;o  correlate  with  phase  front 
curvature.  This  quantity  is  not  considered  useful  for  pointing  error 
correction. 
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A  simple,  but  effective,  method  to  compensate  for  pointing  error 
was  demonstrated.  This  technique  could  be  applied  to  any  radome  with 
good  results. 

This  analysis  has  considered  some  of  the  fundamental  factors  that 
give  insight  into  radome  performance.  Evaluation  of  the  results  pre¬ 
sented  can  only  be  based  on  experiments  that  both  validate  and  indicate 
areas  for  improvement.  Such  experiments  would  certainly  include  com¬ 
parisons  of  this  work  with  measured  pointing  errors  and  with  plots  of 
the  fields  inside  the  radome.  A  very  interesting  area  for  future  work, 
both  analytical  and  experimental,  is  the  effect  of  the  radome  tip.  Tip 
effects  may  well  account  for  many  differences  in  measured  and  theoreti¬ 
cal  results. 


* 
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APPENDIX 


This  appendix  contains  a  complete  listing  of  the  FORTRAN  program 
used  to  compute  radome- induced  radar  pointing  error.  The  program  is 
extensively  documented  and  even  a  cursory  reading  of  the  conment  cards 
reveals  much  about  the  computations  that  are  performed. 

The  program  requires  approximately  three  seconds  for  an  IBM  370, 
Model  155  to  execute  one  aperture  integration  over  an  8Xq  diameter 
aperture  with  *sXQ  spaced  rays,  including  reflected  rays.  Subroutine 
RTMI,  called  by  subroutine  INTRSC,  is  contained  in  the  IBM  Scientific 
Subroutine  Package  (SSP)  [39]  and,  therefore,  is  not  contained  in  this 
listing. 

The  surface  plots  presented  in  Chapter  III  were  generated  by  a 
modified  version  of  a  plotting  program  presented  by  Watkins  [50]. 
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MAIN  1 


THIS  PROGRAM  COMPUTES  RADAR  POINTING  ERROR 
FOR  A  CIRCULAR  POLARIZED  ANTENNA  UNDER  A 
TANGENT  OGIVF  RADGME 

THE  INCIDENT  FIFLO  IS  LINEAR  POLARIZED 
RADAR  POINTING  ERROR  IS  CALCULATED  FDR  VARYING 
ANGLES  C<=  INCIDENCE  AND  POLARIZATION  DIRECTION 

RE*  L  ARG(90) ,SRRRCA(90) , ERRRCB(90) 

REAL  CAAMPll 90) ,CAAMP2( 90)  .CAPHS1  ISO  )  ,CAPHS2(9!)  ) 
COMPLEX  EX!  ,£YI,EZI,EPEX,EPEY,EPEZ,FPAX,:pAY,E-PAZ 
COMPLEX  TPFR,TPAR.RPER.RPAR,CEXP,CMPLX,J/ (0., 1.  >/ 
complex  SUM* , SUMS, D!CA,DIF3 
COMPLEX  EXT  , F YT . FZT  »  FTTM  ,F  XR YR ,  EZR  ,  Z RTM 
COMPLEX  PHASE  .cXITZ.ZYITE,  iZm  ,E  X  ITM  ,E  Y  ITM,  £ZITM 
COMPLEX  PMAGlt  cMAG2fPHAS0R  ,CAA  ,CA3 
COMPLEX  FPANML . EPANX , EPANY , SFAMZ , = RGHT 
REAL  NX.NY.NZ.MONOSFt  3)  , Pi /. 31 41 593E1 / 

INTEGER  LABEU  2  )  /  '  ALPH  • «  'BETA'/ 

eGUI  VALENCE  (  X  P ,  X,M)  ,  1  ZP .  ZM  ) ,  t  COS  ANG,  A  Z  )  ,  (  S  I  NANG,  AY  ) 

COMMON/ RDAT A/  AL2, C , XM 2 . TTHET A , R AD ,A RC 
COMMON/MESFQS/  XM.YM.ZM 
COMMON  /WEIGHT/  APRAD  .IMS. RHOAP 
CDMM3N/TH/THIK 

PHAS  OP.  (  X  )=CMPLX(C  051  X  ) .  S  IN  (  X  )  ) 

OTR  =  P I / 1  BO  . 

RTD*13D./PI 
SORT  2=S0RT( 2.0) 

INITIALIZE  TRANSMISSION  COEFFICIENT  ARRAYS 
CALL  INIT131) 

ALL  DIMENSIONS  IN  MAIN  ARE  NORMALIZED  TO  WAVELENGTH 
BETA K=2.*PI 
C 

C  RADGME  PARAMETERS  (OUTER  S'JEFACc) 

C  RAD  R \0  0MC  RADIUS  AT  BA  Sc 

C  W  P  A  DP  ME  DIAMETER  AT  3AS- 

C  AH C CJUT  OAOOMr  GENERATING  ARC 

C  FINC  SS  RAC  3MZ  HEIGHT/OI A  METER 

C  YG  DISTANCE  OF  GI^BAL  POINT  FRMM  SASt  PLANE 

u 

C  RCaD  RAOOVL'  AND  Afs-RNNA  PARAVET'PS 

READ  (i,rr.  )  R  AO,  <=IUTSS  .  YG.APOIA 
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201  FORMAT!  4 F\ 0.0 I 
RADI N-RAO-THIK 
RA01N2=R4D1N*R AOIM 
W=2  • *P  AO 

ARCOUT*!  =IN=SS*FINrSS*0.251*W 
ARC I N* ARCOUT-T  HIK 
C  »ARC  OUT-P.A  0 

HIGH«S3r°T  (  (  2.  A  ARC  OUT-RAD  )*RAD) 

WRITE  RADCMfe  PARAMETERS 
WR I TS ( 6  »  10  5  1  W.HIGH.FINESS 
105  FORMAT! *  '  « • RAOQME  DATA',/.*  3A S 2= * . F7 .4 , / . 

#•  HE  IGHT  *'  . F7»4.  /»  '  FINENESS  R  ATI  0  **  .  F6»2  * // 1 

READ  COMPUTATION  OPTIONS 

READ (5.202 ) NANG. IMS. I RpP OS  .MAXPNT . NANGc » IP ;FL  * 

#OFL A  NG. CFLZ. SRRMAX 
202  FORMAT! I  2. 51 1 . 3X, 3F10.0) 

APRA  D=APQIA /2. 

06  LX  =0E  L  Z 

N'JM=LZ=APC!  A/0  "LZ*  1.5 

NUMF  L  X=  N'JMf-L  Z 

F I c  S  TZ*r'cL  l*  <N  HFLZ/2  I 

MONOS  FI 1 )- 1 .02  8/ (  1.61 7*A °0 1 A ) 

MOMCSc (2 1=1.2965/ ! 1,911 7 -APOI A ) 

ERRV,  AX*FRRMAX*DTR 

WRITH6»13  3)YG»APDIA»IMS*IRFP0$»IRT:FL 

103  FOP  M  <T  (  '  •  ,  •  A.f.r*ilNA  GIMBAL  P°INT  ='.p4.2,’  LAMDA»»/ 
#i  APORTURr  01  AMtT  Ek  ■«.F+.2,'  L  AMD  A  *  »  /  •  *  !MS«',I2. 
#/,«  I  RE  POS  =  '.12./.*  I  RE  FL  =  ’ , 12 ) 

*  IF  IRrPOS=l  ANTENNA  IS  REPOSITIONED  UNTIL  REA.L(CAI< 

1  ERRMAX.  *.ND  ERROR  IS  GIM3AL  ANGLE-INCIDENT  ANGLE 

Z  IF  IRrPOS=0  A  NTcNN A  IS  NOT  °EPOSI T IGNEC  ANO  ERROR  IS 

C  REAL ( CA )  WHEN  GIM  9AL  ANGLE  IS  AT  INCIDENT  ANGLE 

C 

DANGE=90.0/(NANGE-1»«DTP 

c 

00  700  !  TETM»1  ,NA  IGfc 
ANGS  =  <  ITt:TM-l )  *0A NGE 
SIN'A  N"  =  S  IN  ( A  *-'32 ) 

CO«B  r'-sCOS  (  ANG  El 
DO  400  I AN G= l ♦ NANG 

C  COMPUTE  ANGLE  JF  INCIDENT  «AVE c  RuNT 

*  1G! »Q ) .-( !/ NG-  1 1 "DEL ANG 
WRI  Tr  (6 . 1  C  4 )  A  NG  I 


%*%£*»**»*»  — •• 
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MAIN  3 

104  FORMAT!*  ' ♦•INCIDENT  ANGLE , F6 .2  ,• DEGREE S • ) 

APG(  IANG)  =  <?0.0-ANGI 
ANGI *4NGI*0TR 
TTHE  TA*  T  AM  (  ANG I  ) 

COMPUTE  DIRECTION  COSINES  DF  INCIDENT  RAYS 
OUTWARD  DIRECTION  ASSUMED  POSITIVE 
AX*0.0 
A  Y=S I  N(  A  NG I  ) 

AZ*COS(  ANGI  J 

FIND  DIRECTION  OF  UNIT  AMPLITUDE  INCIDENT  E-eI ELD 
TE  to  y  CASE:  EX  I  =  1. 0 »  3YI=0.0t  EZI=0.0 

TM  TO  Y  CASE:  EXI*0.0.  E YI =  *C OS ( ANGI ) ,  F ZI =-S I N ( ANG II 
EXITEO® 1 *0 

B  YI  TMO*CDS  (  ANG  I ) 

£Z  ITM0=-SI N ( A  NGI ) 

SET  INITIAL  POINTING  ANGLE  TO  INCIDENT  ANGLE 
ALPHA=90.0*DTR -ANGI 
RET t  *0.  0 

DO  BOD  I D0 I NT=  1  • MA  XPNT 

CDSA  LF  =  CHS ( ALPHA ) 

S I N A  L  p  =  S  I N  (  AL  P  HA  ) 

COSRFT=COS( BETA! 

SIN8ET*$IN(BETA) 

TANALFsTAN (ALPHA) 

SUMA  AND  C I  FA  ARE  MJNOPULSH  SIGNALS  IN  THE 
I N- PLANE  DIRFCTION  (ALPHA  SCAN  ANGLE) 

SUMS  AND  CIFB  ARE  MONOPULSE  SIGNALS  IN  THE 
CRDSS-PLAfF  DIRECTION  (SETA  SCAN  ANGLE) 

ZPR3  SUM  AND  DIF 
SUMA  *<0..0.) 

SUMS  =(0.,0.) 

DIFA  = (0 . . 0 .  ) 

DIPS  =<0.,0.) 

A  RC*  A  RCCUT 
AL2=*RC*#RC 

PCFF^P,M  APERT  JR  i  I VTIGRATI  ON  DF  DIRECT  FJLLD 


u'j  \L  I  n  FRACI1DA3LE 
iHlb  r'Aui  *  10BEC  _ _ 


o  n  o  o  o  n  no  non  on  no  on  on  no  non 
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MAIN  4 

00  50  IX*!  *  NUMELX 
XM0=FIRSTZ-< IX-1»*DELX 

(  XM ♦  YM«  Z  M)  ARE  RAOOMg  COORDINATES  OF  INTEGRATION  POINT 
(XM0.0..ZM0)  APERTURE  COORDINATES  OF  INTEGRATION  POINT 

00  150  L  I*  1  tN'JMSLZ 
ZM0*FIRSTZ-(LI— IJ^DELZ 

CHECK  FOR  MEASUREMENT  POINT  OUTSIDE  APERTURE 
RHOAP*SQRT ( ZM0* ZMD+ XM0* XM0  ) 

IF1RH04P.GT.APRAD)G0  TO  150 

ROTATE  ANTENNA  SAMPLE  POINTS  FCR  THE  GIMBALEO  ANTENNA 
XM=C  0S3ET*  XMO 

YM=COSALc*SINS  =T  *XM0-S  I  NAL  c*ZM0-*-YG 
ZM=S  IN\Lc*SINBET*XMO  +  CO$  ALc*ZM0 

FIND  INTERSECTION  OF  RAY  AND  R A  DOME 
ZR1  =  SQRT  (RAD*RAD-XM*XM  ) 

CALL  INTRSCI X, Y.Z.ZM.ZRl  ) 

FIND  SURFACE  NORMAL  AT  FAY-PAOCME  INTERSECTION 
CALL  N0RM( X,Y.Z,NX,MY,NZ> 

DETERMINE  INCIDENT  ANGLE  OF  f AY  AT  TRANSMISSION  POINT 
THETAI=ARCOSUY*NY  +  AZ*NZ  ) 

FIND  BINCRMaL 

8=(  N  X  A  )  / A E S (  N  X  A) 

CALL  CROSS  <NX,NY,NZ»0.0,  AY,AZ»BX,BY,3ZI 

COMPUTE  COMPLEX  INCIDENT  FIELDS  AT  ANTENNA 
PHAS  ==PH.AS0R{8ETA<*(ZM*AZt-YM’-AY  J  ) 

EXlTE=sXlTrO*PHASE 
EYITM=-YITMO*PHASE 
t-'ZITM=E  ZIT  M0*  PHASE 

rXI *CXITF*CCS/ NE 
E YI*E  YI TM*  SI  NA  ME 
EZ!*CZITM*SINtNE 

FIND  PARALLEL  AND  PCRPENDI  CUl  AR  COMPONENTS 
Oc  INCIDENT  FIELD  Ar  TRANSMISSION  POINT 

G  1=3  X+r  Xl+BY*EYI+5Z*0ZI 
"  P  F  X  =  B  X  *  -:MI  GL 
_-0EY  =  3Y*:MAGl 


T'”r'  "• 
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MAIN  5 


SPEZ  *BZ*FMAG1 

EPAX=SXI-EPEX 
EPAY*EYI-EPSY 
£PAZ=EZ!-CPEZ 

FIND  TRANSMISSION  COEFFICIENTS 
CALL  TRAMST(TH5TAI.TP=R.TPAR) 

EXT*EPC  X*TPER  +  cPA  X*TPAR 
EYT»FP£Y*TPCR+EPAY*TPAR 
EZT=EPcZ*TPcR+EPAZ*TPAR 

FIND  tm  COMPONENT  OF  TRANSMITTED  RAY 
ETTM»COSANG*EYT-SINANG*EZT 

COMPUTE  RHC  COMPONENT  AT  ANTLNNA 

ER  GHT  =  (  EXT-J*FTTM»/S0RT2 

SUWA =SJMA*WSUM( ZMO )*ERGHT 
OIFA*DIFA*WDIF(ZMOI*ERGHT 
S  JMB*S  JM9+WSUMI  XM0)*FRGHT 
DTFB  =  0IF<HW0IF( XM0l**PGHT 

150  C1NTINUE 
50  CONTINUE 

WTLL  WALL  REFLECTIONS  BE  CONSIDERED? 

IF  ( IRiFL.NF. 1  )G0  2C00 

COMPUTATION  OF  REFLECTION  DATA  BEGINS  HER;: 

CHECK  eOR  OPPOSITE  WALL  TR ANSM I SS I  ON-RE FL SCI  ION 
TIPSH0*AT/N(HI GH/RAD) 

IFUNGI.ST.TIPSHDIGC  TO  2000 
C 

arc=<r:in 

AL  2=  ARC*  ARC 

ZFIPST=-IFIX(HIGH/TTH-TA/DELZ)*DELZ 
XFTBSt*-IFIX(RAO/DFLX  >*D£LX 
NUMXP3W/0<rLX  +  ‘  .5 
OfcLZ P=D" LZ /S I N  (  AN3I  » 

NUMZP*-Z=IE ST/DSLZP+0.5 
YM=0.0 
C 

DO  105)  IZR  =  I.NUMZP 

zp*z-i  =  sr+n  zR-n*orLZP 


this  page  is 
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MAIN  6 


C 


C 

c 

c 


c 

c 


c 

c 


DO  1050  I XR*1 ♦ M JMXP 
xp»xfirstmixr-ii*delx 


OFTFRMINE  IF  (XP.O.tZP)  IS  INSIDE  RADOME 
I p ( xP*XP+ZP*ZP.LC.RAD*RAD) GO  TO  1050 


DETERMINE  IF  RAY  CAN  INTERSECT  RADGME 
IF(XP*XP.GF.RA0IN2)G0  TO  1050 

IF ( SORT ( A.L2-IC  ♦ABSI  XP  ) )**2 ) .LE.-ZP-TTHETA ) GO  TO 


1050 


FIND  RAY  REFLECTION  POINT 
ZLl=-S0RT(RADIN2-XP*X?) 

CALL  INTRSCIXR.YR.ZP.ZLl. 0.0) 

FIND  OJY WARD  NORMAL  AT  REFLECTION  POINT 
CALL  NORM(XR.YR,ZR.NX.NY,NZI 


C 

c 


c 

c 

I* 

r 

C 

c 

c 


c 


c 

c 


c 

c 


FINC  UNIT  VECTOR  IN  DIRECTION  OF  REFLECTED  RAY 
AD0TN2=2.0M  AY*NY+AZ*NZ  ) 

CX*A00TN2»NX 
CY=A00TN2*NY-AY 
CZ*A P0TN2*NZ-/ Z 

DETERMINE  IF  ANGLE  OF  INCIDENCE  DF  REFLECTED  RAY  IS 
GR-ATM  THAN  90.  DEGREES  AT  >.NTENN.A  APCR'P  JRE 
DIRECTION  OF  INCIDENT  FI=LD  ASSUMED  FOR  ANTENNA  NORMAL 
IF(CY^AY+AZ*CZ.GT.0.0  )  GO  TD  '050 

FINC  REFLECTED  RAY-ANTENNA  APERTURE  PLANE  INTERSECTION 
DIST  IS  REFLECTION  PO INT-I NTS RSECT ION  POINT  DISTANCE 

DI ST-( YG-YR-ZR*TANALF  > / ( CY+CZ*TANALF) 

X«ND*CX*DIST+XR 
YSNC=CY*DI ST ♦Y R 
ZFND=CZ*OI ST+ZR 

DETcRMlNc  IF  I  NTE  F  SECT!  ON  IS  INSIDE  APERT  jRi 
RHOAP=SORT  {  XFN D*XEND+ ZEND* ZEND +<YEND-YG)  **2) 

I  r  (R  HC'AP.Gfc  •  APRAD  I  GO  TO  1050 

FIND  RAY  TRANSMISSION  POINT 
ZPX  »  -ZL’ 

CALL  INTRSCIXT .YT.ZT.O.O.ZRl) 

CEDIN'-  INCIDENT  FIELD  MAGNITUDES 
r  X I =  t l  •  0 . 3  •  0  >  *  COSANr 
CYI  =  CM°L  X  (  EYIT^O ,  ).  )  I -*3  INANE 
r  ZI  *CMPLX(  0  ZIT  MO  .0.0)  *  SI  XA  NE 


tins  -■'■-r.. 

KrJkrM  out  i 


_  .Av  L  XvJ  1-iX.X. 
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MAIN  7 


FIND  INCIDENCE  ANGLE  AT  TRANSMISSION  POINT 
CALL  NORM(XT,YT,ZT,ANXT,ANYT,ANZT) 

ANGTRN*ARCOS ( A Y*ANYT +AZ*  ANZT ) 

FIND  TRANSMISSION  COEFFICIENT  AT  TRANSMISSION  POINT 
CALL  TRANSTi ANGTRN.TPER,  TPAR) 

FINO  81NORMAL  AT  TRANSMISSION  POINT 
8T»( ANT  X  Al/ABSIANT  X  A) 

CALL  CROSS ( AN XT »ANYT ,ANZT,O.Q,AY»AZ,BXT,BYTtBZTI 

FIND  PARALLEL  AND  PERPENDICULAR  COMPONENTS 
OF  INCIDENT  FIELD  AT  TRANSMISSION  POINT 

EMAG  1=8XT*£  XI  -f8YT*£YI  *BZT*EZI 
SPOX«BXT-*EMAG?. 

EPCY*3YT*2MAG1 

EPEZ=8ZY*£MAG1 


£  PAX  =£X1 -e  P£  X 

epay*eyi-epfy 

EPAZ3HZI-£PiZ 

*2IGHT  TRANSMITTED  FIELDS  WITH  INSERTION  TRANSMISSION* 
COEFFICIENTS:  TF.ANSMITT  £  0  FIELDS  BT  COME  INCIDENT 
FIELDS  AT  REFLECTION  POINT 

tXI=CP*X*=TPF7F.  +  EPAXTTPAR 
FYI*£P=Y*TPER+=PA Y*TPAR 
SZI=EPEZ«TPER+ePAZ*TpAR 

FIND  INCIDENCE  ANGLE  AT  REFLECTION  POINT 
ANGP  FF=AP.C  OS (-NY*  AY-NZ*  A  Z) 

FIND  REFLECTION  COEFFICIENT  AT  RAY  REFLECTION  POINT 
CALL  TRANSMANG’C-.FPER.RPAR) 

COMPUTE  COMPLEX  INC ICENT  FIELDS  AT  REFLECTION  POINT 
PHASE  =PHASOF  ( BETAKM  ZR*  A  Z+YP.»  AY )  ) . 

e XI* FXl*PHf SE 
E  Y!  =7  SE 

?Z I  *  £Z I  *PHA  SF 


FINO  8T  NORMAL  AT  ELECTION  POINT 
3P.=  (-N  X  H  >  / A  B  S <  -  M  X  A  I 

CALL  CROSS  (-NX  ,-NY  ,-N 2,0.0  .AY.  AZ.3XR,  3YF.»  3ZR) 

THIS  PAGE  13  Bwl  QUALITY  PRACTICABLE 
VtijU  (JUPi  H  JtoUSliiSD  10  BDC  — — 


o  o  ooo  on  on  o  oo  on 


183 


MAIN  8 


FIND  PARALLEL  AND  PERPENDICULAR  COMPONENTS 
OF  INCIDENT  FIELD  AT  REFLECTION  POINT 

EMAG2*BXR*EXI+BYR*EYI*BZR*iZI 
FPEX»BXR*EMAG2 
FPEY«8YR*EMAG2 
EPEZ«BZP*EMAG2 

EPAX«5XI-FPFX 
EPAY*EYI-=PEY 
EPA  Z  *E  Z I-c  Pc  Z 

FIND  NORMAL  PA RT  OF  INCIDENT  FIELD  PARALLEL  COMPONENT 
EPANML»C  PAX*NX4-EPAY*NY+EPA  Z*NZ 
EPANX»EPANML*NX 
EPAMY=EPANML*MY 
EPANZ*?  PANML*NZ 

WEIGHT  REFLECTED  FIELDS  WITH  REFLECTION  COEFFICIENTS 
E XR=EP£X *RPER+R PA A* (EPAX-2 .0*=PANX» 
EYR*EP“Y',tPP£F+RPARM2PAY-2«0*EPANY) 

F  ZR  *E  PE  Z*R  PE  F  ♦  RPAR-  ( £  PA  Z-2  .0«r£PANZ  > 

COMPJTS  INCIDENT  FIELD  AT  APERTURE  BY  ACCOUNTING  FCR 

phase  shift  fpom  reclection  point 

PHAS  E=°HASOF C -BETAK*DI ST  ) 

^XRsCXR^PHASfc 
<!  YRs?YR*OHAS,I 
EZR*FZR*PHA  SF 

ERTM=F.YR*S  INALF-E  ZR*CQSALF 

FIND  PHC  COMPCNCNT  OF  REFLECTED  FIELD 
FRGHT*IEXR-J*FRTM)/SQRT2 
C 

C  COMPUTE  OBLIQUITY  FACTOR 

OBLI  Q=0.5-0.5* IAY*CY  +  AZ*CZ) 

ZSND*Z  END/CCSA  LF 
WGHT  SA  =  W  SUM  (ZEND  )  *0Bl  IQ 
WGHTDA* WOI F ( Zr  NO )*G6L I Q 
WGHTS?*WSUM( XEND )*CBLIQ 
WGHTD3=WDI Ft XEND) *OBL IQ 
C 

SUMi »<JM\+LRGHT*WGHTSA 
01  FA  sOi  F-.  +  Cr  GH  T*  wGHTv  \ 

RGH  T*wGHTS  3 
01  Ft)  =  rl  FB4,cRGH'r'«WGHTC  B 
1U50  C  jnti'.’u: 

THIS  PAGE  IS  :  *  '  v<rA'  PRACriw  r.j 
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C 

2000  CONTINUE 
C 

PRINT102.SUMA.SUMB,0IFA,CIFB 
102  FORMAT)'  • . 'SUMA*' # 2E 15. 6. 5X,' SUMB*' ,2615.6,/ 

A.'  0IFA*'.2£15.6.5X,'0IF8*',25I5.6) 

C 

C  COMPUTE  COMPLEX  INDICATED  ANGLES 
CAA*-MCNOSF )  I  MS )  *  J*DI  FA/SUMA 
CAB*  +MONOSF  ( IMS ) * J*CI FB/ SUMS 
C 

CAAMPA*C  ABS ( CAA)*RTD 
CAAMPB*CABS ( CA  3 ) *RTD 

CAPH SA* ATAN2 ( A IMAGIC AA » , RE AL (CAA ) I *RTD 
CAPHS8*ATAN2</  IMAG)CA3)  ,  RE  AL  I CA  B )  )  *RTD 
C 

WRIT  5  (6. 10  8)  C  AAMPA.CAPHSA  .CA  AMP6,  C/-PHSB 

108  FORMAT) •  '.'ALPHA  COMPLEX  ANGLE  *  ',!:12.4.'  AT',F6.1, 
A'  OF  GREFS'  »  /  *  '  BETA  COMPLEX  ANGLE  *',E12.4.'  AT', 
AF6.1.'  PCGPF£S  •  I 

C 

CAAMP1)  I  ANG)*C  AAMPA 
CAAMP2(IANG)*CaAMPS 
CAPHSK  IANG )*CAPHSA 
CAPHS2  (  I  1NGJ*CAPFS9 
ALPS  RR=R  :A  L  (CAA) 

BFTERR*RFAL(CAB» 

AOF G*ALF  :RR*RTO 
3DEG*9ETftRR*RTD 
PRINT  109 *A DEG, 3D EG 

109  FORMAT!*  ', 'ALPHA  INDICATED  ERF OR* '» F10  *6 ,' DEGREES '»/ » 
A*  BETA  INDICATED  ERROR* ', F 10.6 , '  0EGRE5S ' I 

C 

IF)IR5PCS)750, 750,751 
C 

751  ALPHA*ALPHA+ALF=RR 
9ETA*B"T  A+BTTb  RR 
ERF  OR  A  *  (  \LPHA-PI/2.0*ANGI l  *FT0 
F RF  OFB*B  *TA*RTC 

I F) ABS) ALFERRI .LT. ERRMAX.AND.ABSI35TERRI.lt. ERF  MAX) 

AGO  TO  350 
C 

300  CONTINUE 
GO  TO  350 
C 

750  PRR  ORA  *A  D’:G 
FRRnRe*2orr, 

3  50  PP  I NT10  6  »T  F  ROF  A .  F  RRORB  ,  I  PC  I  NT 


THIS  PAGE  IS  BEST  QUALITY  PRAr.TIP.il BLE 
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MAIN  10 

106  FORMAT! •  «,  ‘ALPHA  POINTING  ERROR*  *.F8.5,«  DEGRESS ' 
*/,*  SETA  POINTING  ERROR* *. F8,5,  •  OEGRESS,  OBTAINED 
*  ’AFTER  ' .1  3 » '  ITERATIONS* • //I 

C 

ERRRCA! IANGI*FRRORA 
ERRRCB! I ANG) *E  RRORB 
C 

400  CONTINUE 

ANGE*ANGS*RTD 

PRINT107.ANG5 

107  FORMAT! »  *.'E  ANGLE *•  . F6 .3 » •  DEGRESS' I 
C 

CALL  LPLOT ( ARG .ERRRCA*  CAPHS1 .NANG. LASEL !  1 1  » 1 1 
CALL  L°L  QT I ARG  »S  PRRC8 » CAPHS2  .NANG.LA8EL  (  2  I » 1 1 
CALL  LPLOT !APG,CAAMPl,CAPHSltNANG,LA3EL( II ,  l> 

CALL  LPLOT! ARG* CAAMP2 , CAPHS2 , NANG. LABEL  !  2 ) til 
700  CONTINUE 
C 

STOP 

END 


C 


N3PM  17 

S UBR  OUT  I  MS  NOP  M ( X , Y . Z , NX , N Y .  !U  I 

THIS  SUBROUTINE  COMPUTES  THE  OUTWARD  SURFACE 
REAL  NX. NY, NZ. NR 

COMMON/ RDATA/  AL2 , C , XM2 . TTHETA , RAD  , ARC 


XsX+l.'-ZO 

z*z+i.:-2o 

R*SQRT(X*X-t'Z*!Z  I 
NY* Y /ARC 

NR*S QRT( l.-NY*NYI 

NXsX/R^NR 

NZ*Z /R  *N R 

RETURN 

FND 


^SPAGIilq, 
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INTRSC  11 


SUBROUTINE  INTRSCt  X,  Y  ,  Z,  XL  1,  XR1 1 

THIS  SUBROUTINE  COMPUTES  Thfc  RAY-RADQME  INTERSECTION 
POINT,  (X,Y,Z>,  OF  A  PAY  PASSING  THROUGH  THE  POINT, 
(XM.YM.ZM),  AT  THE  SPHERICAL  ANGLE  THETA 
WITH  XL 1<Z <XR 1 . 


EXTERNAL  FCT 
OATA  EPS/2.0F.-4/ 

COMMON/ ROA  T  A/  AL2,C,XM2, TTHETA  *  RAD , ARC 
COMMON/ MTSP OS/  XM.YM.ZM 
XM2*XM*XM 
C 

C  CHFCK  FOR  NORMAL  INCICSNCE 

IFCABS(TTHETA) .GT.l.E03)GQ  TC  SO 
C 

C  FPS  IS  MAXIMUM  ERROR  IN  WAVELENGTHS 

CALL  PTMKROOT.F.eCT.XLI.XRl.EPS,  25.ISPI 

Z*ROOT 

X*XM 

Y=SQCT(AL2-<  SORTC  XM2+Z*Z  )-^C  1  =*-*2  > 

IFdFB.EQ.OlPFTUPN 

PRINT  20.IFR 

23  FQRMATC  •  ,  'SUBROUTINE  RTM I  DID  NOT  FIND  RAY*, 

* 'INTERSECTION  POINT*,/,*  £  HR  UR  CODE  IS  IiR=  ',12) 
STOP 
<50  X*XM 

Z=ZV 

Y=$ORT(  AL2-C  SORT  l XM2+Z*Z  H-C)-**t2  ) 

RETURN 
•  END 


FCT  12 


FUNCTION  FCT(Z) 

FCT  IS  THE  DIFFERENCE  IN  THE  Y-COQRC INATCS  0=  A 

PAY  PASSING  TH50UGH  (XM,YM,ZM)  AND  THL  MDOM£ 

SURFACE  AS  A  FUNCTION  OF  Z 

CO  AMON/ RCATr /  / L2 , C , XM2 , T?HI T A , RAD  , ARC 

COM M ON/ MCS  POS /  XM,YM,ZM 

eCT* YM>  (  Z- ZM )  *  TTHiT A— SORT  {  £L2-  (  S  3P.T  (XM2  +  Z*Z  )  +C  I  **2) 

RrTURN 

FNO 


lsB£S1(lu ALin  PRACUCABLE 

XHlSPA^lb^oI^ 
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SUBROUTINE  INIT(N) 

THIS  SUBROUTINE  INITIILIZGS  THE  MATRICES  USED  FOR 
TRANSMISSION  AND  REFLFCTUCN  COEFFICIENT  INTERPOLATION 
N  EQUISPACED  DATA  POINTS  ARE  USED 

COMPLEX  ATPERt  501 .ATPARI50I .ARPSRC50I .ARP4R150I 
REAL  ANGC50).P 1/3.141593/ 

COMPLEX  TPFR»T  PAP  »RPEP.  *RPAR»  PHAS3P. »CMPLX 
COMMON/TOATA/ATPER » ATPAR ♦ ARPER . ARPAR , ANG . 01 FO.CIFl , 
A0IF2.0CL.M, 00. 01,02 
PHAS  OR(X>*CMPLX(COStX) »  S  IN  ( X  )) 

DEL*PI/2.O/FL0AT(N-l) 

00  10  1*1. N 

ANG<  I  l=OEL*<  I-ll 

CALL  TRANS (ANG(II»TPER»TPAR,ARPeR( I ) » ARPAR II), 
RPERIPD.PARIPD) 

COMPUTE  INSERTION  TRANSMISSION  COEFFICIENTS 
ATPERd  )=CaBS( TPER)*PHASQR(“P5RI PD) 

ATPAR 1 1 ) =CABS <  TPAR )*PHA SQR (— PARI  PD ) 

10  CONTINUE 

00*2  .0  *L'GL*OEL 
Ol=-CEL*OGL 
0  2=0  0 
M=N 

RETURN 

SNO 


WSUM  22 

FUNCTION  WSUM(ZMO) 

SUM  WEIGHTING  FUNCTION 
PEAL  PI/3. 141593/ 

COMMON  /  /FI  GHT  /  APRAO.IMS.R 

IMS* 1  UNIFORM 
I  MS  *  2  COSINE  TAPER 

GO  T0( 10,20). IMS 
C 

10  wSUM=1.0 
RFTURN 
C 

20  WSUM*C3S(PI/2.*R/APRADI 
RCTURN 
cND 


SHIS  PAGE  iC  V.-' "7  rViACTI  . 
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WOIP  23 

FUNCTION  WOIF(ZMO) 

DIFFERENCE  WEIGHTING  FUNCTION 
COMMON  /WEIGHT /  APRAD.IMS.R 
GO  TOU  #2)  .IMS 

1  WDIF*ZM0*2.0/A  PR AO 
RETURN 

2  WDIFaZM0*2.93i06/APRAD*WSUM|ZM0> 
RETURN 

END 


TRANST  14 

SUBROUTINE  TP  AN$T( ARG  »TPEP  I  »TPAR!  ) 

THIS  SUBROUTINE  INTERPOLATES  THE  TRANSMISSION  CAT* 
SUPPLIED  BY  IN  IT  WITH  A  SECOND  CROER  CURVE  THROUGH 
THE  THR  S'*  CLOSEST  DATA  POINTS. 

COMPLEX  ATP5R1 50) ,ATPAR< 50 ) »ARPER( 50) t ARP AR ( 50 ) 
REAL  ANG( 50 ) «AL( 3  ) 

COMPLEX  TPRFI , TP  API .CMPLX, PHASQR 

COMv  CN/TCAT  A/AT PER .AT^AP , ARPER  » APP AR , ANG . DI FO » 

#01  FI  ,OIF?,OEL ,N,00  .01,02 

I *AP  G/DEL+  1.5 
IFU.cO.  1)1*2 
I «={  I  .EG.N) !  *N—  1 
OIFO*APG-ANG(  I-?,  ) 

OIF l«ARG*-ANG { I  ) 

0IF2*APG-ANGl I ♦!) 

ALU  )=OI  -1^0:F2/D0 
AU  2  )=0TF0*0!F2/D1 
ALI  3)=0IcO-^0IFl/02 
TPF°  I=(0.0  .0.0) 

T°AF  I  =  (  0  .0 , 0 .0  ) 

co  io 

TPnPIsTP';RI+iTPFRII-2*KI'‘,AL(K) 

10  TPAF I»TPAF I+ATOAP ( I-2+K )  *AL(  K. ) 

RETURN 

END 


.'.5*  13  best  quality  practicable 
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TRANSR  15 

SUBROUTINE  TRANSR ( A RG ,RP£R  . RPAR J 

THIS  SUBROUTINE  INTERPOLATES  THE  REFLECTION  DATA 
SUPPLIED  BY  INIT  WITH  A  SECCNO  ORDER  CURVE. 

COMPLEX  ARPERI 50) •  A  R  P  A  R  (  50  )  .RPSR.RPAR,ATP£R(50i . 
fATPAR ( 50 ) 

REAL  ANG(50  J » AL(3  ) 

COMMON/TQATA/ATPER .AT PAR  .ARP5R  »ARPAR » ANG  *  DI  FOt 
*OIF1.DIF2,0EL,N,DO»D1.D2 

I-ARG/OSL+1.5 
IFCI.EO. 1) 1*2 
IF  C I .FO.NII-N-1 
DIFO*ARG-ANG(  1-1) 

0!  FI  *ARG-Af'G(  I ) 

CIF2*ARG-ANG(  I *11 
AL ( 1 ) *0 1 Fl»D IF  2/DO 
AL(2»=CIF0*DIF2/01 
AL( 3 )*0I F0*DIF I/O 2 
RPER*<0 .0.0.0) 

RPAR=(Q. 0.0.0) 

DO  10  <*1.3 

RPFR*RP<=R+ARPER(  1-2  ♦<  MALI  Kl 
10  RPAR=RPAR+ARPAR<I-2+K>*AL( <> 

RETURN 

END 


CROSS  16 


SUBROUTINE  CROSS (AX.AY,AZ,3X,3Y,BZ,CX,CY,CZ) 

THIS  SUBPOUTINC  COMPUTES  THE  NORMALIZED  CROSS 
PRODUCT  OF  A  INTO  3 
C*U  x  B ) / ABS l a  X  8) 


CX*A Y*BZ-AZTBY 

CY*AZ*3X-AX*'B  Z 

CZ  *.\  XvBY-AY*3X 

CMA  ;  =  SOPT(CX*C  X*CY*CY+CZ*CZJ 

CX-C  x/C  MAG 

CY*C Y/CMAG 

C7*CZ/CMAG 


of TURN 
END 


'*•  r.^ 


TRANS  18 


SUSP  OUT  INF  TRA NS < THETA , TPE RP .TPARA , RPERP , RPARA , 

A PERI PD.PARIPO) 

C 

C  THIS  SUBROUTINE  COMPUTES  COMPLEX  TRANSMISSION  AND 

C  REFLECTION  COEFFICIENTS  FOR  MULTIPLE-LAYER, 

C  SANDWICHED.  PLANAR  SHEETS. 

C 

C  ALL  (DIMENSIONED  VARIABLES  HAVE  BEEN  DIMENSIONED  FOP 
C  A  MAXIMUM  OF  FIVE  LAYERS. 

C 

C  OIMFNSIONS  IN  THIS  SUBROUTINE  ARE  IN  METERS 

C 

C  TPSPP  COMPLEX  PERPENDICULAR  TRANSMISSION  COEFFICIENT 
C  TPAF.A  COMPLEX  PARALLEL  TRANSMISSION  COEFFICIENT 
C  RPORP  COMPLEX  PERPENDICULAR  REFLECTION  COEFFICIENT 
C  RPARA  COMPLEX  PARALLEL  REFLECTION  COEFFICIENT 
C  THST A  ANGLE  OF  INCIDENCE 

C  NL  NUMBER  OF  LAYERS  SANDWICHED  TOGETHER 
C  BETA  PROPAGATION  CONSTANT  IN  AIR 
C  MU  PFR1EABILITY  Q<=  LAYER  RELATIVE  TO  AIR 

C  EPSILN  DIELECTRIC  CONSTANT  OF  LAYER  RELATIVE  TO  -MR 
C  TVNOEP  DIELS CTFIC  LOSS  TANGE  NT 
C  SIGMA  CONDUCTANCE  OF  LAYEP. 

C  GA-1MAI  PROPAGATION  CONSTANT 
C  YIPERP  PERPENDICULAR  ADMITTANCE 
C  Y I  PARA  PARALLEL  ADMITTANCE 
C  THIK  THICKNESS  OF  LAYER 

C  T 0 T 0  TOTAL  THICKNESS  OF  SANDWICHED  LAYERS 

C  PERI  PD  PERPENDICULAR  INSERTION  PHASE  DELAY 
C  PARIPO  PAFALLf-L  INSERTION  PHASE  DELAY 
C 

COMPLEX  Z.CLXP.CSQRT.CDY.J/tO.O.l.Ol/.CMPLX 
COMPLEX  TPFRP  »  TP  AP.A  »KSP»  GAMM  AI  ,YI  PAR/  ,  YI  PijRP,  =  ,  F 
COMPLEX  aPP,3PP.C?P.DPP, APR, 3PR.CPR.DPR, APRPR.BPRPR 
COMPLEX  CPRPR.DPRPR.APfFR,  BP  AF.R »  CPARR.  DPARR 
COMPL-X  G,H,S.T,RP?RP,RPARA,CCOSH,CSlNH 
COMPLEX  APrRP  4  5 ) • 3PERP ( 5 » , CPFRP ( 5 ) .OPERP ( 5) 

COMPLEX  iPAPA  <  51  .BP'PA  (5).CPARA(5)  ,DPARA  (5) 

DTVENSICN  KMU( 5)  ,KEP<  5) , GAMMA! (5) , Yl PERP ( 5  I , Y I  PARA ( 5  ) 
DIMC  NS  I CN  FPS I LN ( 5 ) »LYM ( 5 ) .THIK ( 5 1 , MU ( 5 ) , S 1 3MA< 5) 
DIM-NSICN  TANO=P<  5)  ,3DY(  5) 

Pc  A  I.  LAM  DA  »  MU »  PI/3. 14  15«3/,K  VJ 
I  N’t  G-  0  LFLA.G/0/ 

SOUI  VALiMCr-iMU.KMU) 

CCMMCN/TH/THIKWL 

r 

CCOSH(Z)=(CcXP(Z)+CiXP(-Z)  >  / 2. 
f  SINHm  »(Cf  XPl  Z)-CCXP(-Z»  )/2. 


TtUS  PAni 
i'KuM  tor  i 
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Trans  19 


SINT*SIN(THFTA  ) 

COS” =CDS {T  H£T A  I 
SINT  2*SI NT*  SINT 

CHECK  IF  CAT A  IS  TO  B=  READ 
I  Ft  LFLAG.lt . 0) GO  TO  100 

500  00  22  I * NL 

GAMMAIt  I  )  =  J*35T/,*CSQRT<KPUtI  )*KEPt  D-SINT2  ) 
YIPERP  <  I  )®CSORT  t  KMU  ( I  )*KEP  (D-SINT2I /KMU  ( I  » 
YI  PARA  1 1  J=KFP(II/CSORT(K-«um*KEP{  II-SINT2) 
ODY (  II-GAMMAllI )*THlKtI  ) 

APERPtI ) =CCOS  H ( OD Y  C I ) ) 

BPEF  Ptl) sCSINHtOOYCI)  ) /YIPERPt I ) 

CPERPd  )=YIPERPtl  t*CSINH{DDYtI  )  ) 

DPERPII  )  =APERP  ( I ) 

APARAt I »=APF  RP  (I  ) 

BPARA ( ! ) =CS INHtOOYt I ) )/YIPARA(I) 

CPAP  A ( T 5  . 1 PARA  tI>*CS INHIDOYt I ) ) 

DPARAt  i  .  'CPERPtl  t 
22  CCNT  I  f  uE 
C 

IF  (  NL.SO.l IGQ  TQ  25 
C 

APP=APERPt  1) 

BPPsBP^RPt 1) 

CPPsCPERPd  ) 

QPP-OP=RPt  1} 

APR=APARA( 1 J 
8PR=  BPAR.lt  1  ) 

CPR*CPARA(  1 ) 

OPR=OPARA< 1) 

C 

HO  3  5  !=2,M 

APRPR*APP*APERP(I >  +SP P*C  PiRP  1 1  ) 

3  PR  P  R  =A  P  P*  B  P.=  R  P 1 1  )  B°  P*  D  P”  R  P  ( I  ) 
CPRPR=CPP*/iPERP(I  )  +0P  P*C  P.^RP  ( !  ) 
OPRPK=CPP^BPi:RPU  )  ♦  CPP*C?::RP(  I  ) 

APARF=APR*  APARAt  IH-BPR*C  PA  RA  (  I  ) 

8PA k R=APP*BP  *  R At  1 1 +8PR “C  PA  RA ( I ) 
CPA'.RsCPRVAPARA  ( I J  *QPR*  CPAFA  (I  ) 

OPA  RR=C  PR*  SPAR  At  I  I  ♦CPR^OPARA(I) 

APP  =  AP  R  PR 
BPP.-*ePR  P® 
rPP=CPRPR 
OPP^OPP  op 

APR-* APu°R  BU?  --Au/.  !>‘- 

iw«  ■  ■  •' 
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TRANS  20 

BPR  *  BPAR  R 
CPR=CPAP.P 
OPR«OPARR 
35  CONTINUE 
30  CONTINUE 
C 

E*APRPR  *8PR PR* COST  ♦CPRPR/CCST+DPRPR 
.  F=APARR+BPARR/C0ST4‘CPARR*CQST«-DPARR 
TPARA*( 2.0.Q.0I/F 
TPEP  P=(2.0,0.0)/S 
C 

G=APRPR+BPRPR*C0ST-CPPPR/C3St-DPRPR 
H*A  PR  PR  ♦9  PR  PR*  COST  +CP RPR /COST +QP RPP 
RPSR  P*G/H 

S=AP  ARR  +  8P ARR/COS T-CP  ARR  *C  GST-DPARR 
T=<*  PARR  ♦BPAR  R/COST+CP  AP.R*CCST  +-DPARR 
RPARA=S/T 
C 

PCCN1=8~TA*TOT0*COST 

PERIP0=-ATAN2(AIMAG(TPHRP)  ,R»iAL  (TPGRP  )  J-PCLN1 
PARI P0=-ATAN2  ( AIM  AG(TPARA) »R2AL ITPAfU ) )-PCONl 
C 

50  I F  {  P  AF  I  P  0,  GI*  «0  •  I  GO  TO  60 
PAPI  PC-PiPIP0  +  2«*PI 
GO  T  0  50 
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